Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,556,207 | 0 | T | C | 43.8% | 62.5 / 13.9 | 48 | S215G (AGT→GGT) | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base T (22/5); new base C (0/21); total (22/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.22e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA > NC_000913/3556059‑3556355 | gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTcc < 1:725141/151‑1 (MQ=255) cGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTcccc < 2:21358/151‑1 (MQ=255) caccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgt > 2:909987/1‑151 (MQ=255) cccgcGGCACTGACGCGTTTTCCACCGACGAAAAAAAAGCGTTG‑GGTGA‑TATTTTCTTTTTCGACTTAAAGGGAGACGTTATTCCCCGCCCCCGGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAgtg < 2:291117/147‑1 (MQ=255) cgcgTTTTCCCCCGACGAAAAAAAGGCGTTC‑GGTGA‑TTTTTCCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATTTTTTGTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTc < 1:541084/148‑1 (MQ=255) ccctcTTTCCCCCGACAAAAAAAAGGCGTCC‑GGTAA‑TTTTTTCCCTTCCGACTCCAAGCGAGAGGTTATCCCCCGCCCCCTGGCCCCGCGGCGGTTTATAAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTc < 2:931736/144‑1 (MQ=255) gctcTTTCCACGGACGAAAAAAAAGCG‑TC‑GGTTATTTTTTTCCCTTTCGACTAAAAGCGAGCCGTTATTCCCCGCCCCCTGGCCGCGGGGAAGTTTAGGAATATTCTTTTCAAGCAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTc < 2:552723/144‑1 (MQ=255) gcgcTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGc > 1:877924/1‑151 (MQ=255) gcTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGAt > 1:874032/1‑142 (MQ=255) tCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCa > 1:700867/1‑150 (MQ=255) cgacgaCAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg > 2:173524/1‑151 (MQ=255) cgacgaCAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg > 1:570266/1‑151 (MQ=255) cgacAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg > 1:367739/1‑151 (MQ=255) acAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg > 1:446185/1‑151 (MQ=255) aaaaaaaaGCGTTC‑GGTGA‑TTTTTCCCCTTTCGACTTAAAGCGAGACGTTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc < 1:105634/148‑1 (MQ=255) aaaaaaGGGGTCGGGTGA‑TTTTTTCCTTTTCGACTAAAAGCGGGACTTTTTCCCCCGCCCCTTGGCCGCGCGGCGGTTTATGAATTTTCTGTTCAGGAAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt < 2:553473/144‑1 (MQ=255) aaaaaGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc > 1:916123/1‑150 (MQ=255) tC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt > 2:593195/1‑150 (MQ=255) ggtattttttccccTTACGCCTTCAAGCGAGACGTTTTTCCCCGCCCCCGGGTCGCGGGGGAGTTTATGAATTTTTTTTTCATGGAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt < 2:726285/135‑1 (MQ=255) tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc > 1:798944/1‑151 (MQ=255) ttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct > 2:1039024/1‑151 (MQ=255) ttccGACTTCAGGGGAGACGGTTTTCCCCGGCCCCAGGCCGCGGGGCAGTTAAGGAATTTTCTGTTCAGGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc < 2:567272/145‑1 (MQ=255) tcggcctaaaggggagaGGGTATCCCCCGGCCCCGGGGCCGGCGGAAGTTTATGATTTTTTTGTTCAGGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg < 1:155110/137‑1 (MQ=255) ttaaaGCGAGAGGTTATACCCCGCCCCCGGGCCGCGCGGAGGGTTATAAATTTTCTGTTGATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg < 1:910028/147‑1 (MQ=255) aaGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac > 2:932443/1‑150 (MQ=255) gagaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt > 2:316510/1‑151 (MQ=255) gaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCGCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTc < 1:158158/151‑1 (MQ=255) acgttatccccgggCCCCGGGGCGCGCGCCTGTTTATAAATTTTTTTTTCTTGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt < 1:587039/136‑1 (MQ=255) gttATTCCCGGCCCCCGGGTCGCGGGGCAGTTTAGGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt < 2:1013311/145‑1 (MQ=255) tcccccGCCCCCGGGTCGCGGGGCAGGTTATAAATATTCTGTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgcg < 1:603725/146‑1 (MQ=255) tACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc > 2:878023/1‑151 (MQ=255) cccGGCCCCTGGTCGCGCGGAAGTTTATGAATTTTCTGTCCAGGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc < 1:515249/151‑1 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgt > 1:894650/1‑151 (MQ=255) gcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt > 2:134258/1‑151 (MQ=255) cgGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt < 1:718035/150‑1 (MQ=255) cgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt > 1:519077/1‑150 (MQ=255) ggaagtttATAAATTTTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa < 1:804101/145‑1 (MQ=255) ggCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa < 2:52851/151‑1 (MQ=255) ggCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa > 1:122788/1‑150 (MQ=255) gcgaGGTAATAAATTTTCTGTTCGGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa < 1:1039073/140‑1 (MQ=255) aaTATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg > 1:437677/1‑151 (MQ=255) tatTCTGTCCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg < 1:721262/150‑1 (MQ=255) ttCTGTTCATGGAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg < 1:617539/149‑1 (MQ=255) tCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt > 2:920011/1‑151 (MQ=255) gaaGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc < 2:298936/146‑1 (MQ=255) ggaaaaac‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg < 2:831869/144‑1 (MQ=255) aaaaa‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccac > 1:733956/1‑151 (MQ=255) cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACcaca < 2:367718/150‑1 (MQ=255) | GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA > NC_000913/3556059‑3556355 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |