Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,578,0680TC53.8% 9.6 / 14.2 13I104V (ATT→GTT) maeBmalic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base C (0/7);  total (2/11)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02

CGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAA  >  NC_000913/2577926‑2578215
                                                                                                                                              |                                                                                                                                                    
cgcAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCg                                                                                                                                               <  2:203867/150‑1 (MQ=255)
                                       aTGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGtt                                                                                                       <  2:194391/151‑1 (MQ=255)
                                                                           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)
                                                                                                                                              |                                                                                                                                                    
CGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAG‑CACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAA  >  NC_000913/2577926‑2578215

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑