| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,434,291 | T→C | 6.8% | intergenic (+43/‑93) | rhlB → / → gpp | ATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,434,291 | 0 | T | C | 6.8% | 165.9 / 2.4 | 59 | intergenic (+43/‑93) | rhlB/gpp | ATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase |
| Reads supporting (aligned to +/- strand): ref base T (37/18); new base C (2/2); total (39/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.98e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01 | |||||||||||
CCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACA > minE/2434222‑2434359 | ccGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATg > 1:1212777/1‑71 (MQ=255)ccGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTCTATCAAGATg > 1:3154847/1‑71 (MQ=255) aTCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCaa > 1:2752797/1‑68 (MQ=255) aTCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCaa > 1:4905170/1‑68 (MQ=255) tcgtcgtTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatg < 1:2014098/71‑1 (MQ=255) ttCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatga < 1:1491227/66‑1 (MQ=255) ttCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatga < 1:1541547/66‑1 (MQ=255) ttCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATccc > 1:2637809/1‑71 (MQ=255) tCAGGTTAAATAAAAAATGTCATACTCAAAATAATTCAAACTGTATCAAGATGGCAAGAgaatga < 1:4206878/65‑1 (MQ=255) gTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaat < 1:1229138/59‑1 (MQ=255) ttAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaa > 1:4194020/1‑57 (MQ=255) taaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGc > 1:4197692/1‑71 (MQ=255) taaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGc > 1:3528107/1‑71 (MQ=255) taaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCACGATGAGc > 1:5042383/1‑71 (MQ=255) aaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatgaat > 1:1404664/1‑60 (MQ=255) aaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCt > 1:4493462/1‑71 (MQ=255) ataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTa < 1:2318112/71‑1 (MQ=255) aaaaaTGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTg > 1:2179355/1‑71 (MQ=255) aTGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag > 1:2023369/1‑70 (MQ=255) tCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtagt < 1:1230803/68‑1 (MQ=255) taTACACAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCt < 1:997666/57‑1 (MQ=9) taTACACAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGt > 1:4293028/2‑71 (MQ=9) aTACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaa < 1:32071/41‑1 (MQ=255) aTACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGt > 1:5601956/1‑70 (MQ=255) tACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGa > 1:195662/1‑71 (MQ=255) tACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGa > 1:3953088/1‑71 (MQ=255) tACACAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGa < 1:4024673/71‑1 (MQ=9) cTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAg > 1:3876870/1‑52 (MQ=255) acAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGATTc > 1:2925141/2‑71 (MQ=12) cAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatgaa > 1:5197303/1‑40 (MQ=255) aaaaTCATTCAAACTGTATCAAGATGGCAAGAgaa > 1:5102755/1‑35 (MQ=255) aaaaTCATTCAAACTGTATCAAGATGGCAAGAgaa > 1:3690673/1‑35 (MQ=255) aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa > 1:3893752/1‑70 (MQ=255) aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa < 1:4193246/70‑1 (MQ=255) aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg > 1:5458564/1‑71 (MQ=255) aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg > 1:2517386/1‑71 (MQ=255) aaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCg < 1:3036918/43‑1 (MQ=255) aaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCt < 1:1694394/50‑1 (MQ=255) aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGAt > 1:4469841/1‑43 (MQ=255) aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGAtt > 1:4116376/1‑65 (MQ=255) aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt > 1:5154368/1‑70 (MQ=255) aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt > 1:192398/1‑70 (MQ=255) aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTg > 1:2951237/1‑59 (MQ=255) aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTg > 1:5533870/1‑59 (MQ=255) aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCACGATGAGCTTACTGTAGTAAGTg > 1:3781978/1‑59 (MQ=255) cAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag > 1:5087141/1‑50 (MQ=255) cAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag > 1:1442162/1‑50 (MQ=255) aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc < 1:4275211/58‑1 (MQ=255) aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc < 1:5226499/58‑1 (MQ=255) aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc < 1:2491927/58‑1 (MQ=255) tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGAtt < 1:55517/52‑1 (MQ=255) tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg < 1:1632874/69‑1 (MQ=255) tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg < 1:3800996/69‑1 (MQ=255) tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGc > 1:1956882/1‑70 (MQ=255) tCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCaa > 1:398427/1‑71 (MQ=255) aaGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAaca > 1:2247963/1‑71 (MQ=255) aaGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAac < 1:80674/70‑1 (MQ=255) aGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAaca > 1:4882327/1‑70 (MQ=255) gATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGaa > 1:3420785/1‑58 (MQ=255) tGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAAcaca > 1:4886431/1‑69 (MQ=255) | CCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACA > minE/2434222‑2434359 |
| 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: ‑ |