| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,171,557 | 0 | C | A | 57.5% | 0.0 / 58.9 | 40 | intergenic (+13/‑80) | rrlB/rrfB | 23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon |
| Reads supporting (aligned to +/- strand): ref base C (9/8); new base A (14/9); total (23/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.49e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAA > NC_000913/4171510‑4171689 | cTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa > 1:879812/1‑139 (MQ=11) aCCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa > 2:1594289/1‑139 (MQ=33) aCCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa > 1:130/1‑139 (MQ=33) aCCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa > 1:66334/1‑139 (MQ=17) ccGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAg < 2:1001917/139‑1 (MQ=17) cGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgc > 1:1446078/1‑139 (MQ=32) tACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg > 2:591440/1‑139 (MQ=33) aCTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCgg < 1:1594289/139‑1 (MQ=33) aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtggt < 2:1446078/139‑1 (MQ=33) aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtggt < 2:1288875/139‑1 (MQ=32) gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTcc < 2:257510/139‑1 (MQ=33) gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGAATAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTGGCCTCGCGGCAGTAGCGCGGTGGTcc < 1:465026/139‑1 (MQ=12) ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAc > 1:1414773/1‑139 (MQ=33) ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa < 2:147459/124‑1 (MQ=18) ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa > 1:147459/1‑124 (MQ=18) ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg < 2:486745/129‑1 (MQ=17) ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg > 1:486745/1‑129 (MQ=14) ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAc > 2:849898/1‑139 (MQ=18) ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAc > 1:197446/1‑139 (MQ=18) gTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTcc < 1:1042090/134‑1 (MQ=18) gTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGATGGTcc > 2:1042090/1‑134 (MQ=17) tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGTGAAGATTTTCTGCCTGATACAGATTAAATCAGAACGCAGAAGCGCTCCGATAAAACAGAATTTGATTGGCGGCAGCACCGAGGTGGTCCTACCTg > 1:1592358/1‑139 (MQ=11) ggCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAcc < 1:1203666/139‑1 (MQ=17) ccTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGcc < 2:392999/139‑1 (MQ=32) ccTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGcc < 1:1190227/139‑1 (MQ=18) aCAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAAc > 2:271492/1‑139 (MQ=32) aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg < 1:1378251/112‑1 (MQ=17) aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg > 2:1378251/1‑112 (MQ=18) aCAACGCCGAAGAAGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAAc > 2:598145/1‑139 (MQ=18) cAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCg > 2:1131277/1‑135 (MQ=32) cAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCg < 1:1131277/135‑1 (MQ=32) cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCTGCCTGATACAAATTAAATCAGAACGCGGAACCGGTCTGATAAAACAGAATGCGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGACGAACTCAg > 2:1022929/1‑139 (MQ=11) ccGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggt > 2:860137/1‑108 (MQ=32) ccGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggt < 1:860137/108‑1 (MQ=33) ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc > 2:1173914/1‑112 (MQ=14) ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc < 1:1173914/112‑1 (MQ=18) ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACt > 2:575472/1‑134 (MQ=17) ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACt < 1:575472/134‑1 (MQ=17) ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGaa > 1:986137/1‑139 (MQ=17) ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGaa > 1:412849/1‑139 (MQ=17) | CTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAA > NC_000913/4171510‑4171689 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |