| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 833,885 | C→G | 14.3% | L781L (CTC→CTG) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 833,885 | 0 | C | G | 14.3% | 137.2 / 13.8 | 56 | L781L (CTC→CTG) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base C (27/21); new base G (4/4); total (31/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01 | |||||||||||
GACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCA > minE/833817‑833953 | gACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa > 1:1941801/1‑71 (MQ=255) cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc > 1:3766115/1‑67 (MQ=255) gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa > 1:628468/1‑70 (MQ=255) gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa > 1:517880/1‑70 (MQ=255) tttgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGa < 1:1693588/65‑1 (MQ=255) attgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACATGCCTGTTCTCCGatat < 1:78509/68‑1 (MQ=255) tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat > 1:3906359/1‑70 (MQ=255) tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat > 1:1794463/1‑70 (MQ=255) gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca < 1:2543467/70‑1 (MQ=255) gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcataat > 1:1805447/1‑71 (MQ=255) aaGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTc > 1:3540995/1‑61 (MQ=255) aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt < 1:2103152/71‑1 (MQ=255) aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt < 1:2738845/71‑1 (MQ=255) aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtt > 1:2705971/1‑70 (MQ=255) ggCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtttt < 1:3019309/71‑1 (MQ=255) gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc < 1:142174/53‑1 (MQ=255) gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg < 1:2905222/60‑1 (MQ=255) gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTg > 1:465022/1‑71 (MQ=255) cAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa > 1:337403/1‑54 (MQ=255) aGCTGGATCTGGTGGTTACCCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa > 1:2885293/1‑52 (MQ=255) gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGa > 1:996399/1‑71 (MQ=255) cTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGAc < 1:2875369/71‑1 (MQ=255) tCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTAt < 1:476899/46‑1 (MQ=255) tCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGa > 1:1051907/1‑71 (MQ=255) cTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc > 1:1257131/1‑71 (MQ=255) tggtggTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc < 1:2172931/41‑1 (MQ=255) tggtggTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGAc > 1:2466527/1‑64 (MQ=255) gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt < 1:3694872/70‑1 (MQ=255) gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg > 1:1699652/1‑71 (MQ=255) tggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg < 1:2478699/71‑1 (MQ=255) ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg > 1:2964215/1‑70 (MQ=255) ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg > 1:2571116/1‑70 (MQ=255) gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg > 1:871411/1‑43 (MQ=255) gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg > 1:2511607/1‑43 (MQ=255) gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGt < 1:1800876/70‑1 (MQ=255) aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTAt > 1:260950/1‑36 (MQ=255) aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg > 1:1969392/1‑65 (MQ=255) aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg > 1:2002814/1‑65 (MQ=255) cGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaa < 1:1824569/71‑1 (MQ=255) gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAc < 1:704334/40‑1 (MQ=255) gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa < 1:1861217/71‑1 (MQ=255) gCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa < 1:2948002/71‑1 (MQ=255) tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa > 1:3614129/3‑69 (MQ=255) tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa > 1:1667589/3‑68 (MQ=255) tGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc > 1:1817918/1‑58 (MQ=255) cgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGACACCTGGTACGaaa > 1:2271443/2‑68 (MQ=255) gACTTCCGTCTGTCGAGCACGTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaat < 1:1768595/70‑2 (MQ=255) gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga > 1:1895780/1‑71 (MQ=255) ttCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca < 1:3573381/37‑1 (MQ=255) ttCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca < 1:3568828/37‑1 (MQ=255) ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt < 1:1625578/55‑1 (MQ=255) ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt < 1:2947024/55‑1 (MQ=255) tGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa < 1:2073247/57‑1 (MQ=255) tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGaa < 1:3617365/68‑1 (MQ=255) tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATAc > 1:1785461/1‑71 (MQ=255) gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACtt > 1:1507433/1‑71 (MQ=255) aGCACCTGTCTCTATTCCGACATCATTTTGCCGACgg > 1:2317658/1‑37 (MQ=255) tGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATg < 1:2630184/57‑1 (MQ=255) gTCTCTATTCCGACATCATTTTGCCGACGGCGAcc < 1:1438584/35‑1 (MQ=255) actctATTCCGACATCATTTTGCCGACGGCTACCTg > 1:3879545/2‑36 (MQ=38) ctctATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCa > 1:1652325/1‑71 (MQ=255) | GACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCA > minE/833817‑833953 |
| 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: ‑ |