| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 832,115 | A→G | 18.7% | P191P (CCA→CCG) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 832,115 | 0 | A | G | 18.7% | 170.3 / 33.2 | 80 | P191P (CCA→CCG) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base A (37/28); new base G (4/11); total (41/39) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.60e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01 | |||||||||||
TGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCT > minE/832048‑832185 | tGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAAtt > 1:1041434/1‑71 (MQ=255) tgccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCgg > 1:935703/2‑71 (MQ=255) cgccgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCg > 1:1298201/1‑70 (MQ=255) cgccgcATCTAACGTTTACACCATCAAAAACTACGGCCCCGACCGTGTTGCTGGTTTCTCGCCAATTCCgg > 1:3104682/1‑71 (MQ=255) gccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCgg < 1:3031776/70‑1 (MQ=255) gcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAAt < 1:2572548/71‑1 (MQ=255) gcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAAt < 1:2560064/71‑1 (MQ=255) cATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATg > 1:3678561/1‑71 (MQ=255) cgcTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGt < 1:3671646/69‑1 (MQ=255) cgcTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGt < 1:636593/69‑1 (MQ=255) cTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc > 1:3169106/1‑70 (MQ=255) aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCg < 1:3300731/68‑1 (MQ=255) ttACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt > 1:1988771/1‑70 (MQ=255) tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGttt < 1:866329/70‑1 (MQ=255) tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTc > 1:3519133/1‑71 (MQ=255) gacCATCAAAACTTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGttt < 1:2066270/67‑1 (MQ=255) gacCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCtt < 1:1990199/70‑1 (MQ=255) cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCtt < 1:1292348/71‑1 (MQ=255) acCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCt < 1:3504019/69‑1 (MQ=255) acCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTa > 1:3564092/1‑71 (MQ=255) ccATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTAc < 1:3482110/71‑1 (MQ=255) cATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACg < 1:2068268/71‑1 (MQ=255) aaaacCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGccg < 1:757389/70‑3 (MQ=255) aaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa > 1:2730605/1‑47 (MQ=255) aaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa > 1:2802201/1‑47 (MQ=255) aaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCgg < 1:3312136/71‑1 (MQ=255) aaaCTACGGCCCGGACCGTGTTGCTGGTTTATCGCCAATTCCGGCaa > 1:942336/1‑47 (MQ=255) aaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCggg < 1:3059410/71‑1 (MQ=255) aaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCgg > 1:371405/1‑70 (MQ=255) cTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCtt < 1:2516922/59‑1 (MQ=255) aCGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTa > 1:204929/1‑58 (MQ=255) cGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGtt > 1:1641533/1‑52 (MQ=255) cGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCa < 1:1010444/70‑1 (MQ=255) ggCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGa > 1:3263008/1‑46 (MQ=255) ggCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGa > 1:2791558/1‑46 (MQ=255) gCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt < 1:2548387/50‑1 (MQ=255) ccGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTa > 1:436276/1‑71 (MQ=255) ggACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtc > 1:3637623/1‑71 (MQ=255) ggACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtc > 1:616363/1‑71 (MQ=255) gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc < 1:2185862/38‑1 (MQ=255) gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATgg > 1:1314347/1‑43 (MQ=255) gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATg < 1:954096/42‑1 (MQ=255) gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATg < 1:1062666/42‑1 (MQ=255) gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATg < 1:2907780/42‑1 (MQ=255) gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct < 1:1695461/71‑1 (MQ=255) gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct < 1:1066680/71‑1 (MQ=255) aCCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctc < 1:14805/71‑1 (MQ=255) cGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctctc < 1:1882145/71‑1 (MQ=255) gtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAt > 1:2504553/1‑63 (MQ=255) ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt > 1:1370585/1‑70 (MQ=255) ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt > 1:3640714/1‑70 (MQ=255) ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt > 1:1338445/1‑70 (MQ=255) ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCAAGCTATCTCTCGCt > 1:3011493/1‑70 (MQ=255) tGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTgc > 1:1516964/1‑70 (MQ=255) gCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt > 1:3524794/1‑71 (MQ=255) ccggTTTCTCGCCGATCCCGGCTATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCtt < 1:1692044/69‑1 (MQ=255) cTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAtt < 1:276185/71‑1 (MQ=255) cTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt < 1:3205353/70‑1 (MQ=255) cggTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTg < 1:1763539/70‑1 (MQ=255) gTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGc < 1:1921210/71‑1 (MQ=255) gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa > 1:688888/1‑66 (MQ=255) gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa > 1:1134028/1‑66 (MQ=255) gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGc < 1:3392997/71‑1 (MQ=255) tttCTCGCCAATTCCGGCAATGTCGATGGTTTCTTa > 1:118393/1‑36 (MQ=255) ttCTCGCCGATCCCGGCGATGTCGATGGTTTCTTa > 1:2212534/1‑35 (MQ=25) tctcGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGATATCTCTCGCTGATTGGCgg > 1:503240/1‑70 (MQ=255) ctcGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTa > 1:34490/1‑71 (MQ=255) gCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCAcct < 1:2778830/71‑3 (MQ=255) gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTg < 1:609973/41‑1 (MQ=255) gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTa > 1:3001070/1‑48 (MQ=255) gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa > 1:754515/1‑59 (MQ=255) gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTAc > 1:3864520/1‑69 (MQ=255) gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTAc > 1:1179512/1‑69 (MQ=255) gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt < 1:604534/71‑1 (MQ=255) gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt < 1:3831561/71‑1 (MQ=255) gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt < 1:741763/71‑1 (MQ=255) gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt > 1:2991987/1‑71 (MQ=255) ccAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTg > 1:2483212/1‑71 (MQ=255) cAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGc > 1:3388541/1‑71 (MQ=255) cAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGc < 1:3323066/71‑1 (MQ=255) aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTgg > 1:3641737/1‑60 (MQ=255) aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCt > 1:1457850/1‑71 (MQ=255) | TGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCT > minE/832048‑832185 |
| 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: ‑ |