Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 394,838 | G→A | 15.7% | P63P (CCG→CCA) | ahpF → | alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 394,838 | 0 | G | A | 15.7% | 121.7 / 13.9 | 51 | P63P (CCG→CCA) | ahpF | alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base G (22/21); new base A (6/2); total (28/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.69e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01 |
TGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGG > minE/394768‑394900 | tGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCg > 1:1171283/1‑71 (MQ=255) ggCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCg < 1:735979/70‑1 (MQ=255) cTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCt > 1:2543070/1‑71 (MQ=255) cTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCt > 1:1098656/1‑71 (MQ=255) tGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCtt > 1:3189008/1‑71 (MQ=255) aaaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCtt < 1:1180880/69‑1 (MQ=255) aaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTcc < 1:931339/71‑1 (MQ=255) aaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTcc < 1:2142631/71‑1 (MQ=255) aGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCa < 1:2176445/71‑1 (MQ=255) aaCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGAt < 1:353972/67‑1 (MQ=255) aaCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCAc > 1:1416119/1‑70 (MQ=255) aCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCa < 1:3943057/71‑1 (MQ=255) cTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGAt > 1:2009776/1‑65 (MQ=255) cTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGAt > 1:2130067/1‑65 (MQ=255) gACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGt > 1:10810/1‑49 (MQ=255) aCAAATTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGt > 1:1651079/1‑71 (MQ=255) aCAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGt > 1:2574983/1‑71 (MQ=255) cAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGtt > 1:1352170/1‑71 (MQ=255) aaaGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAAccc < 1:299581/65‑1 (MQ=255) aaaGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTc < 1:3524168/71‑1 (MQ=255) aGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCt > 1:3827875/1‑51 (MQ=255) ccTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCa < 1:3951080/70‑1 (MQ=255) ccTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAg < 1:507782/71‑1 (MQ=255) cTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTg > 1:1908924/1‑46 (MQ=255) cTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAg < 1:1255016/70‑1 (MQ=255) ttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCaa > 1:2816917/1‑52 (MQ=255) ttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAgggg > 1:2664229/1‑71 (MQ=255) tAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGc > 1:201199/1‑71 (MQ=255) gaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAAc < 1:686568/57‑1 (MQ=255) aTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGttt < 1:2591426/70‑1 (MQ=255) tAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGTGCCACGTTTTg > 1:2968532/1‑71 (MQ=255) tAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg < 1:3230035/71‑1 (MQ=255) aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCAcc < 1:1479219/39‑1 (MQ=255) aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACg < 1:2657327/65‑1 (MQ=255) aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg > 1:2758197/1‑70 (MQ=255) aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg > 1:1291326/1‑70 (MQ=255) aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg > 1:3490388/1‑70 (MQ=255) aCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGc < 1:1776191/70‑1 (MQ=255) aCAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGc < 1:2854796/70‑1 (MQ=255) cAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCa > 1:3297002/1‑54 (MQ=255) cAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAgg > 1:1136182/1‑56 (MQ=255) cAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGCGCCACGTTTTGCa > 1:925430/1‑70 (MQ=255) cTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCt > 1:3939194/1‑71 (MQ=255) gCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAgg < 1:2464742/50‑1 (MQ=255) gTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCt < 1:472128/71‑1 (MQ=255) gCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTgg > 1:2367742/1‑71 (MQ=255) cGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTggg > 1:3784870/1‑71 (MQ=255) cGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTgg > 1:3280324/1‑70 (MQ=255) cGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTgg < 1:734275/70‑1 (MQ=255) tAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAg < 1:697844/40‑1 (MQ=255) aGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTg > 1:2270697/1‑65 (MQ=255) aGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTg > 1:319668/1‑65 (MQ=255) cGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTcccc < 1:2045783/58‑1 (MQ=255) gTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg < 1:181817/47‑1 (MQ=255) | TGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGG > minE/394768‑394900 |
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: ‑ |