| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 834,575 | G→C | 27.0% | P1011P (CCG→CCC) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 834,575 | 0 | G | C | 27.0% | 93.3 / 40.2 | 63 | P1011P (CCG→CCC) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (31/15); new base C (8/9); total (39/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.57e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01 | |||||||||||
GAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATC > minE/834506‑834644 | gAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGa > 1:1145938/1‑71 (MQ=255) aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCAc > 1:3447536/1‑67 (MQ=255) aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctgg > 1:1272132/1‑71 (MQ=255) aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctgg > 1:2264364/1‑71 (MQ=255) aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctgg < 1:946677/71‑1 (MQ=255) aagaagACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctgg > 1:1483700/1‑71 (MQ=255) agaggaCGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctggt > 1:3512271/1‑71 (MQ=255) gaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtc > 1:1206579/1‑71 (MQ=255) 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) ttCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctgg > 1:3934788/1‑51 (MQ=255) tCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa > 1:3157990/1‑68 (MQ=255) tCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac > 1:162863/1‑71 (MQ=255) ccGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac < 1:250499/70‑1 (MQ=255) ccGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac > 1:1261589/1‑70 (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:653439/1‑68 (MQ=255) cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa > 1:3752800/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:3339452/61‑1 (MQ=255) ttATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATa < 1:2956570/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) | GAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATC > minE/834506‑834644 |
| 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: ‑ |