Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,628,6250AC53.8% 5.5 / 14.3 14V105G (GTC→GGC) yfgFcyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base C (0/7);  total (4/10)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-03

GCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAAT  >  NC_000913/2628486‑2628747
                                                                                                                                           |                                                                                                                          
gCTCGCCAGTAATCCTACAAATGCAGAAAACTGGAAAAGTATCAAAAACAGCGTTGCAGGAAACGCAATCTGCCAGAATATGCGCTAGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTGCGCCGGGGAATTAAGGCCCTGTa                                                                                                                     <  2:383409/147‑1 (MQ=255)
     ccAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGcc                                                                                                               >  1:183364/1‑148 (MQ=255)
     ccAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGcc                                                                                                               <  2:183357/148‑1 (MQ=255)
          aaTCCTACAAATGCAGCAAACGGGAAAAGTATCAGAAACAGCGTGGCAGGAAACACAATCTGCCAGAATATACGCGGGGAAATAAACCGGGTACCCCCATGTAAACCGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACag                                                                                                        <  2:199934/150‑1 (MQ=255)
              cTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaat                                                                                                   >  1:160704/1‑151 (MQ=255)
                         gcaAACGGGAAAAATATCAGAAACAGTGTTGCAGGAAACACAATCGCCCAAAATATCCGCTGGGAAATTAAACGGTTATCCCCATGTGAAACGTTATTACGCCGGGAAAAAAAGCCCCAGTAGCCGCCCCAACAGAGAAtgataatgataa                                                                                        <  1:132160/151‑1 (MQ=255)
                          caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATgataatgataat                                                                                       >  1:740199/1‑151 (MQ=255)
                           aaaCGGAAAAAGTATCAGAAACAGCGTGGCAGGAAACACAATCTGCCAGAATATCCGCTGGAAAATAAACCGGGTACCCCCAGGTGAAACGTTTTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGaatgataatg                                                                                            <  2:621714/145‑1 (MQ=255)
                              cTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaa                                                                                    >  1:689076/1‑150 (MQ=255)
                                           aGAAACAGCGTTGAAGAAAACACAATCTGCCAAAATATACGCGGGGAAATTAAACGGGTATCACCATGTAAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTa                                                                         <  1:242477/148‑1 (MQ=255)
                                                  gctttGCAGGAAACACAATCTGCAAAAATATCCGCTGGAAAATTAAACGGTTATCCCCTTGTAAAACTTTTTTCCGCCGGGGATTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc                                                                 <  1:2899/146‑1 (MQ=255)
                                                                            aaTATACGCTGGGAAATTAAACGGGTATCACCATGTAAAACTTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCat                                        <  1:518305/148‑1 (MQ=255)
                                                                                                      ataaCCATGTGAAACTTTGTTCCGCCGGGAAGAAAAGCCCCGGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGAt            <  2:427946/147‑1 (MQ=255)
                                                                                                                 aaaCGTTGTCACGCCGGGGATTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAAt  <  2:701635/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                          
GCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAAT  >  NC_000913/2628486‑2628747

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