Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,556,207 | 0 | T | C | 34.6% | 31.8 / 10.0 | 26 | S215G (AGT→GGT) ‡ | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base T (15/2); new base C (0/9); total (15/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.76e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC > NC_000913/3556074‑3556343 | gCACTGCCGCGCTTTTCCCCGACGACAAAAAGGCGTTCGGTGATATTTTCCCTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGCCGCGCGGCAGGTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGc < 2:23603/140‑1 (MQ=255) gcgtTTTTCCCCGACGAAAAAAAAGCGTTGGGTGATATTTTCACTTTCGATTTCAAGCGAGACGTTATTCCCCGGCCCCGGGCCGCGGGGCAGTTTATGAATTTTCTGTTCATGCGGGAAAAACCCCCCCGCCAGTGTAg < 2:400906/140‑1 (MQ=255) cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGa > 1:98354/1‑140 (MQ=255) ccGACAAAAAAAAAGCGTTGGGTGATTTTTTCTTTTCCGACTTAAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACg < 2:836750/140‑1 (MQ=255) aacaaaaaaaagGTGTTCGGTAATTTTTTCATTTCCGACTTCAAGCGAGCCGTTTTTCCCCGGCCCTGGGCCGCGGGGGAGTTTATAAATTTTTTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCt < 2:196624/134‑1 (MQ=255) acAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc > 2:56471/1‑140 (MQ=255) aggCGTTCGGTATTTTTTTCCCTTTCGCTTAAAAGGGGGACGTTATTCCCCGCCCCCTGGCCGCGCGGCAGTTTATAAATATTTTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg < 2:803073/138‑1 (MQ=255) aGCGTTCGGTGATTTTTTCACTTTCGACTAAAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg < 2:497776/140‑1 (MQ=255) ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg > 1:740979/1‑140 (MQ=255) tCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg > 2:628929/1‑140 (MQ=255) ttttttCTCTTTCCAGTCAAGGCGAGACGTTATCCCCGGGCCCTTGGTCGCGCGGCTGTTTATGAATATTTTGTGCAGGGGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa < 2:98354/138‑1 (MQ=255) tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa > 1:242868/1‑140 (MQ=255) tCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa > 1:1848/1‑140 (MQ=255) aGCGAGCCGTTTTTCCCCGGCCCCTGGCCGCGGGGCAGTTTATGAATTTTCTTTTCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg < 2:347122/140‑1 (MQ=255) gagcCGGTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATAAATTTTTTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt < 2:1848/140‑1 (MQ=255) gaCGGTATCACCCGCCCCCTGGTCGCGCGGCAGGTTATGAATTTTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt < 2:59216/140‑1 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 1:13830/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 1:574374/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 1:455210/1‑140 (MQ=255) ggccgcgggGCAGTTTATAAATATTCTTTTCAGGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg < 1:628929/137‑1 (MQ=255) ggTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg > 1:230835/1‑140 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 2:803006/1‑140 (MQ=255) cAGTTTATGAATATTCTGTCCAGGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 1:411319/140‑1 (MQ=255) cAGGTTATAAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 2:13830/140‑1 (MQ=255) ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa > 1:790090/1‑140 (MQ=255) ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa > 2:581628/1‑140 (MQ=255) tCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 1:542234/1‑140 (MQ=255) gTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc > 2:795297/1‑140 (MQ=255) ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa < 1:795297/140‑1 (MQ=255) cggggaGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg < 2:230835/135‑1 (MQ=255) aaCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc < 2:542234/140‑1 (MQ=255) | GCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC > NC_000913/3556074‑3556343 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |