| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 836,047 | G→A | 14.5% | Q255Q (CAG→CAA) | 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 | 836,047 | 0 | G | A | 14.5% | 137.1 / 13.5 | 55 | Q255Q (CAG→CAA) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base G (18/29); new base A (2/6); total (20/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.96e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01 | |||||||||||
CTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGAC > minE/835981‑836116 | cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGcc < 1:3729428/69‑1 (MQ=255)cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGcc < 1:2234600/69‑1 (MQ=255)cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCg < 1:1050029/70‑1 (MQ=255)cTTCAACTCCAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAACCGGGTCAGcc < 1:224040/69‑1 (MQ=255) cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg > 1:1174327/1‑71 (MQ=255) cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg > 1:39260/1‑71 (MQ=255) cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg > 1:2590286/1‑71 (MQ=255) aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt < 1:514625/71‑1 (MQ=255) aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt < 1:1229407/71‑1 (MQ=255) tGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGtgtg > 1:2997512/1‑71 (MQ=255) gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCt < 1:534380/71‑1 (MQ=255) gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAg < 1:431094/71‑1 (MQ=255) agCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAg > 1:2151562/1‑70 (MQ=255) agCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGa < 1:78065/71‑1 (MQ=255) gCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGa < 1:2602887/70‑1 (MQ=255) cGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgtgt < 1:2625560/62‑1 (MQ=255) caaGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTc < 1:1988407/62‑1 (MQ=255) aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCt < 1:730330/70‑1 (MQ=255) aaGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtg < 1:1467286/71‑1 (MQ=255) aGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgt > 1:1183589/1‑71 (MQ=255) aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc < 1:1451602/71‑2 (MQ=255) aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc < 1:1001021/71‑2 (MQ=255) aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc < 1:57487/71‑2 (MQ=255) gTCTGAGAAGTTCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgtg < 1:3809722/71‑1 (MQ=255) tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgtg < 1:3616713/70‑1 (MQ=255) cTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtc > 1:1855850/1‑71 (MQ=255) cTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtc > 1:2418464/1‑71 (MQ=255) tGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcg > 1:713415/1‑71 (MQ=255) agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtctgt < 1:3299085/71‑1 (MQ=255) aaGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGAcc < 1:3748508/45‑1 (MQ=255) cATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCt < 1:2894151/47‑1 (MQ=255) cATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCt < 1:2889402/47‑1 (MQ=255) aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATc < 1:2031284/69‑1 (MQ=255) cttctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc < 1:2570468/43‑1 (MQ=255) ttctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcggt > 1:386489/1‑60 (MQ=255) ttctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa > 1:552364/1‑71 (MQ=255) ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGtt < 1:3322880/68‑1 (MQ=255) ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa < 1:1342928/69‑1 (MQ=255) tGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt > 1:3160639/1‑71 (MQ=255) cTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATc > 1:1987431/1‑68 (MQ=255) cTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt < 1:2158343/69‑1 (MQ=255) cTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTg > 1:1706436/1‑71 (MQ=255) gTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtg < 1:949474/71‑1 (MQ=255) ttGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg > 1:3530571/1‑71 (MQ=255) ttGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg > 1:2328352/1‑71 (MQ=255) tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGc < 1:3339722/65‑1 (MQ=255) tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGc < 1:525683/65‑1 (MQ=255) gAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttgt < 1:3000923/70‑1 (MQ=255) gggTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgt < 1:3573289/62‑1 (MQ=255) cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGacg > 1:655405/2‑69 (MQ=255) ggTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGc < 1:3796849/71‑1 (MQ=255) ggTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGc > 1:1508657/1‑71 (MQ=255) gTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGcc < 1:3322266/71‑1 (MQ=255) aGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGa > 1:2824103/1‑70 (MQ=255) aGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGAc > 1:3709775/1‑71 (MQ=255) | CTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGAC > minE/835981‑836116 |
| 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: ‑ |