| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 3,422,037 | 0 | A | T | 23.8% | 13.6 / 10.7 | 21 | noncoding (85/1542 nt) | rrsD | 16S ribosomal RNA of rrnD operon |
| Reads supporting (aligned to +/- strand): ref base A (6/10); new base T (0/5); total (6/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.62e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAAT > CP009273/3421891‑3422160 | acaTCCGATGGCAAGAGTCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAg < 1:987713/149‑1 (MQ=21) ccaatcgCAAGAGGCCCGAGCGTCCCCCTCTTTGGTCTTGCCACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTACCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCaa < 2:916798/138‑1 (MQ=11) cGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATGAGCTACCGTTTCCAGTAGTTATCCCCCTCCAGCAGGCAGTTTCCCAGACCTTACTCACCCGTCCGCCACTCGTCCGCCAACCAGCAAg > 2:281817/1‑144 (MQ=11) aTGGCAAGAGGCCCGTGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTCTCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCt < 2:504266/144‑1 (MQ=17) tGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCTTGCTGTTTcg < 1:118771/149‑2 (MQ=12) ggCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCTTGCTGTTTcgc < 1:774691/149‑3 (MQ=12) cAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGttt > 2:177586/1‑144 (MQ=17) aaGAGGCCCGGGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTc < 2:1013772/144‑1 (MQ=17) aaGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGCTTc > 2:709709/1‑144 (MQ=17) cgagGCCCGCGGGACCCCCTCTTTGTTCCTTCGACGTTAGGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCTTGCTGTTTcg < 2:968368/143‑2 (MQ=11) agagGGCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATGAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCTTGCTGTTTcg < 2:116811/144‑2 (MQ=11) agagGCCCGGGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCTTGCTGTTTcg < 2:192026/144‑2 (MQ=12) agagGCCCGAGCGTCCCCCTCTTTGGTCTTGCTACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTcc < 2:336765/144‑1 (MQ=17) gCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCGGTAGTGATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACCCGTCCGCGAAACAGCAAGCTGTTTCCTGTTACCGt > 1:874203/1‑149 (MQ=255) cccGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCCCCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTa > 2:80155/1‑144 (MQ=17) ccGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTc > 1:188195/1‑149 (MQ=32) gAGGGTCCCCCTCTTTGGTGTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGCCATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGa < 1:553860/149‑1 (MQ=17) caGGGTCCCCCTCTTTGGTCTTGCGACGTTACACGGTATTACCTACCGTTTCCAGTAGTTATCCCCCGCCATCAGGCAGTTTCCCAGTCATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGa < 1:79562/148‑1 (MQ=11) aGGGTCCCCCTCTTGGGTCTTGCGACGTTATGCGGTATTAGCCACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGt < 2:6096/144‑1 (MQ=14) caccAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACCCGTCAGCGAAACAGCAAGCTGTTTCCGGTTACCGTTCGACCTGCATGTGTTAGGCCTGCCGCCAGCGTTCATTCTGAGCCATGATCAAACTCGTCAATTTAAAAGttt < 1:72962/146‑1 (MQ=1) catgCAGTTTCCCAGACACTCCTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTTTTCCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGAt < 1:373695/146‑1 (MQ=1) ccccTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAAt < 2:78125/141‑1 (MQ=1) | ACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAAT > CP009273/3421891‑3422160 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |