| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,173,364 | 0 | C | . | 64.7% | 15.5 / 19.6 | 17 | intergenic (‑491/+918) | gatD/gatB | galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component |
| Reads supporting (aligned to +/- strand): ref base C (2/4); new base . (7/4); total (9/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.35e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACG‑ATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATC > NC_000913/2173141‑2173580 | gtgtGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACG‑ATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCTGATc > 1:224841/1‑251 (MQ=255) cccccacccgaccctctcccgttcggATCCACAGGGGGATGCTCAAAAGTAGGACCCCCGAGATTAAGGGGCGGAACAGAGTTCCACG‑AGGCCCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGATTGGGATAAAAAGCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc < 2:16189/225‑1 (MQ=255) ccAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACG‑ATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt > 1:81579/1‑251 (MQ=255) aaTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACG‑ATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc > 1:262312/1‑250 (MQ=255) gATCCCCAGGGTGATGCTCATAATGATGACACTCGAGATTAAGGTGCGGAACAGATTTCCACG‑ATACACGGCGACCGCCATAGCCACGAAGAAGCCGATGGTGGAAAGATCGCCCACCGGCAGCACCTGATTACCCGGCACACAAACCGCAATTAAAATGATGAGTGGGATAAAAATCAGACTTGCCGATAC‑‑CCACGCCGTATTTCCCAGCAGCACTGCAGGATCAAGGCCAATAAGGAACTCCTGAccgc > 2:213609/1‑251 (MQ=255) gATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACG‑ATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc > 1:16189/1‑225 (MQ=255) acaGGGTGAGGCTCATAATGATGCCACCCGAGATTAAGGTGCGGAACAGATTTCCACG‑ATGCACGGCGACCGCCATCGCAACGACGAAGCCGTTGGTGGCAAGCTCGCCATACAGCAGCACCTGATTACCCGGCACACAAAAAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGACGGATCAAGGCCAATCAGGAAATCATGACTGCCGaa > 2:227375/1‑250 (MQ=255) aaattagtcccaccaggaTTAAGGTCCGGAACAGTTTTCCACA‑ATCCCCGGCCCCCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGTTGAGTGGGATAAAAACCAGCCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa < 2:81579/235‑1 (MQ=255) ataTTAAGGTGCGGAAAAGATTTCC‑CGTATGCACGCCGACCTCCATCGCCACGAAGAATCCGATGGTGGAAAGATCGCCAAACGACAGCTCCTGATTACCCGGCACACAAACACCAATTAAAATGGTAAGTGGTATAAAAATCAGGCTGTCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGTCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGAATACGAGCTTGCTTa < 1:183167/249‑1 (MQ=255) gggAACAGATTCCCGCG‑ATGATCGGCACCCGCCTTCGCCCCGAAGAAGCCGATGTGGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTAAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGt < 1:131414/250‑1 (MQ=255) cACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt > 2:156196/1‑251 (MQ=255) cccccATCCCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCGAATTCCCCGGCACACAAACAGCAAATAAAATGGTGAGTGGAATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc < 1:53817/248‑1 (MQ=255) aaaaaagCGGATGGGGGCAGGACCGCCAAAGGGCAGCACCTGATTACCCGGAACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCGGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagca < 2:264444/247‑1 (MQ=255) gaaaGATCGCAAAACGGCAGCCCCTGATTCCCCGGCACACAAACAGCATTTAAAAGGGTGAGGGGAATAAAAATCGGTGTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGt < 1:310717/249‑1 (MQ=255) ccGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGAAAAAAACCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAc < 1:41406/250‑1 (MQ=255) aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc > 1:273109/1‑251 (MQ=255) aGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc > 1:2104/1‑249 (MQ=255) | GTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACG‑ATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATC > NC_000913/2173141‑2173580 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |