| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 834,221 | T→C | 22.9% | Y893Y (TAT→TAC) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 834,221 | 0 | T | C | 22.9% | 109.8 / 32.3 | 61 | Y893Y (TAT→TAC) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base T (34/13); new base C (9/5); total (43/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.39e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 | |||||||||||
GAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAA > minE/834152‑834286 | gAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTAt < 1:1140874/70‑1 (MQ=255) ggCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGAc > 1:296099/1‑71 (MQ=255) ggCGAATGCGATCTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTAc < 1:834024/71‑1 (MQ=255) gCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACt < 1:3167794/71‑1 (MQ=255) cGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACtt < 1:3174812/71‑1 (MQ=255) gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc > 1:92043/1‑70 (MQ=255) tGCGATCTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGa < 1:3664380/71‑1 (MQ=255) gACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc < 1:2486990/66‑1 (MQ=255) ccTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCt > 1:953478/1‑71 (MQ=255) cTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCtt < 1:2009306/71‑1 (MQ=255) tGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCtt > 1:1517072/1‑70 (MQ=255) gATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTa > 1:2449633/1‑71 (MQ=255) attccAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCtt < 1:2850084/65‑1 (MQ=255) attccAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTAc > 1:1821970/4‑71 (MQ=255) aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc < 1:600680/71‑1 (MQ=255) ttccAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTAcc > 1:1029086/3‑71 (MQ=255) ttccAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTAcc < 1:708956/69‑1 (MQ=255) ttccAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTAcc > 1:3436047/3‑71 (MQ=255) tccAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACct > 1:198545/2‑71 (MQ=255) ccAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCgg > 1:1851382/1‑49 (MQ=255) ccAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAAc < 1:223077/61‑1 (MQ=255) cAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctc > 1:1208483/1‑70 (MQ=255) cAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACctcg > 1:1074586/1‑70 (MQ=255) cAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACctcg > 1:2751193/1‑70 (MQ=255) aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGc < 1:1726855/48‑1 (MQ=255) aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa > 1:388248/1‑58 (MQ=255) aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc > 1:3209724/1‑66 (MQ=255) aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa > 1:3951602/1‑71 (MQ=255) aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa > 1:378712/1‑71 (MQ=255) ggTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTAt > 1:3084154/1‑71 (MQ=255) gTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACt > 1:799961/1‑50 (MQ=255) aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg < 1:19573/71‑1 (MQ=255) aaaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCgg > 1:2444626/1‑71 (MQ=255) aaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGc > 1:1980575/1‑71 (MQ=255) aaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGc > 1:186387/1‑71 (MQ=255) aaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgggc > 1:1391442/1‑69 (MQ=255) aCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACt < 1:3739987/45‑1 (MQ=255) aCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGccc < 1:2561972/71‑1 (MQ=255) ccGCGCCACACATTATGGTCGTAGAGCGCGATTATca > 1:1461547/1‑36 (MQ=255) cgcgCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCg > 1:1932590/1‑70 (MQ=255) gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt > 1:1281191/1‑71 (MQ=255) cacacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGAt < 1:3158329/70‑1 (MQ=255) cacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg > 1:2937784/1‑69 (MQ=255) tATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTa > 1:2014898/1‑44 (MQ=255) aTGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAAt > 1:320999/1‑71 (MQ=255) ggTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct < 1:3806373/46‑1 (MQ=255) agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg > 1:603530/1‑71 (MQ=255) agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg > 1:3275776/1‑71 (MQ=255) gagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATgg < 1:3911881/71‑1 (MQ=255) agCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGc > 1:1367208/1‑71 (MQ=255) cgcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCgg > 1:2996914/1‑71 (MQ=255) cgcgACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATg > 1:3911127/1‑51 (MQ=255) cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt > 1:558700/1‑70 (MQ=255) gATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTaa < 1:1628483/71‑1 (MQ=255) gATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTa > 1:1968555/1‑70 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGa > 1:2904507/1‑44 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGa > 1:2842278/1‑44 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg > 1:132678/1‑46 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg > 1:102027/1‑46 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa > 1:306251/1‑60 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa > 1:2718844/1‑60 (MQ=255) aTCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGa > 1:75292/1‑42 (MQ=255) | GAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAA > minE/834152‑834286 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |