| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,591,660 | T→A | 10.0% | intergenic (+119/+70) | yfeS → / ← cysM | conserved hypothetical protein/cysteine synthase B |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,591,660 | 0 | T | A | 10.0% | 158.4 / 9.7 | 60 | intergenic (+119/+70) | yfeS/cysM | conserved hypothetical protein/cysteine synthase B |
| Reads supporting (aligned to +/- strand): ref base T (23/30); new base A (2/4); total (25/35) | |||||||||||
| 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.71e-01 | |||||||||||
TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCC > minE/1591592‑1591729 | tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTg < 1:3094598/71‑1 (MQ=255) gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTg > 1:368045/1‑70 (MQ=255) gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAgg > 1:914177/1‑70 (MQ=255) tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa > 1:3191892/1‑66 (MQ=255) cggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGa > 1:1079344/1‑71 (MQ=255) ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGAt < 1:1043999/71‑1 (MQ=255) ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGAt < 1:3787615/71‑1 (MQ=255) ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGATAATTCAATATATTGCAGAAACCTTGTAGGCCTGa > 1:1335392/1‑70 (MQ=255) cgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCtt < 1:902035/58‑1 (MQ=255) tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGc < 1:3232102/62‑1 (MQ=11) tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGc < 1:1552453/71‑1 (MQ=11) tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCttc < 1:3376450/57‑2 (MQ=255) aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTa < 1:767442/71‑1 (MQ=255) cGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGTAGAAATCATGTAGGCCTGATAAGCGTAg < 1:885947/71‑1 (MQ=11) tATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAGCGCAt > 1:1982966/1‑71 (MQ=11) tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa > 1:2455354/1‑50 (MQ=255) tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:3651780/1‑71 (MQ=255) tATCCGGCCTACAAAATCGTGCTAATTCAATATATT‑CAGAAACCCTGt < 1:2019118/48‑1 (MQ=17) aTCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAGCGCAt < 1:3337832/70‑1 (MQ=11) aTCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAGCGCAt > 1:3311919/1‑70 (MQ=11) aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGt < 1:796068/63‑1 (MQ=255) aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCa < 1:902692/69‑1 (MQ=255) aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt < 1:1139618/70‑1 (MQ=255) tCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:3479931/1‑69 (MQ=255) tCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCa < 1:2904689/71‑1 (MQ=255) tCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCa < 1:2009982/71‑1 (MQ=255) cTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGt > 1:2094927/1‑71 (MQ=255) tACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATa < 1:330227/50‑1 (MQ=255) tACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAg < 1:3370196/65‑1 (MQ=255) aCAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa < 1:1582663/40‑1 (MQ=255) aCAAAATCGTGCTAACTCAATATATTGCAGAAACCTTg < 1:3494534/38‑1 (MQ=255) cAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt > 1:2725506/1‑70 (MQ=255) cAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtt < 1:3343955/69‑1 (MQ=255) aaaaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTa > 1:2587683/1‑53 (MQ=255) aaaaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa > 1:3533124/1‑65 (MQ=255) aaaaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtttt > 1:2226687/1‑70 (MQ=255) aaaaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTg > 1:504686/1‑71 (MQ=255) aaaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGc < 1:2541037/71‑1 (MQ=255) aaaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGc < 1:3787753/71‑1 (MQ=255) aTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGt < 1:2231481/71‑1 (MQ=255) tGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt < 1:1944330/71‑1 (MQ=255) cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca > 1:1131552/1‑71 (MQ=255) cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca > 1:1066242/1‑71 (MQ=255) cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtc < 1:1575277/70‑1 (MQ=255) aaTTCAATATATTGCAGAAACCTTGTAGGCCTGATaa < 1:932603/37‑1 (MQ=255) aaTTCAATATATTGCAGAAACCTTGTAGGCCTGATaa < 1:2299356/37‑1 (MQ=255) aaTTCAATATATTGCAGAAACCTTGTAGGCCTGATaa < 1:2275496/37‑1 (MQ=255) ttCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTa < 1:3204818/40‑1 (MQ=255) ttCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTa < 1:3918802/40‑1 (MQ=255) ttCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAg > 1:3492262/1‑71 (MQ=255) cAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTa > 1:3859369/1‑38 (MQ=255) cAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGttt > 1:852190/1‑61 (MQ=255) aaTATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGtc < 1:2035183/70‑1 (MQ=255) atatatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATc > 1:481255/1‑43 (MQ=255) atatatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGtctc > 1:170672/1‑71 (MQ=255) atatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGttt < 1:2513640/57‑1 (MQ=255) tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGa < 1:992398/71‑1 (MQ=255) tGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGtc < 1:1367725/63‑1 (MQ=255) cAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCt < 1:138625/70‑1 (MQ=255) gAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGtctc > 1:3303897/1‑61 (MQ=40) aCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAAt > 1:2005530/1‑71 (MQ=255) cTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATcc < 1:1053587/71‑1 (MQ=255) | TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCC > minE/1591592‑1591729 |
| 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: ‑ |