| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 833,966 | C→T | 21.6% | H808H (CAC→CAT) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 833,966 | 0 | C | T | 21.6% | 105.8 / 24.2 | 51 | H808H (CAC→CAT) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base C (17/23); new base T (9/2); total (26/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.85e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01 | |||||||||||
CATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCAT > minE/833897‑834031 | catTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAcc > 1:939747/1‑71 (MQ=255)catTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAcc > 1:2515037/1‑71 (MQ=255)catTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAcc > 1:491623/1‑71 (MQ=255) atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAc > 1:3445131/1‑69 (MQ=255) atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAc > 1:3112857/1‑69 (MQ=255) ttttGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCg > 1:3592799/1‑71 (MQ=255) aCCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGt < 1:1994593/71‑1 (MQ=255) ccTGGTACGAAAAACACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTc < 1:336103/71‑1 (MQ=255) tGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa > 1:2608466/1‑71 (MQ=255) ggTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGAt > 1:1187164/1‑71 (MQ=255) ggTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGAt > 1:3901467/1‑71 (MQ=255) gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGAt < 1:2583304/70‑1 (MQ=255) gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATc > 1:1557574/1‑71 (MQ=255) gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa < 1:1155464/69‑1 (MQ=255) gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa < 1:170601/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:3654128/1‑71 (MQ=255) aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTg > 1:1927852/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:2682221/1‑70 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:2320193/70‑1 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:2382236/70‑1 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCg < 1:4019447/71‑1 (MQ=255) tACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCg < 1:3559300/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) | CATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCAT > minE/833897‑834031 |
| 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: ‑ |