| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,858 | T→C | 12.7% | C192C (TGT→TGC) | narH → | nitrate reductase 1, beta (Fe‑S) subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 835,858 | 0 | T | C | 12.7% | 165.9 / 15.4 | 63 | C192C (TGT→TGC) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base T (31/24); new base C (5/3); total (36/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01 | |||||||||||
TAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATA > minE/835789‑835928 | taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg > 1:1669244/1‑71 (MQ=255) gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtgt > 1:2683894/1‑70 (MQ=255) cAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGc > 1:761217/1‑71 (MQ=255) ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg < 1:1640584/64‑1 (MQ=255) ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc < 1:3718816/71‑1 (MQ=255) ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc < 1:904302/71‑1 (MQ=255) tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt > 1:1518126/2‑71 (MQ=255) cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt > 1:2410564/1‑71 (MQ=255) gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt < 1:751524/70‑1 (MQ=255) aaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGT‑GCGACCTGcc > 1:2527627/1‑69 (MQ=255) acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCg < 1:3555968/71‑1 (MQ=255) ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGc < 1:3346729/71‑1 (MQ=255) cTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGc < 1:117428/70‑1 (MQ=255) cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc > 1:47455/1‑57 (MQ=255) cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCgg < 1:1394333/69‑1 (MQ=255) cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCgg < 1:3516571/69‑1 (MQ=255) cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg > 1:3527305/1‑69 (MQ=255) cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg > 1:117734/1‑69 (MQ=255) cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg > 1:3600629/1‑69 (MQ=255) atgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGAGGTGc > 1:3149782/1‑71 (MQ=255) tATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGccc > 1:1951244/1‑56 (MQ=255) tATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAttt < 1:3044142/70‑1 (MQ=255) aTTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg > 1:3106394/1‑42 (MQ=255) tttGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg > 1:1395579/1‑41 (MQ=255) tttGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAtt < 1:1673720/67‑1 (MQ=255) tttGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc > 1:1111159/1‑70 (MQ=255) tttGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc > 1:3463761/1‑70 (MQ=255) ctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTActa < 1:1585897/70‑3 (MQ=255) tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGa > 1:1564231/1‑45 (MQ=255) tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg > 1:1829469/1‑71 (MQ=255) tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg > 1:4045314/1‑71 (MQ=255) gcgcCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt < 1:624428/60‑1 (MQ=255) gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg < 1:2385307/71‑3 (MQ=255) gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg < 1:373081/71‑3 (MQ=255) cgcCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg > 1:1196786/1‑70 (MQ=255) ccTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag < 1:2951842/71‑1 (MQ=255) tgtgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag > 1:3680443/1‑69 (MQ=255) tgtgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag > 1:1926426/1‑69 (MQ=255) tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:2590296/69‑1 (MQ=255) tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:669598/69‑1 (MQ=255) tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag < 1:2782525/70‑1 (MQ=255) cGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg > 1:1010523/1‑71 (MQ=255) cGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg > 1:3143601/1‑71 (MQ=255) aaCACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg > 1:1651483/1‑70 (MQ=255) aaCACTGCCTCAA‑TCCCAGATGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACgg > 1:2350301/1‑70 (MQ=255) acacTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa > 1:2775066/1‑64 (MQ=255) acacTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaa > 1:3369353/1‑61 (MQ=255) acacTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCa < 1:71421/71‑1 (MQ=255) cacTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGa < 1:3839322/46‑1 (MQ=255) acTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag < 1:3717996/60‑1 (MQ=255) tGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:294023/60‑1 (MQ=255) tGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGt > 1:3699780/1‑71 (MQ=255) tGCCTGAA‑‑CCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:39978/59‑1 (MQ=255) ccTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc > 1:4021566/1‑71 (MQ=255) ccTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc > 1:911336/1‑71 (MQ=255) cTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt < 1:2695292/71‑1 (MQ=255) cTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt < 1:1308139/71‑1 (MQ=255) aa‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt > 1:757825/1‑68 (MQ=255) a‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga > 1:1386404/1‑55 (MQ=255) a‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc > 1:3624603/1‑71 (MQ=255) a‑CCCGGCATGTGTGGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga > 1:3363334/1‑55 (MQ=255) cccGGCATGTGTGGCGACCTGCGCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc > 1:4018316/1‑70 (MQ=255) cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc > 1:3070801/1‑70 (MQ=255) ccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:3014250/51‑1 (MQ=255) cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc < 1:4018509/71‑1 (MQ=255) ggCATGTGTGGCGACCTGCCCGAGCGGTGCGAttt > 1:2682848/1‑35 (MQ=255) ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga > 1:3533152/1‑51 (MQ=255) ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga > 1:1090911/1‑51 (MQ=255) tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:3252888/45‑1 (MQ=255) tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:1350903/45‑1 (MQ=255) tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGc < 1:1662162/51‑1 (MQ=255) tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg < 1:2435676/60‑1 (MQ=255) tgtgtgGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGa < 1:1152876/71‑1 (MQ=255) gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt < 1:1035948/58‑1 (MQ=255) gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt < 1:3808084/58‑1 (MQ=255) gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt < 1:35998/58‑1 (MQ=255) tgtgGCGACCTGCCCGAGCGGTTCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa < 1:3843982/71‑1 (MQ=255) tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa > 1:1343375/1‑43 (MQ=255) tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt < 1:3039922/70‑1 (MQ=255) | TAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATA > minE/835789‑835928 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |