| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,763,636 | A→G | 17.1% | noncoding (208/1542 nt) | rrsG ← | 16S ribosomal RNA |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,763,636 | 0 | A | G | 17.1% | 74.8 / 12.0 | 35 | noncoding (208/1542 nt) | rrsG | 16S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base A (14/15); new base G (3/3); total (17/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01 | |||||||||||
TGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCA > minE/1763571‑1763704 | tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc > 1:1931139/1‑71 (MQ=255)tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc > 1:1404769/1‑71 (MQ=40) tGAGCCGTTACCCCACCAACGAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccccc < 1:2868434/71‑1 (MQ=21) aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCTATGGCAAGAGGCCCGAGGGTCCCCct < 1:488838/70‑1 (MQ=255) ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCttt < 1:80402/69‑1 (MQ=255) ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg < 1:1640189/71‑1 (MQ=39) aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc > 1:1053686/1‑71 (MQ=39) aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc > 1:1751742/1‑71 (MQ=39) ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc < 1:2776350/70‑1 (MQ=255) caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc > 1:805036/1‑71 (MQ=255) caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc > 1:2207905/1‑71 (MQ=255) caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc > 1:2733132/1‑71 (MQ=39) ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctct > 1:2034634/1‑59 (MQ=255) caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa < 1:531351/70‑1 (MQ=39) acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc > 1:1184266/1‑49 (MQ=38) ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa > 1:341529/1‑59 (MQ=31) ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGTGGTAtt < 1:384579/71‑1 (MQ=21) ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAtt < 1:825389/71‑1 (MQ=34) ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt < 1:656023/70‑1 (MQ=35) aTCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAg > 1:413244/1‑71 (MQ=34) ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCt < 1:2315421/68‑1 (MQ=35) gCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa < 1:240507/68‑1 (MQ=35) cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGt > 1:2952521/1‑71 (MQ=34) acaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGtt < 1:73962/71‑1 (MQ=34) acaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGtt < 1:1089689/71‑1 (MQ=34) aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc > 1:2709480/1‑71 (MQ=34) cGATGGCAAGAGGCCCGAATGTCCCCCTCTTTGGTCTTGCGACATTATGCgg > 1:2334779/1‑52 (MQ=255) tGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt < 1:292929/52‑1 (MQ=35) ggCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtagt > 1:1993495/1‑71 (MQ=34) agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt > 1:2439691/1‑45 (MQ=38) agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGtt < 1:409788/56‑1 (MQ=35) gCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATCCGGTATTAGCTACCGTTTcc < 1:2692995/57‑1 (MQ=255) cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTcc < 1:145816/71‑1 (MQ=34) cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc > 1:306626/1‑70 (MQ=35) gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa < 1:885129/71‑1 (MQ=34) aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtagt > 1:2646338/1‑58 (MQ=35) | TGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCA > minE/1763571‑1763704 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |