| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 834,584 | T→C | 19.1% | S1014S (TCT→TCC) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 834,584 | 0 | T | C | 19.1% | 96.6 / 16.4 | 47 | S1014S (TCT→TCC) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base T (24/14); new base C (3/6); total (27/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.42e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGG > minE/834516‑834650 | gaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctgg > 1:3154823/1‑71 (MQ=255)gaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctgg > 1:92992/1‑71 (MQ=255) cGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCtggtccggtc > 1:3442243/1‑71 (MQ=255) gagaAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCtggtccggtc < 1:1052559/70‑1 (MQ=255) gaagaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGgtctggtctg > 1:2688955/1‑70 (MQ=255) agaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGa < 1:2507330/70‑1 (MQ=255) gaTTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGa > 1:2016113/1‑69 (MQ=255) gaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGgtct‑gtctgg < 1:1862095/67‑1 (MQ=255) tCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGaa > 1:272111/1‑68 (MQ=255) tCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt > 1:2720303/1‑71 (MQ=255) tCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt > 1:2658065/1‑71 (MQ=255) cGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGa < 1:579043/71‑1 (MQ=255) gCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa < 1:597578/71‑1 (MQ=255) cTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAac < 1:17509/71‑1 (MQ=255) tCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa > 1:3157990/1‑68 (MQ=255) tCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac > 1:162863/1‑71 (MQ=255) ccGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac > 1:1261589/1‑70 (MQ=255) ccGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac < 1:250499/70‑1 (MQ=255) atatTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCt < 1:3872957/71‑1 (MQ=255) tatTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCtt < 1:1375015/71‑1 (MQ=255) ttCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCtggtccggtc < 1:366804/48‑1 (MQ=255) ttCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctg < 1:283827/45‑1 (MQ=255) cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa > 1:3752800/1‑68 (MQ=255) cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa > 1:653439/1‑68 (MQ=255) cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa < 1:1228092/70‑1 (MQ=255) cggCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTAc < 1:550323/67‑1 (MQ=255) aGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAAc > 1:1431713/1‑71 (MQ=255) gcgcaGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACg < 1:3136073/67‑1 (MQ=255) acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGcc < 1:710114/70‑1 (MQ=255) acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg > 1:3223269/1‑71 (MQ=255) acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg > 1:526679/1‑71 (MQ=255) acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg > 1:3656926/1‑71 (MQ=255) gcAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt > 1:715966/1‑41 (MQ=255) gcAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt < 1:3931577/41‑1 (MQ=255) gcAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACAcc < 1:1113689/71‑1 (MQ=255) taaaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCggga < 1:3233394/63‑2 (MQ=255) cAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac > 1:2292212/1‑69 (MQ=255) aaaaTTATCTCCAGCCCGACCTGGTCTGGTCTGGATGATGAACACGTTTCTTACAACGCCGGTTACACCaa > 1:1759730/1‑71 (MQ=255) aaaaTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCtt > 1:3949103/1‑52 (MQ=255) aaaTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAAc > 1:3283398/1‑71 (MQ=255) aaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAAc > 1:3884871/1‑71 (MQ=255) aaTTATCTCCAGCCCGCCCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAAc > 1:2240582/1‑70 (MQ=255) aTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTa > 1:3101921/1‑62 (MQ=255) aTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGt < 1:2394697/71‑1 (MQ=255) aTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGt < 1:2390005/71‑1 (MQ=255) ttATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATa < 1:2956570/61‑1 (MQ=255) ttATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATa < 1:3339452/61‑1 (MQ=255) aTCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAAc > 1:1464103/1‑51 (MQ=255) tctcCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCgg > 1:49419/1‑55 (MQ=255) cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa < 1:3493738/69‑1 (MQ=255) cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc > 1:1885374/1‑71 (MQ=255) aGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCt > 1:898653/1‑39 (MQ=255) ccGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGAt > 1:2931444/1‑71 (MQ=255) cGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATc > 1:2420899/1‑71 (MQ=255) ccTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa < 1:997730/71‑1 (MQ=255) ctggtccggtCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCat < 1:2195101/58‑2 (MQ=255) gtctggtctgGAAGATGAACACTTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa < 1:2066342/67‑1 (MQ=255) gtctggtctgGAAGATGAACACGTTTCTTACAACGc < 1:2911476/36‑1 (MQ=255) gtccggtctgGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTgg < 1:3774248/70‑1 (MQ=255) | GACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGG > minE/834516‑834650 |
| 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: ‑ |