| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,517,435 | 0 | C | G | 34.5% | 47.7 / 11.4 | 29 | A16G (GCA→GGA) | ortT | orphan toxin OrtT |
| Reads supporting (aligned to +/- strand): ref base C (3/16); new base G (10/0); total (13/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.43e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTT > NC_000913/1517370‑1517501 | tCTGGCTGGACGGCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGaa > 1:2119812/1‑68 (MQ=255) tggctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATc < 1:2874558/68‑1 (MQ=255) ggACGCCCGTACATGTCTCTCTCTCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCGATCGcgct > 1:2260955/1‑65 (MQ=255) ggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCGGCAATCGCAtt < 1:1290503/68‑1 (MQ=255) aCGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTc < 1:461299/68‑1 (MQ=255) cGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAAATGCAATTCt > 1:525382/1‑68 (MQ=255) cGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCt < 1:2878715/68‑1 (MQ=255) cGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCt < 1:2597644/68‑1 (MQ=255) cGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCt < 1:1950692/68‑1 (MQ=255) agCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCt < 1:508143/67‑1 (MQ=255) ccGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTAt < 1:177803/68‑1 (MQ=255) gTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCa < 1:3004837/68‑1 (MQ=255) aCATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCAcc < 1:70453/68‑1 (MQ=255) aCATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCAcc < 1:2776754/68‑1 (MQ=255) aTGTCTCTCTATCAACACATGCTTGTTTTTTATGCTGTTATGGCAGCAAACGCATTTCTTTTTAccct > 1:2649437/1‑66 (MQ=255) aTGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAATCGCATTTTTTTTCACCTg > 1:3373528/1‑68 (MQ=255) aTGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAAAGGCGGTTTTTCTCACCta > 1:2556823/1‑67 (MQ=255) aTGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAAACGCCTTTTTTTTTAcccg > 1:2786790/1‑66 (MQ=255) aTGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAAACGCCTTTTTTTTTACCTg > 1:2090224/1‑68 (MQ=255) aTGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAAACGCATTTTTTTTTACCTg > 1:2795020/1‑68 (MQ=255) aTGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTTTTATCAcccg > 1:670695/1‑66 (MQ=255) aTGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAAACGCCTTTTTTTTAACCTg > 1:525793/1‑68 (MQ=255) aTGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAAACGCAATTCTTTTTAcccg > 1:2068781/1‑66 (MQ=255) aTGTCTCTCTATCAACACATGCATGTTTTTTATGCGGTTATGGCAGCAATCGCATTTTTTATTACCTg > 1:2319632/1‑68 (MQ=39) tctctcTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAAAAGCATTTTTTTTTACCTgggt > 1:1821978/1‑66 (MQ=255) tctctcTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTTTTACCTGGtt > 1:959926/1‑68 (MQ=255) tctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAATAGCATTTTTTTTTACCTGGcttt > 1:277717/1‑67 (MQ=255) ctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAAACTCATTTCTTTTCACCTGGtttct > 1:338000/1‑68 (MQ=255) ctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGGAAACGCATTTTTTTTCACCTGGtttct > 1:1425276/1‑68 (MQ=255) aTCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAAACGCCTTTCTTATCACCTCGctttttttt > 1:1433086/1‑66 (MQ=255) cacaTGCTTGTTTTTTATGCGGTTATGGCAGCAAAAGGATTTCTTATTACCTGGGTTTTTTTTTACGa > 1:3199275/1‑68 (MQ=255) caTGCTTGTTTTTTATGCGGTTATGGCAGCGATCGGATTTTTTTTTACCCGGGTTTTTTTTTACGCTa > 1:2193434/1‑68 (MQ=255) caTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCCTTTTTTATCACCCGGTTTCTTTTTCCCGATa > 1:2809435/1‑68 (MQ=37) tGCTTGTTTTTTATGCGGTTATGGCAGGAATCGCATTTTTTATTACCCGGTTTCTTTCTCACGATAAg > 1:292624/1‑68 (MQ=255) cTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGaa < 1:320831/68‑1 (MQ=255) tttttATGCGGTTATGGCAGGAATCGCATTCCTTTTCACCTGGTTTCTTTCTCACGATAAAAAACGCa > 1:122450/1‑68 (MQ=255) tttttATGCGGTTATGGCAGGAAACGCATTTTTTATCACCTTGTTTTTTTCTCACGATAAGAAAAGCa > 1:1797070/1‑68 (MQ=255) tttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCa < 1:1733543/68‑1 (MQ=255) tATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCg < 1:1495044/68‑1 (MQ=255) tGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGtt < 1:1337843/68‑1 (MQ=255) tGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgc < 1:112311/68‑1 (MQ=255) cagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGcc < 1:2379874/68‑1 (MQ=255) cagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGcc < 1:406766/68‑1 (MQ=255) cagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGcc < 1:2212121/68‑1 (MQ=255) cagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGcc < 1:2128217/68‑1 (MQ=255) cagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGcc < 1:1418075/68‑1 (MQ=255) cagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGcc < 1:1394759/68‑1 (MQ=255) gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCtt < 1:3320439/68‑1 (MQ=255) gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCtt < 1:201994/68‑1 (MQ=255) gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCtt < 1:1281948/68‑1 (MQ=255) gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCtt < 1:1055376/68‑1 (MQ=255) | TCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTT > NC_000913/1517370‑1517501 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |