| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 505,440 | 0 | C | A | 50.0% | 11.6 / 14.2 | 20 | A176E (GCA→GAA) | ushA | bifunctional UDP‑sugar hydrolase/5'‑nucleotidase |
| Reads supporting (aligned to +/- strand): ref base C (9/1); new base A (0/10); total (9/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.19e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACC‑GTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGC > NC_000913/505291‑505589 | aCGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACC‑GTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGc < 2:192138/150‑1 (MQ=255) gggCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACC‑GTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGa < 1:157816/151‑1 (MQ=255) ccAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACC‑GTGGGCGCTGTTTAAGCGTCAGGATCTGACAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAAt > 1:87799/1‑150 (MQ=255) ttCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACC‑GTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTc > 2:72644/1‑150 (MQ=255) cGAGCGCCTGTTTAAACC‑GTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTg > 2:185833/1‑151 (MQ=255) ttAAACCGGGGGGCTGTGTTTAAGCGGCGGGATTAAAAAATTG‑CCGTTATTGGCCTGACAACCGTTAACCACGAAAAAATGGGAACCCCGGAAGACTTCACTT‑AT‑T‑TGAATTTTGTTAGGCCCGCCGATGAAGCGAAGCTGGt < 1:248339/143‑1 (MQ=255) gTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGGGATTCAGGAGCTGc > 1:42809/1‑149 (MQ=255) gTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAAcagacaga > 2:1225/1‑150 (MQ=255) tttACGCTTCAGGATATAAAAATTG‑CCTTTATTGGCTTAAAACCCGAGGAAACAGAAAAAATTGGTAGCCCGGAAAACTTCACTG‑ATAA‑CGATTTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAacagacag < 1:50444/148‑1 (MQ=255) aaGCGTCAGGATTTAAAAATTG‑CCTTTATTGGGCTAACACCCGATGACCCAGAAAAAATTGGTACCCCGAAATCCTTCACTG‑ATAT‑CGATTTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGaaaa < 2:290445/149‑1 (MQ=255) aaGCGTCAGAATCTAAAAATTG‑CCGTTATGGGGCTGACACCCAATGACACAGAAAAATTTGGTACCCCGAAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGc < 2:368862/151‑1 (MQ=255) aagtgtcaggtgataaaaTTGCCCTTTTTTGGGCTGACACCGAAACACAAAAAAAAAATGGGTAACCAGAAATATTTC‑CTGAATAT‑GGATTTTCG‑AAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGc < 1:216343/137‑1 (MQ=255) gtaaGGAACTAAAATTTG‑CCTTAATTGGGCAGACAAACGATGAAGCAAAAAAACTAGGTACCGGAGAATTGTTCGCGA‑ATAT‑GGAATTTAG‑AAAGCCCGCCGATGAAGCGAAGCTGATGATTCAGGAGCTGCAACAGAGAGAAAATCCAGa < 1:72625/148‑1 (MQ=255) aGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACat > 2:203669/1‑151 (MQ=255) ctaaaaaTTG‑CTTTTATGGGGCTGACACCCAATGACACAGAAAAAATTGGTACCCCGAAATACTTCAC‑G‑ATTTGCGAATTTTG‑AAGGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACattat < 2:184395/146‑1 (MQ=255) tGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATtattat > 2:328735/1‑151 (MQ=255) tGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATtattat > 1:24581/1‑151 (MQ=255) aaaaTTG‑CCGTTATTGGGCTGACAACCGATGACACAGAAAAAAGTGGTAACCAGGAATACTTCACTA‑ATAT‑CGAATTTCA‑AAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATcg < 1:67362/151‑1 (MQ=255) gggCTAACACCCGATGACACACAAAAAATTGATCCCCCGTAATCTTCCACTA‑ATAT‑AATTTTTCG‑TATGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATg < 2:201170/150‑1 (MQ=255) acccGAAGACACAGAAAAAATTGGAACCCCGGAATACTTCACTG‑ATAT‑CAATTTTCG‑AAGGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTa < 2:208183/148‑1 (MQ=255) gACACAGCAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATg > 2:351477/1‑151 (MQ=255) cAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGc > 1:357766/1‑150 (MQ=255) | ACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACC‑GTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTG‑ATAT‑CGAATTTCG‑TAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGC > NC_000913/505291‑505589 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |