| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 834,827 | G→A | 34.6% | 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 | 34.6% | 85.1 / 65.9 | 78 | K1095K (AAG→AAA) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (28/23); new base A (12/15); total (40/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.77e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCT > minE/834760‑834897 | aagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCT‑AACTTCCTGACGCCGCACCAGAAGTggg > 1:4676859/1‑71 (MQ=255) gTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTgggg < 1:13723/71‑1 (MQ=255) gTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTgggg < 1:1587241/71‑1 (MQ=255) cAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactc > 1:2093182/1‑71 (MQ=255) aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCAcc > 1:4653409/1‑70 (MQ=255) cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTa > 1:5375499/1‑69 (MQ=255) aaCGGCTTCCCGGAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATggggca < 1:173406/55‑3 (MQ=255) cGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAg > 1:2368751/1‑69 (MQ=255) cGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGc > 1:4637303/1‑70 (MQ=255) ggCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGa > 1:2051331/1‑71 (MQ=255) cAACCAGGAAAAAGCGCTCAACTTTCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa < 1:3852092/71‑1 (MQ=255) aaCCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa < 1:2695136/70‑1 (MQ=255) ccggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAcc < 1:4368715/68‑1 (MQ=255) cAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACaa > 1:1402299/1‑68 (MQ=255) cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc < 1:3157486/67‑1 (MQ=255) cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc < 1:567457/67‑1 (MQ=255) cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc < 1:189225/67‑1 (MQ=255) aGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctg > 1:3144394/1‑71 (MQ=255) gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcca > 1:2401833/1‑47 (MQ=255) gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct < 1:3603640/71‑1 (MQ=255) gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct < 1:1125150/71‑1 (MQ=255) gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct < 1:916451/71‑1 (MQ=255) gAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTAGAGTGAAAACct > 1:2600314/1‑68 (MQ=255) gAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAgtg < 1:2542890/61‑3 (MQ=255) aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctga > 1:4887880/1‑71 (MQ=255) aaaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCAttcact > 1:2436063/1‑46 (MQ=21) aaaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaat > 1:4634689/1‑71 (MQ=255) aaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc < 1:666464/41‑1 (MQ=255) aaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatg < 1:3896675/71‑1 (MQ=255) aaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatg < 1:5477393/71‑1 (MQ=255) aGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcca > 1:3706156/1‑42 (MQ=255) aGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcca > 1:1472731/1‑42 (MQ=255) aGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcca > 1:4597606/1‑42 (MQ=255) aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc > 1:305566/1‑71 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac > 1:5217717/1‑47 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac > 1:5211438/1‑47 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac > 1:839167/1‑47 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac > 1:1098225/1‑47 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactcca > 1:1030869/1‑46 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct < 1:4230928/71‑1 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct < 1:3647857/71‑1 (MQ=255) gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct < 1:778149/71‑1 (MQ=255) cgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct > 1:3873271/1‑70 (MQ=255) cttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAc > 1:132462/4‑71 (MQ=255) cttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAc > 1:3654334/4‑71 (MQ=255) ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACcttct < 1:2322660/59‑1 (MQ=255) ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcg > 1:2673952/3‑70 (MQ=255) ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcg > 1:2037702/3‑70 (MQ=255) aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGCCtt > 1:456030/1‑70 (MQ=255) aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt < 1:2071783/71‑1 (MQ=255) aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt < 1:3100376/71‑1 (MQ=255) aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt < 1:5208427/71‑1 (MQ=255) aaCTTCCTGACGCCGCACCAGAAATGGGGCAttcactc > 1:236988/1‑38 (MQ=25) aaCTTCCTGACGCCGCACCAGAAATGGGGCAttcactc > 1:1503412/1‑38 (MQ=25) aaCTTCCTGACGCCGCACCAGAAATGGGGCAttcactc > 1:5025560/1‑38 (MQ=25) aaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct > 1:270191/1‑68 (MQ=255) aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt > 1:1076048/1‑69 (MQ=255) aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt > 1:1695595/1‑69 (MQ=255) cTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAg > 1:5416809/1‑71 (MQ=255) tCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGGCa < 1:2015197/48‑1 (MQ=255) tCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGt > 1:3679281/1‑71 (MQ=255) tCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGaaa < 1:4457032/48‑3 (MQ=19) cTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTcg < 1:1511231/71‑1 (MQ=255) gACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgctaa < 1:3006596/53‑3 (MQ=255) aCGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcgg < 1:4046506/71‑1 (MQ=255) aCGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct < 1:1805862/62‑4 (MQ=255) cgccgcACCAGAAATGGGGCATTCAGTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct < 1:1920594/61‑4 (MQ=255) gccgcACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgct < 1:1119767/48‑1 (MQ=19) gccgcACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgct < 1:1606080/48‑1 (MQ=19) ccgcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtgg < 1:4735454/71‑1 (MQ=255) cgcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggt > 1:2559416/1‑71 (MQ=255) gcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtc > 1:4439869/1‑71 (MQ=255) aCCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccg < 1:368621/71‑1 (MQ=255) cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg < 1:5671764/45‑1 (MQ=255) cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg < 1:3101428/45‑1 (MQ=255) cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg < 1:434018/45‑1 (MQ=255) aGAAGTAGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcgg < 1:3939067/61‑1 (MQ=255) gAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAgg > 1:2987148/1‑54 (MQ=255) gTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCt > 1:1825984/1‑71 (MQ=255) | AAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCT > minE/834760‑834897 |
| 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: ‑ |