Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,465,7870GA48.4% 27.4 / 11.2 31R132C (CGC→TGC) treBtrehalose‑specific PTS enzyme: IIB and IIC component
Reads supporting (aligned to +/- strand):  ref base G (15/1);  new base A (0/15);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 5.66e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTT  >  NC_000913/4465651‑4465935
                                                                                                                                        |                                                                                                                                                    
cAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCa                                                                                                                                        >  2:322651/1‑151 (MQ=255)
      aTACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAAccg                                                                                                                                   >  1:477521/1‑150 (MQ=255)
                       gaagaagaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCggg                                                                                                                  >  2:476820/1‑150 (MQ=255)
                                              aCCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACcg                                                                                           >  2:190365/1‑150 (MQ=255)
                                                        aTCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag                                                                                >  2:588266/1‑151 (MQ=255)
                                                               aTCGTTTGCAGGGAAGGGTACATTGGCGCCAGCGTCTGCCCGTTGCTCAGGGGAAAATCGCCGATCACATTGCAAAACCCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAAc                                                                          <  2:9320/150‑1 (MQ=255)
                                                                tCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCaa                                                                           >  2:192153/1‑148 (MQ=255)
                                                                tCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCAAAACCCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTg                                                                        <  2:184742/151‑1 (MQ=255)
                                                                   tttGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCt                                                                      >  2:207430/1‑150 (MQ=255)
                                                                      gCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATg                                                                  >  1:618185/1‑151 (MQ=255)
                                                                         ggtaaGGATACATTTGCGCCGGCGTCTGACCTTTGCTCGTGGTCAACTCGCCGTTCAATTTGCAAAACCCGAGGAAAAATCCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGcc                                                                <  2:39563/147‑1 (MQ=255)
                                                                         gggAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCAAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGcc                                                                <  2:113211/150‑1 (MQ=255)
                                                                              gggTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCa                                                          >  2:95124/1‑151 (MQ=255)
                                                                                gTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCata                                                        >  2:657326/1‑151 (MQ=255)
                                                                                 tACATTTGCGCCAGCGTCTGCCCGTTGCTCATGGGAAAATCGCCGATCACATTGCAAAACCCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCatat                                                       <  2:530440/151‑1 (MQ=255)
                                                                                 tACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCAAAACCCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCatat                                                       <  2:462624/151‑1 (MQ=255)
                                                                                 tACATTGGCGCCAGCGTCTGACCGTTGCTCATGGGAAAACCGCCGATCACATTGAAAAACCCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCatat                                                       <  1:341508/151‑1 (MQ=255)
                                                                                 tACATTGCCGCCAGCGTCTGACCGTTCCTCAGGGCAAAACCCCCGATCACATTGCAAAAACCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCatat                                                       <  1:481080/151‑1 (MQ=255)
                                                                                      ttGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAACCCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGc                                                  <  2:132096/151‑1 (MQ=255)
                                                                                                          tGCTCATGGGAAAATCGCCGATCACATTGCAAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTg                              <  2:457715/151‑1 (MQ=255)
                                                                                                           gCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTg                              >  1:340723/1‑150 (MQ=255)
                                                                                                                   gaaaaCCGCCGATCACATTGCAAAAACCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTAt                      <  2:231106/148‑1 (MQ=255)
                                                                                                                   gCAAATCGCCGATCACATTGCAAAAACCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTAt                      <  2:209072/150‑1 (MQ=255)
                                                                                                                    cAAATCGCCGATCAAATTGCAAAACCCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATc                     <  2:92763/150‑1 (MQ=255)
                                                                                                                        tCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAAc                  >  1:452067/1‑149 (MQ=255)
                                                                                                                            cGATCACATTGCAAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTgcgc            <  2:274086/151‑1 (MQ=255)
                                                                                                                              aTCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTgcgc            >  2:261737/1‑149 (MQ=255)
                                                                                                                              aTCACATTGCAAAACCCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGcc           <  2:54110/150‑1 (MQ=255)
                                                                                                                              aTCACATTGCAAAAACCGAGAACAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGcc           <  1:454193/150‑1 (MQ=255)
                                                                                                                                      gCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGtt  >  2:382403/1‑151 (MQ=255)
                                                                                                                                      gCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGtt  >  2:35692/1‑151 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
CAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTT  >  NC_000913/4465651‑4465935

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