| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 834,296 | C→G | 20.9% | A918A (GCC→GCG) ‡ | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 834,296 | 0 | C | G | 20.9% | 139.9 / 30.0 | 67 | A918A (GCC→GCG) ‡ | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base C (22/31); new base G (10/4); total (32/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.05e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01 | |||||||||||
ACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGA > minE/834228‑834363 | aCTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGG‑AAAGGGATTGCCTgg < 1:419840/71‑1 (MQ=255) aCGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGaa < 1:291860/70‑1 (MQ=255) aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg < 1:267508/65‑1 (MQ=255) aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc > 1:4518783/1‑71 (MQ=255) gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg < 1:3178416/62‑1 (MQ=255) cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGCTAATGGCGGTAAAGGGATTGCCTGGAACACCCaga < 1:761213/71‑1 (MQ=255) tACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc < 1:3115838/70‑1 (MQ=255) tACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATACCCaggat < 1:4847850/70‑4 (MQ=255) aCCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCg < 1:5329022/70‑1 (MQ=255) ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga < 1:2660199/70‑1 (MQ=255) ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga < 1:3160474/70‑1 (MQ=255) ctctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga < 1:4396547/71‑1 (MQ=255) ctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc < 1:210801/65‑1 (MQ=255) tATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATg > 1:1766222/1‑70 (MQ=255) aTCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGa > 1:894272/1‑71 (MQ=255) aTCGGCCCGCTGATGGAGAAAATAGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGa > 1:4372801/1‑71 (MQ=255) tCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg < 1:4481599/50‑1 (MQ=255) tCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGAt > 1:3797728/1‑71 (MQ=255) cGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATgg < 1:802142/68‑1 (MQ=255) cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg < 1:2551366/44‑1 (MQ=255) cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg < 1:4663557/44‑1 (MQ=255) cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg < 1:2845189/44‑1 (MQ=255) cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga < 1:5435845/58‑1 (MQ=255) cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga < 1:3581008/58‑1 (MQ=255) gCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg < 1:1579095/43‑1 (MQ=255) cTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgt < 1:4322964/42‑2 (MQ=255) gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa > 1:5430540/1‑42 (MQ=16) gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa > 1:2891107/1‑42 (MQ=16) gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa > 1:5658164/1‑42 (MQ=16) gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa > 1:935884/1‑42 (MQ=16) gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGa > 1:198260/1‑41 (MQ=16) gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGa > 1:1536115/1‑41 (MQ=16) gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc > 1:515111/1‑46 (MQ=255) gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc > 1:2966700/1‑46 (MQ=255) gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc > 1:2030043/1‑46 (MQ=255) tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc < 1:4887448/45‑1 (MQ=255) tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc < 1:2677410/45‑1 (MQ=255) tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc < 1:1775653/45‑1 (MQ=255) agaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCt > 1:998717/1‑70 (MQ=255) agaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCt > 1:5344887/1‑70 (MQ=255) gaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCa < 1:4824569/71‑1 (MQ=255) aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc < 1:205472/71‑1 (MQ=255) aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc > 1:1411784/1‑71 (MQ=255) aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc < 1:792991/71‑1 (MQ=255) aaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga > 1:3745241/1‑47 (MQ=255) aaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga > 1:2346196/1‑47 (MQ=255) aaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga < 1:3094082/45‑1 (MQ=255) aTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTa < 1:5024182/71‑1 (MQ=255) aTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCAAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTa > 1:665690/1‑71 (MQ=255) ggTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTacac > 1:3989957/1‑71 (MQ=255) aTGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga > 1:4239199/1‑39 (MQ=255) aTGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCa > 1:3679147/1‑69 (MQ=255) aTGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCa > 1:2032368/1‑69 (MQ=255) ggCGGTAAAGGGATTTCGTGGAATACCCAGGATGAAATCGATTTCCTCGGCAAACTCAATTACACCaagc < 1:1823331/70‑3 (MQ=255) ggCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTaca > 1:1385535/1‑64 (MQ=255) gCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGc < 1:3435614/70‑1 (MQ=255) aaaGGGATTGCCTGGAACACCCAGAGCGAGAGGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAgg > 1:4798302/1‑71 (MQ=255) aaGGGATTGCCTGGAACACCCAGAGCGAGATGGAt > 1:1149982/1‑35 (MQ=255) gggATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTcc < 1:3989496/71‑1 (MQ=255) gggATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTcc > 1:761562/1‑71 (MQ=255) aTTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGaa > 1:2290267/1‑63 (MQ=255) tGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAg > 1:2595004/1‑62 (MQ=255) tGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAg > 1:1857751/1‑62 (MQ=255) tGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCg > 1:1114380/1‑70 (MQ=255) gCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGa < 1:5102534/70‑1 (MQ=255) gCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGa < 1:1887641/70‑1 (MQ=255) gCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGa < 1:827078/70‑1 (MQ=255) | ACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGA > minE/834228‑834363 |
| 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: ‑ |