| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,591,637 | T→A | 17.8% | intergenic (+96/+93) | 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,637 | 0 | T | A | 17.8% | 181.2 / 30.5 | 73 | intergenic (+96/+93) | yfeS/cysM | conserved hypothetical protein/cysteine synthase B |
| Reads supporting (aligned to +/- strand): ref base T (24/36); new base A (8/5); total (32/41) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.19e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 | |||||||||||
TTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCA > minE/1591570‑1591706 | ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa < 1:2719/70‑1 (MQ=255)taagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa < 1:873801‑M1/46‑2 (MQ=255) aagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat > 1:2248164‑M1/25‑69 (MQ=255) aagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat > 1:2856372‑M1/25‑69 (MQ=255) aagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat > 1:1298180‑M1/25‑69 (MQ=255) aagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat < 1:3521988‑M1/47‑3 (MQ=255) aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat < 1:3225617‑M1/47‑3 (MQ=255) agcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACgctaaatc > 1:351329‑M1/24‑68 (MQ=255) agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat > 1:541617‑M1/24‑68 (MQ=255) agcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACgctaaatc > 1:800747‑M1/24‑68 (MQ=255) aaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa > 1:312931/1‑71 (MQ=35) gtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAAtcagg > 1:279019‑M1/21‑69 (MQ=255) tcgcatcaggcgttgatgtCGGATGCGGCATTAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc > 1:1712698‑M1/20‑68 (MQ=255) tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc < 1:120748/64‑1 (MQ=35) catcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa > 1:2714985‑M1/17‑65 (MQ=255) atcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa > 1:189275‑M1/16‑60 (MQ=255) gcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc < 1:2956338/71‑1 (MQ=39) caggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag < 1:1515700‑M1/50‑2 (MQ=255) caggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag < 1:573430‑M1/50‑2 (MQ=255) caggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag < 1:3023827‑M1/50‑2 (MQ=255) ttGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa < 1:1793520/60‑1 (MQ=35) tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa > 1:178762/1‑71 (MQ=11) gTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGaaa > 1:3629449/1‑71 (MQ=255) ttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat > 1:421014‑M1/8‑52 (MQ=255) ttATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAcc > 1:2415884/1‑71 (MQ=255) tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAcc < 1:2674011/70‑1 (MQ=255) atgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa > 1:2038921‑M1/5‑53 (MQ=255) tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTg < 1:3094598/71‑1 (MQ=255) gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTg > 1:368045/1‑70 (MQ=255) gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaatt > 1:7191/1‑60 (MQ=11) gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATAGTGCAAATTCAATACATTGCAGaatt > 1:521196/1‑60 (MQ=11) gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAgg > 1:914177/1‑70 (MQ=255) tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa > 1:2794054/1‑58 (MQ=11) 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) tAAACGCCTTATCCGGGCTACAAAATCGTGCTAATTCAATATATTGCAGAAAcc < 1:194189/54‑1 (MQ=255) tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGc < 1:3232102/62‑1 (MQ=11) tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGc < 1:1552453/71‑1 (MQ=11) tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAata < 1:1871677/41‑1 (MQ=255) tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTg < 1:1078776/46‑1 (MQ=255) tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAcc < 1:822651/54‑1 (MQ=255) tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCttc < 1:3376450/57‑2 (MQ=255) aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTc > 1:2590870/1‑35 (MQ=40) aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc < 1:1044093/44‑1 (MQ=255) aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTa < 1:767442/71‑1 (MQ=255) cGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGTAGAAATCATGTAGGCCTGATAAGCGTAg < 1:885947/71‑1 (MQ=11) ttATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGa < 1:3654295/42‑1 (MQ=255) 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:3311919/1‑70 (MQ=11) aTCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAGCGCAt < 1:3337832/70‑1 (MQ=11) aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGt < 1:796068/63‑1 (MQ=255) aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCa < 1:902692/69‑1 (MQ=255) aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt < 1:1139618/70‑1 (MQ=255) aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa < 1:3096282/40‑1 (MQ=11) tacGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa > 1:1756942/3‑40 (MQ=11) tCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:3479931/1‑69 (MQ=255) tCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCa < 1:2009982/71‑1 (MQ=255) tCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCa < 1:2904689/71‑1 (MQ=255) tCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa > 1:3711414/1‑40 (MQ=11) 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) | TTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCA > minE/1591570‑1591706 |
| 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: ‑ |