Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 833,966 | C→T | 48.7% | H808H (CAC→CAT) | narG → | nitrate reductase 1, alpha subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 833,966 | 0 | C | T | 48.7% | 5.3 / 105.2 | 80 | H808H (CAC→CAT) | narG | nitrate reductase 1, alpha subunit |
Reads supporting (aligned to +/- strand): ref base C (12/29); new base T (16/23); total (28/52) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.50e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01 |
CATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCA > minE/833897‑834030 | catTTTGCCGACGGCGACCTGGTACGAAAAAGAAGACATGAATACTTCGGATATGCATCCGTTTATTCAcc > 1:2473935/1‑71 (MQ=255) caatttGCCGACGGCGACCTGGTACGAAAAAGACTACATGAATACTTCGGATATGCATCCGTTTATTCAcc > 1:1915154/4‑71 (MQ=255) atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAccc > 1:1279436/1‑71 (MQ=255) gCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg < 1:3934286/70‑1 (MQ=255) cGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcg > 1:182574/1‑71 (MQ=255) gCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc < 1:4778080/71‑1 (MQ=255) gCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc < 1:2071345/71‑1 (MQ=255) cGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg < 1:1610348/62‑1 (MQ=255) cTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTc > 1:3764463/1‑70 (MQ=255) cTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCg < 1:495935/71‑1 (MQ=255) gTACGAAAAAGACTATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa < 1:2753154/69‑1 (MQ=255) gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCttt > 1:2928721/1‑54 (MQ=255) gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa < 1:774894/69‑1 (MQ=255) gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATc > 1:5188932/1‑71 (MQ=255) gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCTGTCGAt < 1:5336978/70‑1 (MQ=255) gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt < 1:4590972/70‑1 (MQ=255) gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt < 1:2241915/70‑1 (MQ=255) gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt < 1:3132333/70‑1 (MQ=255) tACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCt < 1:1476946/53‑1 (MQ=255) aCGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCg > 1:93619/1‑71 (MQ=255) gAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAcc < 1:533759/46‑1 (MQ=255) aaaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGcc > 1:159557/1‑71 (MQ=255) aaaaGACGATATGAACATCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGcgt < 1:5196116/71‑3 (MQ=255) aaaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGcgt < 1:1165135/71‑3 (MQ=255) aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGcc < 1:2971644/70‑1 (MQ=255) aaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTg < 1:527492/71‑1 (MQ=255) acgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc < 1:4214709/52‑1 (MQ=255) cgatATGAACACCTCGGATATGCATCCGTTTATTCATCCGCttt > 1:796807/1‑42 (MQ=25) cgatATGAACACCTCGGATATGCATCCGTTTATTCATCCGCttt > 1:1265174/1‑42 (MQ=25) cgatATGAACACCTCGGATATGCATCCGTTTATTCATCCGCttt > 1:2625660/1‑42 (MQ=25) gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg < 1:2757281/70‑1 (MQ=255) gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg < 1:673802/70‑1 (MQ=255) gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg < 1:5605469/70‑1 (MQ=255) ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg < 1:4542792/46‑1 (MQ=255) ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg < 1:3226024/46‑1 (MQ=255) ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg < 1:4017036/46‑1 (MQ=255) ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg < 1:3953829/46‑1 (MQ=255) acATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa < 1:295201/55‑1 (MQ=255) cATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg > 1:4451928/1‑40 (MQ=255) aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcggcg > 1:1377363/1‑48 (MQ=255) aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcggcg > 1:761349/1‑48 (MQ=255) aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcagcg > 1:1760693/1‑48 (MQ=255) aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCg < 1:52347/52‑1 (MQ=37) aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa < 1:3625764/53‑1 (MQ=37) tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATc < 1:2116227/54‑1 (MQ=255) tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa < 1:2070046/69‑1 (MQ=255) gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGaa > 1:4454508/1‑69 (MQ=255) aaTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg > 1:1726496/1‑64 (MQ=255) aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:3035158/71‑1 (MQ=255) cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa < 1:2661546/46‑1 (MQ=255) cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa < 1:317352/46‑1 (MQ=255) cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa < 1:642867/46‑1 (MQ=255) ctCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGaa < 1:2503919/57‑1 (MQ=255) ctCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCAGCGTGGGaa < 1:2231905/57‑1 (MQ=255) tCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATc > 1:3815887/1‑46 (MQ=255) cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:5187583/1‑71 (MQ=255) ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTggg < 1:1846389/71‑1 (MQ=255) ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:3689816/1‑70 (MQ=255) tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGTGaaa < 1:1481207/71‑1 (MQ=255) tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGaaa < 1:2915590/71‑1 (MQ=255) tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGa < 1:3322576/69‑1 (MQ=255) tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCgg < 1:3744128/42‑1 (MQ=38) tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg < 1:2517544/48‑1 (MQ=37) tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg < 1:4323523/48‑1 (MQ=37) tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg < 1:1551834/48‑1 (MQ=37) tGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGc < 1:2124322/52‑1 (MQ=255) gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGa > 1:2439283/1‑60 (MQ=255) cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCt < 1:3401193/70‑1 (MQ=255) cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa < 1:4377260/71‑1 (MQ=255) cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa < 1:771430/71‑1 (MQ=255) tCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:78578/55‑1 (MQ=255) cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaaa < 1:1469837/71‑1 (MQ=255) cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACa > 1:4287315/1‑69 (MQ=255) cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt > 1:2225820/1‑39 (MQ=255) cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACa > 1:3522696/1‑69 (MQ=255) gTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAAt < 1:4112380/63‑1 (MQ=255) tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa < 1:2210751/46‑1 (MQ=255) tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:5515961/51‑1 (MQ=255) ttATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAgg > 1:4407565/1‑69 (MQ=255) tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:5665279/1‑55 (MQ=255) tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:5369519/1‑55 (MQ=255) tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggta > 1:4662936/1‑68 (MQ=255) tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa > 1:3288003/1‑71 (MQ=255) tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa > 1:5064259/1‑71 (MQ=255) | CATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCA > minE/833897‑834030 |
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: ‑ |