| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 833,993 | C→G | 33.3% | A817A (GCC→GCG) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 833,993 | 0 | C | G | 33.3% | 48.1 / 40.7 | 48 | A817A (GCC→GCG) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base C (13/19); new base G (4/12); total (17/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.50e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01 | |||||||||||
AAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGT > minE/833924‑834058 | aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt < 1:2567963/71‑1 (MQ=255)aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt > 1:3946166/1‑71 (MQ=255) aaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTg < 1:3714343/71‑1 (MQ=255) aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTg > 1:2116689/1‑71 (MQ=255) aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTg < 1:942542/71‑1 (MQ=255) cgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGaa < 1:3126815/70‑1 (MQ=255) cgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg < 1:404217/71‑1 (MQ=255) ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGaa < 1:2987246/66‑1 (MQ=255) aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTgg < 1:1288601/59‑1 (MQ=255) aCTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGa < 1:3113469/59‑1 (MQ=255) tCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTg > 1:2156777/1‑71 (MQ=255) cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:3257707/1‑71 (MQ=255) cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:681132/1‑71 (MQ=255) ggATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTgg < 1:1921612/70‑1 (MQ=255) gATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGa < 1:1009121/71‑1 (MQ=255) atatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaa < 1:110971/71‑1 (MQ=255) tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGaaa < 1:2340387/71‑1 (MQ=255) tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaaa < 1:2129344/71‑1 (MQ=255) cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa < 1:2271006/70‑1 (MQ=255) cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa < 1:429712/70‑1 (MQ=255) ttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc > 1:583890/1‑70 (MQ=255) tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggta > 1:1205943/1‑68 (MQ=255) tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggt < 1:1102742/70‑3 (MQ=255) tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGaa > 1:2600892/1‑40 (MQ=255) ttCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATc > 1:1281837/1‑71 (MQ=255) tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggt < 1:1713221/63‑3 (MQ=255) tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca < 1:2048209/70‑3 (MQ=255) tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca < 1:781429/70‑3 (MQ=255) tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca < 1:292052/70‑3 (MQ=255) tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca < 1:921612/70‑3 (MQ=255) cGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca > 1:3470978/1‑67 (MQ=255) cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg > 1:2037592/1‑68 (MQ=255) cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg > 1:1838794/1‑68 (MQ=255) gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag > 1:1536389/1‑70 (MQ=255) gTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATc < 1:2554677/49‑1 (MQ=255) tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg < 1:2575807/63‑1 (MQ=255) tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg < 1:2063502/63‑1 (MQ=255) cggTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGaa < 1:1211621/71‑1 (MQ=255) gTCGATACGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAg > 1:210192/1‑70 (MQ=255) tCGATCCGGCCTGGGAAGCGAAACGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAAt < 1:1941501/60‑1 (MQ=255) tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaaa < 1:3122900/36‑1 (MQ=255) tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCg < 1:1616871/66‑1 (MQ=255) tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtg > 1:2782389/1‑71 (MQ=255) cGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg < 1:2342961/52‑1 (MQ=255) aTCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCcaaa < 1:3020136/53‑2 (MQ=255) tCCGGCCTTGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg < 1:945227/71‑1 (MQ=255) ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt > 1:3054163/1‑71 (MQ=255) ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt > 1:1433009/1‑71 (MQ=255) cTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATtct > 1:3820388/1‑53 (MQ=255) | AAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGT > minE/833924‑834058 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |