| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,763,636 | 0 | A | G | 35.1% | 66.7 / 77.2 | 74 | noncoding (208/1542 nt) | rrsG | 16S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base A (28/20); new base G (12/14); total (40/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.39e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATC > minE/1763567‑1763706 | gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa < 1:1001247/70‑1 (MQ=255)gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa < 1:1878346/70‑1 (MQ=255)gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa < 1:1677902/70‑1 (MQ=255)gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa < 1:557538/70‑1 (MQ=255)gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg > 1:3856221/1‑71 (MQ=255) cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGt < 1:313376/71‑1 (MQ=255) ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc > 1:2616534/1‑71 (MQ=255) tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc > 1:581392/1‑71 (MQ=255) tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc > 1:1736038/1‑71 (MQ=255) tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc > 1:1052480/1‑71 (MQ=40) ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccc > 1:3246388/1‑71 (MQ=255) ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccc > 1:184395/1‑71 (MQ=255) ggTGAGCCGTTACCCCACCAAAAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg > 1:3235574/1‑66 (MQ=255) tGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccccc < 1:1012355/71‑1 (MQ=255) aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc < 1:1033662/71‑1 (MQ=255) aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc < 1:3350105/71‑1 (MQ=255) aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc < 1:4096439/71‑1 (MQ=39) aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc < 1:4091749/71‑1 (MQ=39) aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct < 1:106070/70‑1 (MQ=40) gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc < 1:2517417/70‑1 (MQ=255) gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctct > 1:1359693/1‑71 (MQ=39) gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctct > 1:732055/1‑71 (MQ=39) ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCtt < 1:2957496/71‑1 (MQ=255) ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCtt < 1:2744803/71‑1 (MQ=255) ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCtt < 1:3524526/71‑1 (MQ=39) cGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt < 1:715685/71‑1 (MQ=39) gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg > 1:1277200/1‑71 (MQ=255) ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCATTgg > 1:132760/1‑71 (MQ=255) aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa > 1:3105337/1‑54 (MQ=255) aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa > 1:2367444/1‑54 (MQ=255) aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc > 1:3784956/1‑71 (MQ=255) aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc > 1:1359163/1‑59 (MQ=255) aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc > 1:1077384/1‑59 (MQ=255) aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt < 1:2725771/67‑1 (MQ=40) ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc < 1:3530224/70‑1 (MQ=255) ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt < 1:432687/71‑1 (MQ=255) ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt < 1:712105/71‑1 (MQ=255) ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc < 1:3249292/57‑1 (MQ=255) caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc < 1:3235062/71‑1 (MQ=255) accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc > 1:2869419/1‑70 (MQ=255) accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCg < 1:2869063/71‑1 (MQ=255) accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCg < 1:332584/71‑1 (MQ=255) accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCg < 1:2469072/71‑1 (MQ=255) accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt > 1:2745423/1‑52 (MQ=255) ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc > 1:2713356/1‑69 (MQ=255) ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGa > 1:915783/1‑71 (MQ=255) caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAc > 1:2961804/1‑71 (MQ=255) aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccccc > 1:3969825/1‑54 (MQ=255) aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccccc > 1:1017711/1‑54 (MQ=255) aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccccc < 1:1231749/54‑1 (MQ=38) acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt < 1:749620/71‑1 (MQ=38) aaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa > 1:238799/1‑71 (MQ=35) tAATCCCATCTGGGCACATCCGATGGCAAGAGGCACGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCg > 1:1615283/1‑71 (MQ=21) aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg > 1:770310/1‑71 (MQ=34) aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg > 1:216827/1‑70 (MQ=35) aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg > 1:1388601/1‑70 (MQ=35) aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCg < 1:2213894/69‑1 (MQ=35) tCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa > 1:1964700/1‑71 (MQ=34) ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt > 1:1081209/1‑70 (MQ=35) ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt > 1:3057453/1‑70 (MQ=35) ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt < 1:2842782/70‑1 (MQ=35) ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa > 1:2063828/1‑65 (MQ=35) ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc < 1:1391161/67‑1 (MQ=35) ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc < 1:800696/67‑1 (MQ=35) ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAc < 1:1505412/70‑1 (MQ=35) cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt < 1:3786492/61‑1 (MQ=35) cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGt > 1:978229/1‑71 (MQ=34) acaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGtt > 1:2112972/1‑71 (MQ=34) aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc > 1:2906346/1‑71 (MQ=34) gATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt < 1:3647074/71‑1 (MQ=34) gCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtagt > 1:2401790/1‑70 (MQ=35) gCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt < 1:423333/71‑1 (MQ=34) aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTa > 1:461235/1‑70 (MQ=35) aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt < 1:3817253/71‑1 (MQ=34) agagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATc > 1:2166823/1‑71 (MQ=34) gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc < 1:751958/50‑1 (MQ=38) gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcc > 1:3015907/1‑71 (MQ=34) agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc < 1:1878144/49‑1 (MQ=38) gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa > 1:1178584/1‑71 (MQ=34) gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa > 1:372999/1‑71 (MQ=34) aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAt < 1:512645/71‑1 (MQ=34) aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAt < 1:1782649/71‑1 (MQ=34) aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAt > 1:3008085/1‑70 (MQ=35) aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATc < 1:600604/71‑1 (MQ=34) | GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATC > minE/1763567‑1763706 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |