| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 833,867 | T→C | 13.3% | R775R (CGT→CGC) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 833,867 | 0 | T | C | 13.3% | 179.7 / 17.4 | 75 | R775R (CGT→CGC) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base T (38/27); new base C (5/5); total (43/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.36e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01 | |||||||||||
AGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACA > minE/833801‑833934 | agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTg > 1:5303677/1‑70 (MQ=255)agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt > 1:4526262/1‑71 (MQ=255)agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt > 1:4808825/1‑71 (MQ=255)agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt > 1:4832198/1‑71 (MQ=255)agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt > 1:2941702/1‑71 (MQ=255)agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt < 1:5625951/71‑1 (MQ=255) tGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt < 1:1380689/66‑1 (MQ=255) tGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt < 1:3412403/66‑1 (MQ=255) aCAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGt > 1:3478109/1‑50 (MQ=255) aCAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGt > 1:5556828/1‑50 (MQ=255) aCAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGt > 1:1780554/1‑50 (MQ=255) aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt < 1:443498/52‑1 (MQ=255) aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTc > 1:1998311/1‑71 (MQ=255) aTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTcc < 1:5286047/71‑1 (MQ=255) ggTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGc > 1:750484/1‑54 (MQ=255) gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAc < 1:4024818/71‑1 (MQ=255) attgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGatat < 1:2025996/68‑1 (MQ=255) tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg < 1:663837/49‑1 (MQ=255) tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg < 1:3298717/49‑1 (MQ=255) tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg < 1:2182199/66‑1 (MQ=255) ggAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGc < 1:4466279/49‑1 (MQ=255) ggAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGc < 1:311481/49‑1 (MQ=255) ggAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca > 1:5251164/1‑71 (MQ=255) gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa > 1:5042738/1‑65 (MQ=255) gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatcat > 1:5199391/1‑71 (MQ=255) gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca < 1:1612153/70‑1 (MQ=255) gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTcca > 1:557420/1‑44 (MQ=255) gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGAtatcg < 1:4205290/70‑2 (MQ=255) aaGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtt > 1:2428842/1‑71 (MQ=255) gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTg < 1:5113514/71‑1 (MQ=255) cAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCAc > 1:796553/1‑45 (MQ=255) aGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg < 1:2521653/40‑1 (MQ=255) aGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACGTGTCTCTATTCCGACATCAttt > 1:4358453/1‑66 (MQ=255) aGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGcc > 1:4173647/1‑70 (MQ=255) gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca > 1:110470/1‑62 (MQ=255) gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGa < 1:4925303/71‑1 (MQ=255) gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGa > 1:5493622/1‑71 (MQ=255) cTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGa > 1:3182038/1‑37 (MQ=255) cTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGa > 1:4072121/1‑37 (MQ=255) ggATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTACGAcatca > 1:5185177/1‑59 (MQ=255) gATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGc < 1:2478558/71‑1 (MQ=255) gATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGc < 1:5279210/71‑1 (MQ=255) aTCTGGAGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGc > 1:728416/1‑70 (MQ=255) cTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa > 1:1669180/1‑50 (MQ=255) cTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa > 1:785298/1‑50 (MQ=255) gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGActt < 1:3872112/70‑3 (MQ=255) gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtt < 1:983705/51‑1 (MQ=255) gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtt < 1:3865070/51‑1 (MQ=255) gtACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGcc < 1:3076230/54‑1 (MQ=255) aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg < 1:1473821/70‑1 (MQ=255) aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg < 1:2530617/70‑1 (MQ=255) aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg > 1:2414707/1‑70 (MQ=255) aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg < 1:4562829/70‑1 (MQ=255) cGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa < 1:774887/71‑1 (MQ=255) gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa < 1:5376024/71‑1 (MQ=255) gCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa < 1:3798898/71‑1 (MQ=255) tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaaa > 1:1865983/3‑70 (MQ=255) cgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa > 1:2054885/2‑67 (MQ=255) gACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaa < 1:227333/69‑1 (MQ=255) gACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaa < 1:2869982/69‑1 (MQ=255) gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaaaa < 1:2516655/69‑1 (MQ=255) gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga > 1:3497426/1‑71 (MQ=255) gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga > 1:3490634/1‑71 (MQ=255) aCTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGcc > 1:3902865/1‑46 (MQ=255) ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAgacga > 1:2829301/1‑71 (MQ=255) ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc > 1:4541119/1‑51 (MQ=255) ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa > 1:5451873/1‑71 (MQ=255) ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa > 1:252478/1‑71 (MQ=255) ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa > 1:3290762/1‑71 (MQ=255) ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa > 1:3333857/1‑71 (MQ=255) ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa > 1:3799207/1‑71 (MQ=255) ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa > 1:1729023/1‑71 (MQ=255) ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGCCGACCTGGTACGAAAAAGACGACa > 1:3451766/1‑71 (MQ=255) ccGTCTCTCGAGCACCTGTCTCTATTCCGAcatca > 1:4583525/1‑35 (MQ=255) cGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa > 1:4675448/1‑61 (MQ=255) cGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa > 1:590478/1‑70 (MQ=255) gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaag > 1:5404182/1‑59 (MQ=255) gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaa > 1:1978027/1‑61 (MQ=255) | AGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACA > minE/833801‑833934 |
| 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: ‑ |