Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,434,301 | A→G | 9.8% | intergenic (+53/‑83) | 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,301 | 0 | A | G | 9.8% | 135.4 / 5.3 | 51 | intergenic (+53/‑83) | rhlB/gpp | ATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase |
Reads supporting (aligned to +/- strand): ref base A (30/16); new base G (2/3); total (32/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.48e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01 |
TCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTG > minE/2434233‑2434370 | 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:3953088/1‑71 (MQ=255) tACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGa > 1:195662/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:4193246/70‑1 (MQ=255) aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa > 1:3893752/1‑70 (MQ=255) aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg > 1:2517386/1‑71 (MQ=255) aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg > 1:5458564/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:192398/1‑70 (MQ=255) aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt > 1:5154368/1‑70 (MQ=255) aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTg > 1:5533870/1‑59 (MQ=255) aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTg > 1:2951237/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:5226499/58‑1 (MQ=255) aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc < 1:2491927/58‑1 (MQ=255) aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc < 1:4275211/58‑1 (MQ=255) tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGAtt < 1:55517/52‑1 (MQ=255) tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg < 1:3800996/69‑1 (MQ=255) tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg < 1:1632874/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) cggCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCaga < 1:3762691/70‑1 (MQ=2) gCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAgaga < 1:769099/71‑1 (MQ=255) agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg < 1:4304180/71‑1 (MQ=255) agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg < 1:2893665/71‑1 (MQ=255) agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg < 1:5298159/71‑1 (MQ=255) gagaATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCaa > 1:279177/1‑59 (MQ=255) gaatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc > 1:4161099/1‑55 (MQ=255) aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg < 1:2961155/49‑1 (MQ=255) aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg < 1:4740567/49‑1 (MQ=255) aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGAtt > 1:4042398/1‑70 (MQ=255) aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg < 1:5089868/71‑1 (MQ=255) atgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg < 1:1837188/66‑1 (MQ=255) atgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg < 1:1370517/70‑1 (MQ=255) | TCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTG > minE/2434233‑2434370 |
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: ‑ |