| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 832,469 | G→A | 7.6% | A309A (GCG→GCA) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 832,469 | 0 | G | A | 7.6% | 228.6 / 6.2 | 79 | A309A (GCG→GCA) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (20/53); new base A (3/3); total (23/56) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.50e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01 | |||||||||||
CTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTC > minE/832399‑832539 | cTGAATTCGCCAAAGTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCg < 1:2510370/71‑1 (MQ=255) tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa > 1:4293305/1‑71 (MQ=255) aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa < 1:30606/69‑1 (MQ=255) aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt < 1:2656352/70‑1 (MQ=255) aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg > 1:4135535/1‑71 (MQ=255) aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg < 1:4550437/71‑1 (MQ=255) aTCGCCAAACTGTGCGATCTGTGGCTGGCCCCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggt > 1:1041489/1‑70 (MQ=255) aTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggg > 1:1459041/1‑71 (MQ=255) gCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCa > 1:5162204/1‑71 (MQ=255) ccAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCat < 1:4956967/71‑2 (MQ=255) cAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCatg > 1:972398/1‑69 (MQ=255) aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt < 1:3292348/71‑1 (MQ=255) aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt < 1:4691403/71‑1 (MQ=255) aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt < 1:647448/71‑1 (MQ=255) tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt > 1:1595232/1‑70 (MQ=255) cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa > 1:5161605/1‑49 (MQ=255) tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg < 1:4393558/67‑1 (MQ=255) tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa < 1:3415584/68‑1 (MQ=255) ggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg < 1:5047364/66‑1 (MQ=255) gctggcACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATg > 1:3183114/1‑46 (MQ=255) aCCGAAACAGGGCTCCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa < 1:545220/60‑1 (MQ=255) ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt < 1:3193421/47‑1 (MQ=255) ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt < 1:2098150/47‑1 (MQ=255) ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt < 1:689026/47‑1 (MQ=255) ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc < 1:2653088/69‑1 (MQ=255) cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGa > 1:991128/1‑70 (MQ=255) cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc < 1:2245183/71‑1 (MQ=255) gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtg > 1:1505257/1‑49 (MQ=255) gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc > 1:2577609/1‑70 (MQ=255) aaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACaa > 1:5637108/1‑70 (MQ=255) aCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa < 1:4273547/54‑1 (MQ=255) cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc < 1:2510651/66‑1 (MQ=255) cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAccc < 1:3708554/71‑1 (MQ=255) aGGGCACCGATGCGTCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc < 1:3677352/65‑1 (MQ=255) ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc < 1:2955774/46‑1 (MQ=255) ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc < 1:4930043/46‑1 (MQ=255) ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc < 1:3481310/46‑1 (MQ=255) ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc < 1:5364858/46‑1 (MQ=255) gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc < 1:903663/71‑1 (MQ=255) cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg < 1:249938/47‑1 (MQ=255) cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg < 1:3710010/47‑1 (MQ=255) cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg < 1:3703196/47‑1 (MQ=255) ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg < 1:116302/45‑1 (MQ=255) ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg < 1:1822337/45‑1 (MQ=255) ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg < 1:860809/71‑1 (MQ=255) ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg < 1:1412590/71‑1 (MQ=255) ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGATGCGTGAATTCCACCTCGACAACCCAAGCCAg < 1:1297774/71‑1 (MQ=255) cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg > 1:2784953/1‑70 (MQ=255) cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGt < 1:4327533/71‑1 (MQ=255) cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGt < 1:5093007/71‑1 (MQ=255) gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc < 1:4033357/54‑1 (MQ=255) gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg < 1:5112555/69‑1 (MQ=255) aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt < 1:1572162/71‑1 (MQ=255) aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt < 1:481577/71‑1 (MQ=255) aTGCGGCAATGGCGCTGGCGATGGGCCACCTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt < 1:2553456/71‑1 (MQ=255) ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc < 1:2450968/49‑1 (MQ=255) ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc < 1:4976215/49‑1 (MQ=255) ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa < 1:4598708/71‑1 (MQ=255) ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa < 1:400012/71‑1 (MQ=255) cAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAttt > 1:4337798/1‑67 (MQ=255) cAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCAAGTAttt > 1:4649254/1‑67 (MQ=255) aaTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc < 1:868249/49‑1 (MQ=255) aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCaa < 1:3529823/55‑1 (MQ=255) aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc < 1:805254/58‑1 (MQ=255) aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGa < 1:1758198/71‑1 (MQ=255) tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa < 1:1772831/48‑1 (MQ=255) tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa < 1:160473/48‑1 (MQ=255) tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa < 1:2770755/48‑1 (MQ=255) tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCTAc > 1:4360562/1‑71 (MQ=255) gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACt < 1:5683136/70‑1 (MQ=255) gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTa < 1:2298415/71‑1 (MQ=255) gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTa < 1:1635669/71‑1 (MQ=255) gcgcTGGCGATGGGCCAAGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAc > 1:4827672/1‑65 (MQ=255) cgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAt < 1:5069339/71‑1 (MQ=255) tGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAtgtg > 1:5176594/1‑71 (MQ=255) ggCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc > 1:1222603/1‑71 (MQ=255) gCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCg > 1:5610630/1‑71 (MQ=255) gCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCg > 1:3295924/1‑71 (MQ=255) cGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCTACTATGTGCGt < 1:1811330/71‑1 (MQ=255) cGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGt > 1:210657/1‑71 (MQ=255) gATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTc < 1:4591265/71‑1 (MQ=255) | CTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTC > minE/832399‑832539 |
| 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: ‑ |