| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 834,515 | A→G | 47.8% | E991E (GAA→GAG) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 834,515 | 0 | A | G | 47.8% | 2.1 / 58.0 | 46 | E991E (GAA→GAG) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base A (15/9); new base G (14/8); total (29/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCT > minE/834450‑834584 | aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAagaaga > 1:5674681/1‑68 (MQ=255)aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGAc > 1:2851789/1‑69 (MQ=255)aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGAc > 1:1786923/1‑69 (MQ=255) aaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACgag < 1:5475302/71‑1 (MQ=255) cTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCAGGCGCTGAATAAAGAAGACGAGAAGATCCGc > 1:4356105/1‑71 (MQ=255) cTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCcg > 1:2001859/1‑71 (MQ=255) gAGCGAATTTACCGGTCGTGACCATACGCATCTGGGGCTGAATAAAGAAGACGAGAAGATCCGCTTCcgc > 1:3195853/1‑70 (MQ=255) aGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa < 1:2626879/71‑1 (MQ=255) aGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa < 1:1801173/71‑1 (MQ=255) aGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa > 1:1934889/1‑71 (MQ=255) aaTTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGatat > 1:3871637/1‑70 (MQ=255) aaTTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATAtt < 1:3332867/71‑1 (MQ=255) aTTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATAtt > 1:3796856/1‑70 (MQ=255) taacGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCa > 1:4472993/4‑69 (MQ=255) tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACgaga > 1:2373374/1‑47 (MQ=255) tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCa > 1:1068834/1‑69 (MQ=255) cGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGata < 1:266052/62‑1 (MQ=255) cGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg > 1:5275372/1‑69 (MQ=255) cGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg > 1:571306/1‑69 (MQ=255) ggTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcaca < 1:1678312/71‑1 (MQ=255) cGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGata > 1:2690224/1‑58 (MQ=255) aacATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgcg > 1:5178234/3‑70 (MQ=255) acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGc > 1:650790/2‑65 (MQ=255) acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGc > 1:2943682/2‑65 (MQ=255) acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgcgt > 1:4818750/2‑69 (MQ=255) acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgcgt > 1:1830055/2‑69 (MQ=255) acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgcgt > 1:363581/2‑69 (MQ=255) cATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCaa < 1:5646884/71‑1 (MQ=255) cATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTaa < 1:4841285/71‑1 (MQ=255) cATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTaa < 1:4857332/71‑1 (MQ=255) aTACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCg < 1:2931192/40‑1 (MQ=255) tACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTaaaa > 1:4450702/1‑71 (MQ=255) cccATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAAtt < 1:1796395/69‑1 (MQ=255) cGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCcgc > 1:1926493/1‑63 (MQ=255) ccATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTa < 1:1309694/70‑1 (MQ=255) tCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAg < 1:4613337/49‑1 (MQ=255) tGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGatat > 1:2538345/1‑43 (MQ=19) tGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGat > 1:4701327/1‑41 (MQ=19) tGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGat > 1:576717/1‑41 (MQ=19) cgcTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAAtt < 1:4708290/62‑1 (MQ=255) gAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGa > 1:2512401/1‑71 (MQ=255) cacaaaGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCAc > 1:2355412/4‑71 (MQ=255) aaTAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCg < 1:1233782/69‑1 (MQ=255) aaTAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGa < 1:4425568/70‑1 (MQ=255) acaaaGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCAcc < 1:5475291/69‑1 (MQ=255) tAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCcgc > 1:632745/1‑47 (MQ=255) aaaGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTAtc > 1:3525364/1‑57 (MQ=255) aagaagACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCaaaa > 1:2592900/1‑51 (MQ=255) aagaagACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCaaaa > 1:3180791/1‑51 (MQ=255) agaggaCGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctggt > 1:1276809/1‑71 (MQ=255) gaagaCGAGAAGATCCGCTTCCGCGATATTCAGGc < 1:4392049/35‑1 (MQ=255) gaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCaaaa < 1:2206360/49‑1 (MQ=255) gaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctg > 1:1353679/1‑68 (MQ=255) gaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACct > 1:5152866/1‑67 (MQ=255) aagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtct > 1:3823411/1‑71 (MQ=255) | AAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCT > minE/834450‑834584 |
| 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: ‑ |