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

CGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAA  >  NC_000913/2577926‑2578215
                                                                                                                                              |                                                                                                                                                    
cgcAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCg                                                                                                                                               <  2:203867/150‑1 (MQ=255)
                                       aTGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGtt                                                                                                       <  2:194391/151‑1 (MQ=255)
                                                                           aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGc                                                                   <  1:16686/150‑1 (MQ=255)
                                                                                           cAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTaga                                                  >  2:108656/1‑151 (MQ=255)
                                                                                            aaaaaaTTTGTCCGGGTCGATTTCGTCAACTTCAATGTCAAATCCATCAACCCCGGCGAATTTCTAAAACAGACCGCCCTTGCCTTCCATCCCCGGTTTGCCGGCCGGCGCGCCAATGTTGCCTACCCCCAGCACCGCCGTACCGTTagag                                                 <  2:276482/148‑1 (MQ=255)
                                                                                            aaTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAACCCCGGCGAATTTCTAAAACAGAACGCCCTTGCCTCCCACCACCGGTTTGCCGGCCAGCGCGCAAATGTTGCCTACCCCCAGCACCGCCGTACCGTTagag                                                 <  2:99780/151‑1 (MQ=255)
                                                                                             aTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTAAAATACATCAACCCCGGCGAATTTCTAAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTACCCCCAGCACCGCCGTACCGTTAgaga                                                <  1:204125/151‑1 (MQ=255)
                                                                                                                                ttaaaaTACATCAACCCCGGGGAATTTCTAAAACAGAACGCCCTTTCCTTCCATCCCCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg             <  1:224107/148‑1 (MQ=255)
                                                                                                                                ataaaatccacaaaTCCCGGCGAATTTCTAAAAAAAAACCCCCTTCCCTCCAATCCCGGGTTTGCCTGCGCGCGCGCCAATTTGCCCAACCCCCGAAACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg             <  1:106577/139‑1 (MQ=255)
                                                                                                                                 taaaaTACATCAACCCCGGCGAATTTCTTAAACAGACCGCCCTTCCCTTCCATCCCCGTTTTGCCTGCCAGCGCGCCAATGTTGCTTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg             <  1:149473/148‑1 (MQ=255)
                                                                                                                                 taaaaTACACCACCCCGGGCAATTTTCTAAAAAAGACCCCCTTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCCCCAATGTTCCCTACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg             <  2:186701/148‑1 (MQ=255)
                                                                                                                                 tCAAATCCATCAACCCCGGCGAATTTCTAAAACAGACCGCCCTTGCCTCCCACCCCGGTTTTGCCGGCCGGCGCCCCAATGTGCCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg             <  2:269837/150‑1 (MQ=255)
                                                                                                                                            aaTCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTaa  >  1:151929/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAA  >  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: