| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 834,800 | G→A | 9.8% | A1086A (GCG→GCA) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 834,800 | 0 | G | A | 9.8% | 217.2 / 11.0 | 71 | A1086A (GCG→GCA) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (24/40); new base A (4/3); total (28/43) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.22e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01 | |||||||||||
TCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGC > minE/834731‑834867 | tCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAgcgc > 1:4139931/1‑71 (MQ=255) aTCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACtt < 1:833509/71‑1 (MQ=255) tCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGTGCTCAACTTc < 1:1367901/71‑1 (MQ=255) gACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACtt < 1:602333/68‑1 (MQ=255) cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTgc > 1:2390213/1‑70 (MQ=255) cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGa > 1:1291459/1‑71 (MQ=255) cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGa > 1:3660710/1‑71 (MQ=255) cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGa > 1:3792761/1‑71 (MQ=255) cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGa > 1:565879/1‑71 (MQ=255) cccGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAgcg < 1:2735635/56‑1 (MQ=255) ccGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGAcgc < 1:503402/71‑1 (MQ=255) cGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt < 1:1703389/65‑1 (MQ=255) ttCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACt < 1:2147493/59‑1 (MQ=255) cGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCAc > 1:2487440/1‑71 (MQ=255) tGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTc < 1:1608274/50‑1 (MQ=255) tGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAg < 1:2722209/71‑1 (MQ=255) tGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAg < 1:2717274/71‑1 (MQ=255) gAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAgcgc > 1:1303748/1‑47 (MQ=255) aaaGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACtt < 1:160227/53‑1 (MQ=255) aaaGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAg > 1:3056317/1‑69 (MQ=255) aaaGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa < 1:2131026/70‑1 (MQ=255) agaagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTc > 1:1664772/1‑52 (MQ=255) agaagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTc > 1:4116999/1‑52 (MQ=255) agaagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTc > 1:2259209/1‑52 (MQ=255) aagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGAcgccgc < 1:839932/60‑1 (MQ=255) aagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGAcgc > 1:3208786/1‑57 (MQ=255) aagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGAcg > 1:2007626/1‑56 (MQ=255) agTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTggg < 1:186127/71‑1 (MQ=255) aTAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt < 1:2478335/47‑1 (MQ=255) aTAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt < 1:2482981/47‑1 (MQ=255) aGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACt < 1:2603046/41‑1 (MQ=255) aGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAt > 1:1689175/1‑69 (MQ=255) ggCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt < 1:1525185/44‑1 (MQ=255) gCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccac < 1:182999/71‑1 (MQ=255) ccAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa > 1:1519947/1‑56 (MQ=255) cAGAAATCCAACGGCAACCAGGAAAAAGCGCTCaa > 1:3092901/1‑35 (MQ=255) cAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTc < 1:912205/39‑1 (MQ=255) aGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt < 1:4116763/40‑1 (MQ=255) gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt < 1:3454035/39‑1 (MQ=255) gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTggg < 1:2273653/59‑1 (MQ=255) gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc < 1:2474456/64‑1 (MQ=255) gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc < 1:3055214/64‑1 (MQ=255) aaaTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCCGAAGTGGGGTAtccactcca < 1:3052458/70‑1 (MQ=255) aaTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCAcc > 1:1555038/1‑71 (MQ=255) aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCAc < 1:1179471/46‑1 (MQ=255) aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAg < 1:2562332/52‑1 (MQ=255) aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAt < 1:1423086/60‑1 (MQ=255) aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCt < 1:583933/71‑1 (MQ=255) tCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa < 1:615640/49‑1 (MQ=255) tCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa < 1:990397/49‑1 (MQ=255) caaaCGGCTTCCCGGAAAAAGCACTTAACTTCCTGACGCCGCACCAGa < 1:1225982/46‑1 (MQ=19) cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGaa > 1:2073744/1‑48 (MQ=255) cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTa > 1:2180752/1‑69 (MQ=255) cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTAc < 1:1052559/70‑1 (MQ=255) aaCGGCTTCCCGGAAAAAGCACTTAACTTCCTGACGCCGCACCAGaaa > 1:3225533/1‑47 (MQ=21) aaCGGCTTCCCGGAAAAAGCACTTAACTTCCTGACGCCGCACCAGaaa > 1:3818476/1‑47 (MQ=21) aCGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcc > 1:988776/1‑57 (MQ=255) aCGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGc > 1:3779071/1‑71 (MQ=255) ggCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGa > 1:2348567/1‑71 (MQ=255) gCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAg < 1:2806355/44‑1 (MQ=255) gCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAg < 1:4141499/44‑1 (MQ=255) cAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa < 1:1441309/71‑1 (MQ=255) cccggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAcc > 1:1675368/4‑71 (MQ=255) aGGAAAAAGCGCTCAACTTCCTGACGCCGAACCAGAAGTgg > 1:2864673/1‑41 (MQ=255) ggAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa < 1:1840370/35‑1 (MQ=255) gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct < 1:1600553/71‑1 (MQ=255) aaaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa < 1:1109918/63‑1 (MQ=255) aaaaaGCACTTAACTTCCTGACGCCGCACCAGAAATggggca > 1:1583382/1‑40 (MQ=25) aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgctg > 1:2656707/1‑70 (MQ=255) aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctga > 1:3863982/1‑71 (MQ=255) aaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct < 1:1595360/68‑1 (MQ=255) aGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa < 1:1057953/59‑1 (MQ=255) aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc < 1:259185/71‑1 (MQ=255) aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc < 1:1257276/71‑1 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct < 1:2611229/65‑1 (MQ=255) | TCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGC > minE/834731‑834867 |
| 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: ‑ |