| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 834,827 | G→A | 10.0% | K1095K (AAG→AAA) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 834,827 | 0 | G | A | 10.0% | 118.5 / 3.0 | 40 | K1095K (AAG→AAA) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (16/20); new base A (2/2); total (18/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGG > minE/834761‑834894 | agTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTggg < 1:186127/71‑1 (MQ=255) aGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAt > 1:1689175/1‑69 (MQ=255) gCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccac < 1:182999/71‑1 (MQ=255) gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTggg < 1:2273653/59‑1 (MQ=255) gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc < 1:2474456/64‑1 (MQ=255) gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc < 1:3055214/64‑1 (MQ=255) aaaTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCCGAAGTGGGGTAtccactcca < 1:3052458/70‑1 (MQ=255) aaTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCAcc > 1:1555038/1‑71 (MQ=255) aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAg < 1:2562332/52‑1 (MQ=255) aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAt < 1:1423086/60‑1 (MQ=255) aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCt < 1:583933/71‑1 (MQ=255) cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTa > 1:2180752/1‑69 (MQ=255) cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTAc < 1:1052559/70‑1 (MQ=255) aCGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcc > 1:988776/1‑57 (MQ=255) aCGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGc > 1:3779071/1‑71 (MQ=255) ggCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGa > 1:2348567/1‑71 (MQ=255) gCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAg < 1:4141499/44‑1 (MQ=255) gCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAg < 1:2806355/44‑1 (MQ=255) cAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa < 1:1441309/71‑1 (MQ=255) cccggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAcc > 1:1675368/4‑71 (MQ=255) aGGAAAAAGCGCTCAACTTCCTGACGCCGAACCAGAAGTgg > 1:2864673/1‑41 (MQ=255) gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct < 1:1600553/71‑1 (MQ=255) aaaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa < 1:1109918/63‑1 (MQ=255) aaaaaGCACTTAACTTCCTGACGCCGCACCAGAAATggggca > 1:1583382/1‑40 (MQ=25) aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgctg > 1:2656707/1‑70 (MQ=255) aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctga > 1:3863982/1‑71 (MQ=255) aaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct < 1:1595360/68‑1 (MQ=255) aGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa < 1:1057953/59‑1 (MQ=255) aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc < 1:259185/71‑1 (MQ=255) aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc < 1:1257276/71‑1 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct < 1:2611229/65‑1 (MQ=255) cTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc < 1:350975/36‑1 (MQ=255) cTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCg < 1:2377482/52‑1 (MQ=255) cAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAAc > 1:1835006/1‑55 (MQ=255) aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt < 1:3982902/71‑1 (MQ=255) aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt > 1:3583517/1‑70 (MQ=255) cTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt < 1:1965311/68‑1 (MQ=255) ttCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAgg > 1:2110842/1‑71 (MQ=255) tGACGCCGCACCAGAAGTGGGGTATCCACTCCAcc < 1:2010715/35‑1 (MQ=255) aCGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcgg > 1:3810759/1‑71 (MQ=255) cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccggt > 1:3534572/1‑71 (MQ=255) gAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccggtgg > 1:3496604/1‑71 (MQ=255) aGTGGGGTATCCACTCCACCTACAGCGACAACctgc > 1:4118560/1‑36 (MQ=255) | AGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGG > minE/834761‑834894 |
| 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: ‑ |