Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 835,792 | C→G | 8.1% | S170R (AGC→AGG) | narH → | nitrate reductase 1, beta (Fe‑S) subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 835,792 | 0 | C | G | 8.1% | 173.0 / 3.9 | 61 | S170R (AGC→AGG) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
Reads supporting (aligned to +/- strand): ref base C (28/26); new base G (3/2); total (32/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGT > minE/835722‑835860 | aCCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGc < 1:684421/71‑1 (MQ=255) gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa < 1:9350/71‑1 (MQ=255) gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGACAGTTCGaa < 1:2735343/71‑1 (MQ=255) gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa > 1:1350532/1‑71 (MQ=255) aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa < 1:3080244/70‑1 (MQ=255) aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa > 1:2475322/1‑70 (MQ=255) aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa < 1:1817962/70‑1 (MQ=255) aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAac < 1:828789/71‑1 (MQ=255) aCAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa < 1:1257208/62‑1 (MQ=255) aCAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCtt < 1:3736150/70‑1 (MQ=255) cAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCa > 1:3185386/1‑71 (MQ=255) cTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatgat > 1:2773951/1‑71 (MQ=255) gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa > 1:1860346/1‑56 (MQ=255) gAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAtt < 1:3874513/71‑1 (MQ=255) aaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt < 1:1911298/70‑1 (MQ=255) aaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGACAGTTCGAAAACACCTTCATGATGt > 1:2040788/1‑66 (MQ=255) aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCa > 1:482701/1‑42 (MQ=255) aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAg < 1:3372591/43‑1 (MQ=255) aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc < 1:3845997/71‑1 (MQ=255) aCAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGtt > 1:3337209/1‑43 (MQ=255) cAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc > 1:2973612/1‑71 (MQ=255) aaGAACTTCGACAACATCCAGAAGGCGATGTATAGCCa > 1:3595282/1‑38 (MQ=255) aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc > 1:2139478/1‑66 (MQ=255) cTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtgt < 1:2042434/71‑1 (MQ=255) cTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtgt > 1:2842541/1‑71 (MQ=255) ttCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTc < 1:1021256/48‑1 (MQ=255) tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc > 1:2295280/1‑71 (MQ=255) cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatga > 1:2748015/1‑50 (MQ=255) gACAACATCCAGAAGGCGATGTATAGGCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg < 1:3809266/70‑1 (MQ=255) gACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg < 1:327566/70‑1 (MQ=255) caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat < 1:549767/45‑1 (MQ=255) caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat < 1:1667375/45‑1 (MQ=255) aacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGcc < 1:668319/61‑1 (MQ=255) caTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatga > 1:2481738/1‑44 (MQ=255) caTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatga > 1:1318590/1‑44 (MQ=255) caTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatga > 1:8274/1‑44 (MQ=255) aTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGt < 1:220207/46‑1 (MQ=255) tCCAGAAGGCGATGTATAGCCAGTTCGAAAACAcc < 1:1186687/35‑1 (MQ=255) ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAacac < 1:48615/67‑1 (MQ=255) caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt > 1:1609970/4‑71 (MQ=255) caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt > 1:3979858/4‑71 (MQ=255) caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt > 1:1437971/4‑71 (MQ=255) cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatg > 1:448204/1‑39 (MQ=255) cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAaca > 1:1616683/1‑65 (MQ=255) aaGGAGATGTACGGGCAGTTTGAAAACACCTTC‑‑‑ATGTACCTGCCGCGCCTGTGCGAac < 1:801635/58‑1 (MQ=255) aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAaca > 1:3820804/1‑61 (MQ=255) ggCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc < 1:3816693/71‑1 (MQ=255) gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTTCGAac > 1:1861873/1‑58 (MQ=255) gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGACCACTGCCTGAACCCg > 1:3526373/1‑70 (MQ=255) gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCg > 1:2089863/1‑70 (MQ=255) aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGcc < 1:2474934/46‑1 (MQ=255) aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGcc < 1:2674594/46‑1 (MQ=255) aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac < 1:1813714/55‑1 (MQ=255) aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg > 1:2079504/1‑70 (MQ=255) tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa > 1:2888282/1‑51 (MQ=255) tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtg < 1:1048690/71‑1 (MQ=255) ataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc > 1:3322785/1‑47 (MQ=255) ataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtg < 1:3191144/70‑1 (MQ=255) taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtgt < 1:2586662/44‑1 (MQ=255) taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa > 1:422970/1‑49 (MQ=255) taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg > 1:1669244/1‑71 (MQ=255) gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt > 1:2683894/1‑70 (MQ=255) | ACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGT > minE/835722‑835860 |
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: ‑ |