| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,556,207 | 0 | T | C | 51.8% | 16.3 / 36.0 | 27 | S215G (AGT→GGT) | rtcA | RNA 3'‑terminal phosphate cyclase |
| Reads supporting (aligned to +/- strand): ref base T (11/2); new base C (0/14); total (11/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.98e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-04 | |||||||||||
GCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT > NC_000913/3556074‑3556345 | gCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGa > 2:485717/1‑151 (MQ=255) acgctttttcCACCAACGACAAAAAAGCGTGCGGAGATATTTCCCATGCCGACTAAAAGCGAGCGGGTATCCCCCGCCCCGTGGTCGCGCGGCAGTTTATAAATTTTCTTTTC‑TGGGGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCa < 1:19229/144‑1 (MQ=255) cctcTTTTCACCGACGACAAAAAAGCGTGCGGTGTTATTTTCTCTTTCGACTTAAAGCGAGCCGTTATTCCCGGCCCCCGGGTCGCGCGGCAGTTTATAAATATTCTTTTAAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGc < 1:499903/148‑1 (MQ=255) cTTTTCACCGACGACAAAAAAGCGTTAGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAGTATTCCGTTCATGCAGGGAAAAA‑CTCCCCGCCAGGGTAGCGATTTCACGCTCa > 2:455172/1‑151 (MQ=255) ccccGACAAAAAAAAAGCTTTGGGGGATTTTTCCGCTTTCGACTTCAAGGGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTAAGAATTTTCTGTTCATGAAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc < 2:68770/149‑1 (MQ=255) accacAAAAAGGCGTTCGGTGATTTTTTCCCTTCCGACTTAAAGCGAGACTGTATTCCCCGGCCCCTGGCCGCGCGGCGGTTTATGAATTTTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg < 1:233178/147‑1 (MQ=255) aaaaaaaaaGGTTTGGGTGATATTTTCACTTCCACCTACAGGCGAGACGTTATCCCCCGCCCCCTGGCCGCGCGGCAGGTTATGAATTTTCTTTTCATGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc < 2:552742/148‑1 (MQ=255) gcgtggggtttttttttCCCTTTCGACGTAAAGCGAGACGTTATTCCCCGGCCCCTGGCCGCGCGGCGGGTTATGAATTTTCTTTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc < 1:607934/137‑1 (MQ=255) ttgggTGATTTTTTCCCTTCCGATTTAAAGCGAGACGTTTTTCCCCGCCCCTGGGTCGCGCGGAGGTTTATAAATTTTCAGTTCATGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt < 2:450924/148‑1 (MQ=255) ttCGGGGATTTTTTCCTTTCCCACTTAAAGCGAGCCGGTATTCCCCGGCCCCGGGGCGCGCGGGGGGTTAAAAATTTTCTGTTCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt < 2:213120/151‑1 (MQ=255) cGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa > 2:486626/1‑151 (MQ=255) aattcaaaGGGAAACGGTATTCCCCCCCCCCGGTCCCCGGGGCAGGTTATAAATTTTATTTCCATGGGGAGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATc < 2:594440/145‑1 (MQ=255) aaaacgggtgggTTTTCCCGGGCCCCTTGCGGCCGGGCGGTTTAAGAATGTTCTGTTAAGGCAGGAAAACC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac < 2:600101/141‑1 (MQ=255) gTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg > 1:498682/1‑150 (MQ=255) gTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg > 2:237176/1‑150 (MQ=255) atttcccccGCCCCCTGGTCGCGCGCCAGTTTATAAATTTTCTGTTACTGCGGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg < 2:511887/144‑1 (MQ=255) ccgcgcgGCGGTTTATAAATTTTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAgg < 2:529549/150‑1 (MQ=255) gtttATAATTATTTTGTTCA‑GGAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg < 1:587395/148‑1 (MQ=255) ggTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa > 1:123670/1‑151 (MQ=255) gTTATGAATCTTCTTTTCATGGAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc < 2:147166/148‑1 (MQ=255) gTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGTTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTCCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 2:158219/1‑151 (MQ=255) gTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 2:421299/1‑151 (MQ=255) gTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 2:115117/1‑151 (MQ=255) tataAATTTTCTGTCCAAGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc < 2:353128/151‑1 (MQ=255) tGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGTCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCACGCCAAGTTGCAAGGTGTTACACGATGCCa > 1:88540/1‑151 (MQ=255) gTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgaga > 1:501512/1‑150 (MQ=255) cgggAAAAACCCCCCCGCCAGGGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt < 2:156539/149‑1 (MQ=255) | GCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT > NC_000913/3556074‑3556345 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |