Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,608,4950GT53.6% 3.2 / 17.7 28intergenic (+18/‑211)tenA/codBAminopyrimidine aminohydrolase/Cytosine permease
Reads supporting (aligned to +/- strand):  ref base G (0/13);  new base T (0/15);  total (0/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AACTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAATGGGGTTGGACGG‑GACCGTGTAAA‑CCGCCGTGGGCTGGGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGCGAAC  >  ACB143/3608362‑3608616
                                                                                                                                       |                                                                                                                         
aaCTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAATGGGGTTGGACGG‑GACCGTGTAAA‑CCGCCGTGGGCTGGGTGa                                                                                                                          <  6:30679/135‑1 (MQ=255)
aaCTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAATGGGGTTGGACGG‑GACCGTGTAAA‑CCGCCGTGGGCTGGGTGa                                                                                                                          <  6:27933/135‑1 (MQ=255)
    ggCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGTATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAATGGGGTTGGACGG‑GACCGTGTAAA‑CCGCCGTGGGCTGGGTGATtact                                                                                                                     <  7:281433/136‑1 (MQ=255)
                    aCCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGCCCAGGCCTCCCGCCGGGAGGGGGATTTCTGGCAATTGGGGTTGGACGG‑GACCGTGTAAA‑CCGCCGTGGGCTGGGTGATTACAACCCAGCCCAGGACCg                                                                                                     <  2:224924/136‑1 (MQ=255)
                                             cAGGGATTTTTGGCCAGGCCTCCTGCCTGGTGGCGGATTTCTGGCAATGGGGGTTGGCCGG‑GCCCGTGTAAA‑CCGCCGTGGGCTGGTTTATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAAt                                                                            <  1:247396/136‑1 (MQ=255)
                                                       gcgaccggCCTCCCCCGTGGTGGCGGATTTCTGTCAATTGGGGTGGGAGGC‑GAGCGTCAAAC‑CCCCCGGGGGCTGGGTTATTACTACCCACCCCAGGACCGCTTTGGTGCATGGGCGCCCCAAAATAAGTCGTccc                                                                  <  6:110600/130‑1 (MQ=255)
                                                        gccccggCCTCCAGCCTGGTGGGGTATTTCTGGAAATTGGGGTGGGCCGG‑GCCCGTTAAAC‑CCGCCGTGGGCGGGGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTccc                                                                  <  7:210605/130‑1 (MQ=255)
                                                                            tggggttttCCGGCAAATGGGGGTGGGCGG‑GCCCCGGTAAACCCGCCGTGGGCGGGGTTATAACAACCCCCCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGa                                              <  5:147025/130‑1 (MQ=255)
                                                                            ggggtttttCTGGAAATTGGGGTGGCCCGG‑GACCGTGAAAA‑CCGCCGGGGGCGGGGTGATTACAACCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGa                                              <  1:120688/129‑1 (MQ=255)
                                                                                gATTTCTGGCAAATGGGGTTGGACGG‑GACCGTGTAAA‑CCGCCGTGGGCTGAGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCa                                         <  1:62752/136‑1 (MQ=255)
                                                                                          atttGGGTTTGCCCGG‑GACCGTGTAAA‑CGGCCGTGGGCTGGGTTATAACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                               <  7:74733/133‑1 (MQ=255)
                                                                                          aattGGGGTTGCACGG‑GACCGTGTAAA‑CCGCGGTGGGCTGGGTTATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                               <  8:187701/133‑1 (MQ=255)
                                                                                          aatgggggTGGGAGGG‑TCCCGGTTAAC‑CCGCGGGGGGCGGGTTTATACCAACCCACCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                               <  1:93634/132‑1 (MQ=255)
                                                                                           attGGGGTTGGCCGG‑GACGGTTTAAA‑CCGCCGTGGGCTGGGTTATAACTCCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgtg                              <  5:141679/134‑1 (MQ=255)
                                                                                           attGGGGTGGGACGG‑GCCCGTTAAAC‑CCGCGGTGGGGGGGGTGATTACTACCCCCCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgtg                              <  8:342083/134‑1 (MQ=255)
                                                                                           actGGGGTTGGCCGG‑GCCCGTGTAAC‑CCGCCGGGGGCGGGGTTATAACTCCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgtg                              <  8:350452/134‑1 (MQ=255)
                                                                                           acgggggTTGGCGGG‑GCCCGTGAAAC‑CCGCGGTGGGCGGGGTTATAACTCCCCCCCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                               <  3:320020/132‑1 (MQ=255)
                                                                                           aagggggTTGGACGG‑GACGTTGTAAC‑CCGCGGGGGGCGGGGTGATTACAACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgt                               <  3:270110/132‑1 (MQ=255)
                                                                                           aaTGGGGTGGGCCGG‑GACGGTGAAAC‑CCGCGGTGGGCTGGGTGATTACTACCCTGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgtg                              <  1:119292/136‑1 (MQ=255)
                                                                                            tttggggTGGGCGG‑GGCCGGTTAAACCCGCCGTGGGCGGGTTTATTACTCCCCACCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgtg                              <  5:337877/133‑1 (MQ=255)
                                                                                            ttGGGGTGGGCCGG‑CACTTTTTAAC‑CGGCGGTGGGCTGTGTGATAACACCCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgtg                              <  4:65928/134‑1 (MQ=255)
                                                                                            tgggggTGGGACGG‑GA‑CGTGTAAACCCGCGGGGGGCGGGGTGATTACTACCCACCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGc                             <  2:263669/134‑1 (MQ=255)
                                                                                            tgggggGTGGCCGG‑GACCG‑GTAAACCCGCCGGGGGGTGGGTTATAACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGc                             <  2:295582/134‑1 (MQ=255)
                                                                                            aaggggGTGGTCGGCGACCTTTTAAA‑CCGCCGGGGGGTGGGTTATTACTACCCCGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAAtgtg                              <  1:94142/134‑1 (MQ=255)
                                                                                                 gtggcccGG‑CACCTTGAAAC‑CCGCGGTGGGCTGGGTTATTAAAACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGt                        <  1:252738/130‑1 (MQ=255)
                                                                                                    ggccGG‑CACCGTTTAAC‑CCGCGGGGGGGGGGGTTATAACACCCCCCCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCa                     <  7:335118/133‑1 (MQ=255)
                                                                                                                     acccgcggTGGGCGGGGTGATAACTACCCCCCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGc      <  6:287661/134‑1 (MQ=255)
                                                                                                                         gccgGGGGCTGGGTTATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGCGAAc  <  1:251078/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                         
AACTGGCCGCGCGCAGCATGACCGAGCAGCGGTTCGCCGAACTGGCAGGGATTTTTGGCCAGGCCTCCAGCCTGGAGGCGGATTTCTGGCAAATGGGGTTGGACGG‑GACCGTGTAAA‑CCGCCGTGGGCTGGGTGATTACTACCCAGCCCAGGACCGCTTTGGTGCATGGGCGCCTCAAAATAAGTCGTCCCGCACACTCCAGATCGTCAGATCGCAGCTTTAATGTGCACAGTGCATATGGATTTCTGGGCGAAC  >  ACB143/3608362‑3608616

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑