Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,435,822 | 0 | A | T | 26.7% | 100.0 / 43.5 | 60 | L50Q (CTG→CAG) | PP_1256 | alpha‑ketoglutarate semialdehyde dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (23/21); new base T (9/7); total (32/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCAT > NC_002947/1435673‑1435970 | gAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCa < 2:935625/150‑1 (MQ=255) gAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAAGGCCGCCCAGGCCa < 2:743503/150‑1 (MQ=255) gcgcCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTgcgcgc > 1:156290/1‑150 (MQ=255) gcCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTgcgcgctg > 2:323960/1‑148 (MQ=255) cTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGc > 2:891317/1‑150 (MQ=255) gggCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTg < 1:49692/150‑1 (MQ=255) gCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTggg < 1:323960/150‑1 (MQ=255) gCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGcggc < 2:156290/150‑4 (MQ=255) ccGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACg < 1:557060/150‑1 (MQ=255) cGGTTTCGATCAGGGCGCGGTCGATCAGCGCATCGCCGAGGGCCTCGATCTGTGTGGCGATTGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCACCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACgt < 2:735027/150‑1 (MQ=255) cGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACgt > 1:1136861/1‑150 (MQ=255) tttCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGAGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGtt > 1:707135/1‑150 (MQ=255) ggccgttCATGGCATCGCCTAGTGCCTCGTTCTGTGTGTCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCTAATGCCGCCCAGGCCATCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGa < 1:495058/144‑1 (MQ=255) ggTCGATGAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGa < 2:506281/150‑1 (MQ=255) acAGGGCCTCGCCGAGGGCCTCGATCTGTGTGGCGATGGCGTCAAGGAACTCTGCTCGTTGTTCGAGCGATTTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAgg > 2:317932/2‑150 (MQ=255) cAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGc < 1:408821/150‑1 (MQ=255) gggCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCg > 2:662349/1‑150 (MQ=255) tCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCt > 2:552130/1‑109 (MQ=255) tCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCt < 1:552130/109‑1 (MQ=255) aGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTc < 2:567204/150‑1 (MQ=255) ccTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGAc < 1:486775/150‑1 (MQ=255) cTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGTCCAGCGCGAAGTCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGAcc > 2:171747/1‑150 (MQ=255) cTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGAcc > 1:301252/1‑150 (MQ=255) aTCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCggtg > 2:687037/1‑150 (MQ=255) tgGCGATGGCTTCAAGGCCCCCGGCTCGTTGGTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGt < 1:386107/150‑1 (MQ=255) tgGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGt > 1:1179767/1‑150 (MQ=255) gATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGAt > 2:793229/1‑150 (MQ=255) gATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGAt < 1:662349/150‑1 (MQ=255) ttCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACg > 1:1079570/1‑150 (MQ=255) tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACgg > 1:62089/1‑150 (MQ=255) aaGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgcgc > 1:474759/1‑66 (MQ=255) aaGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgcgc < 2:474759/66‑1 (MQ=255) tctGCTCACTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTcg > 1:1058583/1‑150 (MQ=255) ctGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTcgc < 1:891317/150‑1 (MQ=255) ctGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTcgc < 1:687037/150‑1 (MQ=255) gttgttCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGc > 1:376649/1‑150 (MQ=255) gttCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCgg < 1:793229/150‑1 (MQ=255) ttCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGt > 1:1234857/1‑150 (MQ=255) gAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGtt < 2:449746/112‑1 (MQ=255) gAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGtt > 1:449746/1‑112 (MQ=255) gCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTggg < 1:668611/54‑1 (MQ=25) gCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTggg > 2:668611/1‑54 (MQ=25) tttCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCg > 2:784310/1‑144 (MQ=255) tttCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCATGTGCCGTCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCg < 1:784310/144‑1 (MQ=255) gTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGAt < 2:176173/150‑1 (MQ=255) cATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagca > 1:494610/1‑150 (MQ=255) aTCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagcag > 1:377438/1‑150 (MQ=255) tCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAgg > 2:80035/1‑150 (MQ=255) tGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGtgcca < 2:1179767/150‑5 (MQ=255) gccgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCt < 2:376649/150‑1 (MQ=255) gccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCtgt > 1:1218952/1‑149 (MQ=255) gccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCtgt < 2:1218952/149‑1 (MQ=255) ccAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCgg < 1:727543/57‑1 (MQ=37) ccAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCgg > 2:727543/1‑57 (MQ=37) ccAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTc > 2:731195/1‑84 (MQ=255) ccAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTc < 1:731195/84‑1 (MQ=255) ccAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGa < 1:938807/131‑1 (MQ=255) ccAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAg > 1:309631/1‑150 (MQ=255) ccAGGCCTGAGCGCTGGCATGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGa > 2:938807/1‑131 (MQ=255) cAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGCTGGGAACGGGTCGGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAgg > 2:980127/1‑150 (MQ=255) cAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAgg < 1:80035/150‑1 (MQ=255) cAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCAt > 2:1133970/1‑150 (MQ=255) | GAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCAT > NC_002947/1435673‑1435970 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |