| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 4,604,824 | G→A | 57.1% | S16S (AGC→AGT) | yjjA ← | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,604,824 | 0 | G | A | 57.1% | 3.1 / 24.0 | 21 | S16S (AGC→AGT) | yjjA | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (4/5); new base A (4/8); total (8/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.73e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGA > W3110S.gb/4604774‑4604887 | gagaGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCACTGGCGGCAATGgcacagca > 1:965361/1‑71 (MQ=255) cATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGgcacagcaca > 1:262324/1‑66 (MQ=255) cATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCACTGGCGGCAATGGcacagcacag > 1:1371526/1‑67 (MQ=255) tCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCACTGGCGGCAATGGCACAGCACAGAAGATGttt > 1:1480719/1‑70 (MQ=255) tCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCACTGGCGGCAATGGCACAGCACAGAAGATGttt < 1:24731/70‑1 (MQ=255) tCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCACTGGCGGCAATGGCACAGCACAGAAGATGtt < 1:1054235/69‑1 (MQ=255) cAGGACGCAGCATGCACCCCGGTGGAAATTAATGCACTGGCGGCAATGGCACAGCACAGAAGATGTTTAAc < 1:207511/71‑1 (MQ=255) ggACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAg > 1:1291672/1‑71 (MQ=255) cGCAGCATGCACCCCGGTGGAATTTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGttt > 1:316096/1‑71 (MQ=255) cagcaTGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTtc < 1:1155583/71‑1 (MQ=255) ccccGGTGGAAATTAATGCACTGGCGGCAATGGCACAGCACagaa < 1:222068/45‑1 (MQ=255) ccccGGTGGAAATTAATGCACTGGCGGCAATGGCACAGCACagaa < 1:299288/45‑1 (MQ=255) cccGGTGGAAATTAATGCACTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTa < 1:621595/71‑1 (MQ=255) tGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTtc < 1:140358/56‑1 (MQ=255) ggAAATTAATGCGCGGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTtc < 1:474507/55‑1 (MQ=255) ggAAATTAATGCACTGGCGGCAATGGCACAGCACagaa > 1:725314/1‑38 (MQ=255) ttAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGa < 1:195441/71‑1 (MQ=255) ttAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGa < 1:261698/71‑1 (MQ=255) tAATGCACTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCa < 1:5158/63‑1 (MQ=255) tAATGCACTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCa < 1:1359482/63‑1 (MQ=255) aaTGCGCTGGCGGAAATGGCACAGCACAGAAGATGTTTAACAGTTTtcatcatca > 1:147946/1‑55 (MQ=255) | GAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGA > W3110S.gb/4604774‑4604887 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |