Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,578,143 | 0 | A | C | 42.6% | 94.4 / 15.4 | 54 | L79V (TTA→GTA) | maeB | malic enzyme: putative oxidoreductase/phosphotransacetylase |
Reads supporting (aligned to +/- strand): ref base A (28/3); new base C (0/23); total (28/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.39e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC > NC_000913/2577996‑2578291 | gttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAAccc > 2:335726/1‑151 (MQ=255) aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGc > 1:846111/1‑150 (MQ=255) aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGc > 1:750329/1‑150 (MQ=255) gcggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCAccg < 2:443064/151‑1 (MQ=255) cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCAccgc < 1:726333/151‑1 (MQ=255) acacCTTAAATAAATTTGTCCGGGTCGAGTTCGTCAACTCC‑AATTACAAATACATCAACCCCGGCAAA‑TTTCTTAAA‑CAGAC‑CCCCCTTGCCTC‑CCACCCCCGGTTTGCCTG‑CC‑AGCGCCCCATTGTTCCCTACCCCCAGCACCGCCGTa < 2:153142/150‑1 (MQ=255) aCAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTAc > 1:565437/1‑151 (MQ=255) cTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGtt > 2:261675/1‑151 (MQ=255) aaTAAATTTGCCCGGGTCGATTTCGTCACCTTA‑AATGTAAAATACATCAACCCCGGCAAA‑TTTCTAAAA‑CAGAC‑CCCCCTTCCCTC‑CCACCCCCGTTTTGCCTG‑CC‑AGCGCGCAAATGTTGCCTACCCCCAGCACCGCCGTACCGTTagag < 2:126491/151‑1 (MQ=255) tAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTag > 2:656915/1‑147 (MQ=255) ccGGGTCGATTTCGTCATCTTC‑AATGCAAAATACATCAA‑CCCCGCGAATTTTCTTAAA‑AAAAC‑CGCCTTTGCCTT‑CCACCCCGGGTTTGCCGG‑CC‑AGCCCCCCAATGTGGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGcc < 1:203108/149‑1 (MQ=255) tCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGt > 1:182941/1‑151 (MQ=255) ttCGTCAACTTC‑AATTTAAAATACATCAACCCCGGCGAT‑TTTCTAAAA‑CAAAC‑CGCCTTGCCCTC‑CCACCCCCGGTTTGCCTG‑CC‑AGCGCGACAATGTTCCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTa < 1:766575/148‑1 (MQ=255) tCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGCTACCTc > 2:61988/1‑151 (MQ=255) cgtaaCTTCAAATGCAAAATACATCACCCCCGGCGAA‑TTTCTAAAA‑CAAAC‑CGCCCTTCCCTC‑CCACCCCGGTTTTGCCTG‑CC‑AGCGCCCAAATGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTc < 2:251908/147‑1 (MQ=255) cAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc > 1:329515/1‑151 (MQ=255) aaCTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCg > 1:527912/1‑151 (MQ=255) taaaaTGTCAAATACACCAA‑CCCGGGCAA‑TTTCCTAAACCAGAC‑CGCCCTTCCTTC‑CCACCCCCGGTTTGCCTC‑CC‑GGCGGCCCAGTGTTCCCACCCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCgg < 2:263828/146‑1 (MQ=255) actctaaaaTACAACAATCCGGGGGAT‑TTTCTAAAA‑CAGAG‑CCCCCTTGCCTC‑CCATCACTGTTGTGCCGA‑CG‑AGGGCGCTAATGTTGCCTACCCCCATCTCCGCTGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat < 2:106143/143‑1 (MQ=255) gTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg > 1:471370/1‑151 (MQ=255) cAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt > 1:610310/1‑150 (MQ=255) aaaTACATCAACCCGGGCAAT‑TTCCTAAAA‑CAGAA‑CCCCCTTCCCTC‑CCACCCCGGGTTGGCCTG‑CC‑GGCGCCCCAATTTTCCCAACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt < 2:278147/147‑1 (MQ=255) aaTACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAgg > 1:550240/1‑151 (MQ=255) acccccGGCAAT‑TTTCTAAAC‑CAGAA‑CGCCCTTCCCTC‑CCACCCCCGTTTTGCCGC‑CC‑AGCGCGCCAATGTTGCCTACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt < 2:53558/144‑1 (MQ=255) aaccccGGCAAA‑TTTCTTAAA‑AAAAC‑CCCCTTTGCCTC‑CCACCCCCGTTTTGCCTG‑CC‑AGCGCGCAAATGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt < 1:397668/145‑1 (MQ=255) aTCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg > 1:70957/1‑151 (MQ=255) ccccGGCAAA‑TTTCTAAAA‑CAG‑A‑CCCCCTTTCCTTCCCATCCCCGTTTTCCCTC‑CC‑AGCGCGCAATTTTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACgg < 2:89602/150‑1 (MQ=255) cccGGCAAA‑TTTCTAAAA‑CAAAC‑CGCCCTTGCCTT‑CCACCCCCGGTTGGCCTG‑CC‑AGCGCCCAAATGTTGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg < 1:586182/151‑1 (MQ=255) ccGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGt > 2:136254/1‑151 (MQ=255) ccaattttATTAAA‑CATAC‑CCCCCTCTCCTT‑CCATCACCTTTTTCCCGC‑CC‑AGCGGGCAATCGTAGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTc < 2:352239/144‑1 (MQ=255) caa‑TTTTTTAAA‑CAAAACCCCCCTTCCCTT‑CCA‑CCCCGGTTTTCCCGCCC‑GGCCCGCAAATGTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt < 1:535790/149‑1 (MQ=255) attttCTTAAA‑CAGAC‑CGCCCTTCCCTC‑CCACCCCCGTTTTGCCTG‑CC‑AGCGCGCCATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt < 2:597988/148‑1 (MQ=255) tttttaaaccaaaa‑CGCCGTGGCCTC‑CCACCCCCGG‑TTGCCTT‑CCGGGGGCGCCAATTTTGCCAACCCCCAGCTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt < 1:649779/137‑1 (MQ=255) ttCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg > 1:831417/1‑150 (MQ=255) aaa‑CAGAC‑CGCCTTTGCCTC‑CCACCCCCGTTTTGCCTC‑CC‑AGCGCGCAAATTTTGCCACCCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa < 2:629442/151‑1 (MQ=255) a‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa > 1:377796/1‑149 (MQ=255) aGAC‑CGCCCTTCCCTT‑CCATCACCGTTTTGCCTG‑CC‑AGCGCGCCAATGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGa < 2:344352/150‑1 (MQ=255) aacccccctgcccTT‑CCATCCCCGTTTTGCCTC‑CC‑AGCGCGCAAATGTTCCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa < 1:531329/140‑1 (MQ=255) a‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa > 1:760732/1‑149 (MQ=255) a‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa > 2:891161/1‑149 (MQ=255) gCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg > 2:501072/1‑150 (MQ=255) cccTTCCCTC‑CCACCCCCGTTTTGCCTG‑CC‑AGCGCCCCAATGTTGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt < 1:359629/150‑1 (MQ=255) gCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGc > 2:324371/1‑150 (MQ=255) ccTC‑CCACCCCCGGTTGCCCCC‑CC‑AGCGCGCAAATGTTCCATACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGc < 2:597208/149‑1 (MQ=255) cATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCt > 1:478516/1‑151 (MQ=255) ccGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg > 1:717913/1‑150 (MQ=255) cGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg > 2:914997/1‑149 (MQ=255) ggTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt > 2:696202/1‑151 (MQ=255) gTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt > 2:70602/1‑150 (MQ=255) ttGCCTC‑CC‑AGCCCGCCATTGTTCCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAg < 2:283862/150‑1 (MQ=255) ccggcccgcaattGTTCCCAACCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcg < 2:410431/138‑1 (MQ=255) gcgcgcCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc > 2:431551/1‑149 (MQ=255) attGTTGCCTACCCCCACCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCgcgc < 2:921807/149‑1 (MQ=255) cTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGc > 2:649723/1‑151 (MQ=255) aaCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGcc > 1:365470/1‑150 (MQ=255) | GTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTC‑AATGTCAAATACATCAATCCCGGCGAA‑TTTCTTAAA‑CAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTGCCTG‑CC‑AGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC > NC_000913/2577996‑2578291 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |