| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,578,143 | 0 | A | C | 55.6% | 9.2 / 19.8 | 18 | L79V (TTA→GTA) | maeB | malic enzyme: putative oxidoreductase/phosphotransacetylase |
| Reads supporting (aligned to +/- strand): ref base A (5/3); new base C (0/10); total (5/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.54e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-03 | |||||||||||
AGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAG > NC_000913/2578001‑2578293 | aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAG‑c < 1:16686/150‑1 (MQ=255) cAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTaga > 2:108656/1‑151 (MQ=255) aaaaaaTTTGTCCGGGTCGATTTCGTCAACTTCAATGTCAAATCCATCAACCCCGGCGAATTTCTAAAACAGACCGCCCTTGCCTTCCATCCCCGGTTTGCCGGCCGGCGCGCCAATGTTGCCTACCCCCAG‑CACCGCCGTACCGTTagag < 2:276482/148‑1 (MQ=255) aaTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAACCCCGGCGAATTTCTAAAACAGAACGCCCTTGCCTCCCACCACCGGTTTGCCGGCCAGCGCGCAAATGTTGCCTACCCCCAG‑CACCGCCGTACCGTTagag < 2:99780/151‑1 (MQ=255) aTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTAAAATACATCAACCCCGGCGAATTTCTAAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTACCCCCAG‑CACCGCCGTACCGTTAgaga < 1:204125/151‑1 (MQ=255) ttaaaaTACATCAACCCCGGGGAATTTCTAAAACAGAACGCCCTTTCCTTCCATCCCCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg < 1:224107/148‑1 (MQ=255) ataaaatccacaaaTCCCGGCGAATTTCTAAAAAAAAACCCCCTTCCCTCCAATCCCGGGTTTGCCTGCGCGCGCGCCAAT‑TTGCCCAACCCCCGAAACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg < 1:106577/139‑1 (MQ=255) taaaaTACATCAACCCCGGCGAATTTCTTAAACAGACCGCCCTTCCCTTCCATCCCCGTTTTGCCTGCCAGCGCGCCAATGTTGCTTACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg < 1:149473/148‑1 (MQ=255) taaaaTACACCACCCCGGGCAATTTTCTAAAAAAGACCCCCTTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCCCCAATGTTCCCTACCCCCAG‑CGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg < 2:186701/148‑1 (MQ=255) tCAAATCCATCAACCCCGGCGAATTTCTAAAACAGACCGCCCTTGCCTCCCACCCCGGTTTTGCCGGCCGGCGCCCCAATGTGCCCTACCCCCAG‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg < 2:269837/150‑1 (MQ=255) aaTCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTaa > 1:151929/1‑150 (MQ=255) attttCTTAAACAGACCCCCCTTCCCTCCCACCCCCGTTTTCCCTGCCAGCGCCCCATTGTTCCCTACCCCCAG‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt < 2:221714/149‑1 (MQ=255) ttttGCCTCCCAGCGCGCCATTGTTCCCTACCCCCAC‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt < 1:123984/149‑1 (MQ=255) gcctgccAGCGCGCCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGc > 1:251962/1‑150 (MQ=255) gccAGCGCGCCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcg > 1:197876/1‑151 (MQ=255) cgcgcCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGa > 2:173659/1‑151 (MQ=255) cgcCAATGTTCCCTACCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGAt < 1:159050/150‑1 (MQ=255) cgcAAATGTTCCCTACCCCCAC‑CCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGAt < 2:306520/150‑1 (MQ=255) aCCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCag > 1:153579/1‑151 (MQ=255) | AGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAG > NC_000913/2578001‑2578293 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |