Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,556,207 | 0 | T | C | 48.6% | 30.4 / 18.0 | 37 | S215G (AGT→GGT) | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base T (18/1); new base C (0/18); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.13e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA > NC_000913/3556058‑3556355 | tgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGGTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTc < 1:197252/151‑1 (MQ=255) gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCg > 1:276360/1‑151 (MQ=255) cctttttcccccGACAAAAAAAATGCTTTGGGGGATATTTCCTTTTCCGATTAAAAGCGAGCGGGTATTCCCCGGCCCC‑TGGCCGCGCGGCTGTTTATAAATTTTTT‑GTTAAGGCGGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGc < 2:362857/141‑1 (MQ=255) cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCa > 2:692181/1‑151 (MQ=255) aaagGCGTTGGGTGTTTTTTTCTTTTCCGCCTTCAAGCGAGGGGTTATTCCCCGCCCCC‑GGGTCGCGCGGCAGTTTATGAATTTTTT‑TTTCATGCGGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac < 2:276463/151‑1 (MQ=255) cGTCCGGTGATGTTTTCACTTCCGCCTACAGGCGAGAGGATATT‑CCCGGCCCCTTGGAGCCGCGGCAGGTTATAAATATTCT‑GTCCAGGGAGGAAAAAT‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg < 1:680498/151‑1 (MQ=255) ggTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCCGCTAAt > 1:288669/1‑151 (MQ=255) taatatTTTCACTTCCGCTTAAAAGGGGGCCGGT‑TTACCCCGCCCCTTGGCCGCGCGGCAGGTTATGAAT‑TTCTGGTCCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg < 1:410476/149‑1 (MQ=255) gATATTGTCCCTTTCAACTACAAGCGAGACGGTATAACCCGGCCCC‑GGGTCGCGGGGCAGTTTAAGAATATTCT‑GTACAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg < 1:457910/150‑1 (MQ=255) aattttttCACTTTCGACTTCAAGCGGGACGTTATTCCCCGCCCCT‑TGGTCGCGCGGCAGTTTAGAATTTTTCT‑TTTCA‑GGAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg < 1:657223/146‑1 (MQ=255) attttttCCTTTTCGACTTCAAGCGGGGGGGTATTCCCCGGCCCC‑GGGCCGCGCGGCGGGTTATAAATTTTTT‑GTTCATGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa < 2:346717/148‑1 (MQ=255) tatTTTCACTTTCGACTTCAAGCGAGACGGTATTAACCGGCCCC‑TGGTCGCGCGGAAGGTTAAGAATAATCT‑GTTAATGCAGGGACAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa > 1:512394/1‑150 (MQ=255) tatTTCCACTTCCACCTTAAAGAGAGCCGTTATTACCCGCCCCC‑GGGTCCCGCGACAGTTTATAATTATTCT‑GTTAAAGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAGGGCGCGGCACACCAGCTAATAGAAc < 1:688293/151‑1 (MQ=255) aCTTCAAGCGAGCCGGTATTCCCCGCCCCC‑TGGTCGCGCGGCAGTTTATAAATTTTCT‑GTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATc < 2:314951/150‑1 (MQ=255) aaaGCGAGACGGTATTCCCCGCCCCG‑TGGTCGCGGGGAGGTTTATGAATTTTCT‑GTTCATGAAGGGAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac < 1:356594/150‑1 (MQ=255) ggTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg > 2:688208/1‑151 (MQ=255) gTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg > 1:659140/1‑150 (MQ=255) gTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg > 2:277500/1‑150 (MQ=255) gTATTACCCCGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGCTTTCCCGCTCACCAACCTGGCGCGGCACACCAGCTAATAGCACTTCTCCACGCATATGCACAATGTTCCcg > 1:35722/1‑150 (MQ=255) gCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATCTTCCCGCGCTCGCCAAg > 1:696678/1‑151 (MQ=255) ccgcgcgGCAGGTTATAAATTTTCT‑TTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAgg < 2:503936/150‑1 (MQ=255) ggcgggGCAGTTTATAAATTTTGT‑TTATAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAgg < 1:33992/149‑1 (MQ=255) ctgtttATGAATATTCT‑ATTCATACAGAGAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc < 1:659409/147‑1 (MQ=255) aGGTTATAATTTTTTT‑GTACATGGGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc < 2:151853/150‑1 (MQ=255) ggTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa > 1:100984/1‑151 (MQ=255) tttataaattctgt‑GTTCAGCCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc < 1:397718/135‑1 (MQ=255) tttagaaaTATTCT‑GTTCAGGGAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc < 2:282857/142‑1 (MQ=255) gTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 2:747165/1‑151 (MQ=255) aTGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc > 2:342446/1‑150 (MQ=255) tGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa > 2:481152/1‑151 (MQ=255) tGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGCGTTAAACGATGCCa > 2:153167/1‑151 (MQ=255) ttCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACt < 2:561453/151‑1 (MQ=255) tccacaggggAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACt < 1:743027/144‑1 (MQ=255) gggAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTg > 1:792473/1‑150 (MQ=255) gAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACCACGCTCCCGCCCTCGCCAAGTTGCAAGGCGTTAAACGATGCCACGGTCGAGACTTCCGTTGcaa > 1:517747/1‑149 (MQ=255) gAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGcca > 2:159625/1‑151 (MQ=255) cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACcaca > 1:671720/1‑150 (MQ=255) | TGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA > NC_000913/3556058‑3556355 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |