| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,578,143 | 0 | A | C | 51.9% | 46.3 / 82.5 | 52 | L79V (TTA→GTA) | maeB | malic enzyme: putative oxidoreductase/phosphotransacetylase |
| Reads supporting (aligned to +/- strand): ref base A (24/1); new base C (0/27); total (24/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.86e-14 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-09 | |||||||||||
TGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGC > NC_000913/2577994‑2578290 | tggttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTaa < 1:101113/150‑1 (MQ=255) tCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAAcccc > 2:808694/1‑150 (MQ=255) gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTGTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGc > 2:106680/1‑151 (MQ=255) cgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGt > 1:342207/1‑151 (MQ=255) aCAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCCCCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTACCCCCAGCACCGCCGt < 2:394911/149‑1 (MQ=255) aCAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTAc > 2:5116/1‑151 (MQ=255) aaCTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCg > 1:515816/1‑151 (MQ=255) ccGGGTCGATTTCGTCAACTTAAATGTAAAATACATCAACCCCGGCGAATTTCTTAAACAGAC‑CCCCTTTGCCTT‑CCATCCCCGTTTTG‑CCTGCCAGCGCGCAATTTTTGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGcca < 1:865800/150‑1 (MQ=255) cGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGccacc > 1:14530/1‑151 (MQ=255) ggttCGATTTCTCCAACTACAATGTCAAATACATCAACCCCGGCAATTTTCTAAAACAGAC‑CG‑CCTTGCCTTCCCACCCCCGGTTTG‑CCTGCCGGCGCCCAAGTGTTCCCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGccacca < 2:539919/148‑1 (MQ=255) ggTGGATTTCGTAAACTAAAATGTAAAATACACCACCCCCGGCGTTTTTCTAAAACAGAC‑CGCCCTTCCCTC‑CCACCCCGGGTTTG‑CCGGCCAGCGCGCAAATGTTCCCTACCCCAAGCACCGCCGTACCGTTAGAGATCACCGCCACCAg < 2:280115/151‑1 (MQ=255) aGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTAc > 2:546429/1‑151 (MQ=255) ttCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCt > 1:695723/1‑151 (MQ=255) cgacAACTTCAATGTAAAATCCATCACCCCGGGCGATTTTCTAAAACAGAC‑CCCCCTTCCCTC‑CCACCCCCGTTTTG‑CCGGCCAGCGCCCAATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCg < 2:100947/148‑1 (MQ=255) aaCTTAAATGTAAAATACATCAATCCGGGCGAATTTCTAAAACAGAA‑CGCCCTTGCCTC‑CCACCACCGGTTTG‑CCGGCCAGCGCGCAAATGTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCggg < 1:919812/149‑1 (MQ=255) aatgcaaaaTACATCACCCCGGGCAAATTTCTAAAACAGAC‑CCCCTTTCCCTC‑CCACCCCGGTTTTC‑CCGGCCGGCCCGCAAATGTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat < 2:233461/143‑1 (MQ=255) aaTGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat > 1:947109/1‑151 (MQ=255) gaaaaaTACATAACCCCCGGCAAATTTCTTAAACAACA‑CTCCCTCGCCGT‑CCACCACCGTTTCG‑CCTGCCAGCCCCCTAATGTTCCCTACCCCAACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt < 2:684259/147‑1 (MQ=255) cAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt > 2:448319/1‑150 (MQ=255) cATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt > 1:229130/1‑151 (MQ=255) aTCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt > 2:543755/1‑151 (MQ=255) tCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt > 2:907866/1‑150 (MQ=255) acAACCCCCGCCAATTTTTTAAAAAAAACCGCCCTTCCCTC‑CCACCCCCGTTTTG‑CCTCCCAGCGCGCAAATGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt < 2:695043/149‑1 (MQ=255) cAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCAGGCGGTATATTTGTAGGCTTTTa > 2:591314/1‑150 (MQ=255) acccccGGCAATTTTCTTAAACAAAC‑CCCCCTTCCCTC‑CCCTCACCGGTTTG‑CTTGCCCGCGCGCCATTGTTGCCTACCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtt < 1:722897/143‑1 (MQ=255) acccccGGC‑AATTTTTTAAACAAAA‑CCCCCTTTCCCT‑CCCACCCCCGTTTGCCCTCCAACCGCCCAATTTTTCCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt < 2:575686/145‑1 (MQ=255) aaTCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc > 1:588487/1‑151 (MQ=255) ccccGGCAATTTTCTAAAACAGAC‑CCCCTTTCCCTC‑CCTTCCCCGTTTTC‑CCTGCCGGCGCGCAAATTTTGCCTACCCCAAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt < 2:62523/144‑1 (MQ=255) cccGGCGATTTTCTAAAACAAAC‑CGCCCTTCCCTC‑CCATCCCCGTTTTG‑CCTCCCAGCGCCCCAATTTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg < 2:386333/151‑1 (MQ=255) attttCTAAAACAAAG‑CGCCTTTGCCTT‑CCAACCCCGGTTTG‑CAGGCCAGCACCCAATTTTTTCCTACCCCCACCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt < 1:717438/149‑1 (MQ=255) aaTTTCTAAAACAGAC‑CGCCTTTGCCTC‑CCATCCCCGGTTTG‑CCGGCCAGCCCGCCAATGTTCCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt < 2:376146/150‑1 (MQ=255) tttCTTAAACAGAC‑CGCCCTTGCTTC‑CCATCACCGTTTTG‑CCTCCCAGCGCGCCATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc < 1:540815/150‑1 (MQ=255) aagaaGAA‑CCCCCTCGCCAC‑CCAGCCCCGGTTGC‑TCTGCCGCCGACCCATGGTTGCTTAACCCCAGCTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt < 2:316417/143‑1 (MQ=255) aaacaaac‑CCCCTTTGCCTC‑CCATCCCGGTTTTG‑CCTGCCAGCGCGCCAATGTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa < 1:555256/151‑1 (MQ=255) aaacaaac‑CCCCCTTCCCTC‑CCACCCCCGTTTTG‑CCTGCCAGCGCGCAATTTTTCCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt < 1:877500/147‑1 (MQ=255) aaaCAGAC‑CGCCCTTGCCTT‑CCATCCCCGGTTTG‑CCGGCCGGCGCGCCATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa < 1:448706/151‑1 (MQ=255) aaaCAGAC‑CCCCCTTGCCTC‑CCACCCCCGTTTTG‑CCTGCCAGCGCGCAAATGTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa < 1:527604/151‑1 (MQ=255) aaaCAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc > 1:968829/1‑150 (MQ=255) aaaCAGAA‑CGCCCTTCCCTT‑CCATCCCCGTTTTG‑CCTCCCGGCGCGCCAATGTTCCCTACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAtt < 2:869260/148‑1 (MQ=255) gAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc > 1:210223/1‑150 (MQ=255) accGCCCTTGCCTT‑CCACCCCCGGTTTG‑CCGCCCAGCGCCCCATTGTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa < 2:914409/148‑1 (MQ=255) ccaaccccGGTTTG‑CTCGCCGGCCCGCCAAGTTTCCATACCCCGAGCACCGCCGTGCCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGcc < 1:793099/145‑1 (MQ=255) tCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTgg > 2:428210/1‑151 (MQ=255) ccGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg > 2:275623/1‑150 (MQ=255) ccGGTTTC‑CCTCCCGGCGCGCAAATGTTCCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg < 2:567691/150‑1 (MQ=255) ttgg‑cctgccAGCGCGCAAAGGTGGCCTCCCCCCACAGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTa < 2:451108/147‑1 (MQ=255) g‑cctgccAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc > 2:376151/1‑151 (MQ=255) ccAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcgc > 2:543145/1‑151 (MQ=255) cgcCAATTTTGCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGAt < 2:681198/150‑1 (MQ=255) ccATTGTTGCCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGa < 1:5116/147‑1 (MQ=255) aaTGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCgcgc > 2:129348/1‑151 (MQ=255) tGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCtg > 2:997235/1‑151 (MQ=255) cTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGc > 1:1004616/1‑151 (MQ=255) | TGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAA‑CGCCCTTGCCTT‑CCATCACCGGTTTG‑CCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGC > NC_000913/2577994‑2578290 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |