| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 833,972 | G→T | 23.2% | L810L (CTG→CTT) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 833,972 | 0 | G | T | 23.2% | 144.2 / 24.1 | 56 | L810L (CTG→CTT) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (16/27); new base T (10/3); total (26/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.40e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01 | |||||||||||
ACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAAT > minE/833913‑834042 | aCCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGt < 1:1994593/71‑1 (MQ=255) ccTGGTACGAAAAACACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTc < 1:336103/71‑1 (MQ=255) tGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa > 1:2608466/1‑71 (MQ=255) ggTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGAt > 1:3901467/1‑71 (MQ=255) ggTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGAt > 1:1187164/1‑71 (MQ=255) gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGAt < 1:2583304/70‑1 (MQ=255) gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATc > 1:1557574/1‑71 (MQ=255) gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa < 1:170601/69‑1 (MQ=255) gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa < 1:1155464/69‑1 (MQ=255) tACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATcc > 1:3615375/1‑71 (MQ=255) aaaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGcgt > 1:419502/1‑69 (MQ=255) aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa < 1:2617856/63‑1 (MQ=255) aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTg > 1:1927852/1‑71 (MQ=255) aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTg > 1:3654128/1‑71 (MQ=255) aaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTgg < 1:3488708/71‑1 (MQ=255) aGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTggg < 1:3506043/71‑1 (MQ=255) acATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCg < 1:1667585/71‑1 (MQ=255) gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGaa > 1:831315/1‑64 (MQ=255) gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGaa > 1:3700578/1‑64 (MQ=255) aaTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGaaaa < 1:1804756/70‑1 (MQ=255) aaTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAg > 1:2528268/1‑71 (MQ=255) aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCg < 1:3236200/48‑1 (MQ=255) aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc > 1:3514075/1‑71 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCg < 1:2408507/52‑1 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTggg < 1:1392849/59‑1 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:2382236/70‑1 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:2320193/70‑1 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc > 1:2682221/1‑70 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCg < 1:3559300/71‑1 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCg < 1:4019447/71‑1 (MQ=255) cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGAc < 1:3628093/71‑1 (MQ=255) ttCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACt > 1:2570644/1‑71 (MQ=255) tCGGATATGCATCCGTTTATTCACACGCTGTCTGcggcgg > 1:1644662/1‑40 (MQ=255) cGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTgg > 1:712848/1‑71 (MQ=255) ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTgg > 1:2046069/1‑52 (MQ=38) ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTgg > 1:2041163/1‑52 (MQ=38) atatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaa < 1:3058670/71‑1 (MQ=255) tGCATCCGTTTATTCACCCGCTTTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATc > 1:848939/1‑71 (MQ=255) tGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:3124740/59‑1 (MQ=255) tGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:647183/59‑1 (MQ=255) gCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATcc > 1:1868350/1‑39 (MQ=37) gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATc < 1:111319/70‑1 (MQ=255) cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAAt > 1:909006/1‑64 (MQ=255) cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGa > 1:3973903/1‑55 (MQ=255) tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGaa < 1:2501540/47‑1 (MQ=255) tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc < 1:1345412/71‑1 (MQ=255) ttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc > 1:271761/1‑70 (MQ=255) tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGcc < 1:1142873/70‑1 (MQ=255) aTTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCAt > 1:1229310/1‑71 (MQ=255) aCCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAAt > 1:33336/1‑55 (MQ=255) ccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa < 1:2234632/59‑1 (MQ=255) gCTGTCTGCGGCGGTCGATCCTGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaaga < 1:3680703/71‑1 (MQ=255) gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaaga < 1:1948288/71‑1 (MQ=255) gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaaga < 1:1511280/71‑1 (MQ=255) cTGGCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaaga < 1:128519/70‑1 (MQ=255) tGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGaaa > 1:3088100/1‑71 (MQ=255) gTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAAt > 1:2487114/1‑71 (MQ=255) | ACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAAT > minE/833913‑834042 |
| 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: ‑ |