| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 837,193 | C→T | 27.8% | Q126* (CAA→TAA) | narJ → | molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 837,193 | 0 | C | T | 27.8% | 52.3 / 21.0 | 36 | Q126* (CAA→TAA) | narJ | molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1 |
| Reads supporting (aligned to +/- strand): ref base C (9/17); new base T (2/8); total (11/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.88e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01 | |||||||||||
GCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCT > minE/837126‑837257 | gCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGTaa < 1:3363114/70‑1 (MQ=255)gCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCaa < 1:1246978/70‑1 (MQ=255) cAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAg < 1:3761396/71‑1 (MQ=255) aaCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGc < 1:1688783/71‑1 (MQ=255) aCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGcc > 1:3468125/1‑71 (MQ=255) aGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAg < 1:1324160/62‑1 (MQ=255) gcgAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGa < 1:2209957/70‑1 (MQ=255) aaTTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCa < 1:1179693/52‑1 (MQ=255) aaTTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGt > 1:2792818/1‑71 (MQ=255) aaTTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGt < 1:1941141/71‑1 (MQ=255) aTTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGtt < 1:2864906/71‑1 (MQ=255) ttGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGtt > 1:403507/1‑70 (MQ=255) gCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTg > 1:1679842/1‑70 (MQ=255) gACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCg < 1:1575939/49‑1 (MQ=255) tCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGttt < 1:167588/60‑1 (MQ=255) tATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTc < 1:892506/70‑1 (MQ=255) tATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTc < 1:28987/70‑1 (MQ=255) aTCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTg > 1:2219118/1‑71 (MQ=255) tCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGcc < 1:3331738/37‑1 (MQ=255) tCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGcc < 1:1258813/37‑1 (MQ=255) tCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATcgc > 1:2358735/1‑60 (MQ=255) cTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTgg > 1:1380675/1‑70 (MQ=255) tGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCa > 1:2129851/1‑71 (MQ=255) gAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTgg < 1:280494/67‑1 (MQ=255) gAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCAt < 1:2182645/70‑1 (MQ=255) tACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGata < 1:3450850/49‑1 (MQ=255) tACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGata < 1:1041071/49‑1 (MQ=255) tACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGata < 1:1036526/49‑1 (MQ=255) cTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGatat < 1:1457310/47‑1 (MQ=255) cTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGatat < 1:903128/47‑1 (MQ=255) ggCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGa > 1:1926539/1‑38 (MQ=255) aGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCt < 1:988042/52‑1 (MQ=255) aGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGt > 1:1306452/1‑71 (MQ=255) aGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGt < 1:579914/71‑1 (MQ=255) cTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCt > 1:1112303/1‑71 (MQ=255) cTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCt < 1:3410863/71‑1 (MQ=255) | GCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCT > minE/837126‑837257 |
| 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: ‑ |