Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,578,1430AC42.6% 94.4 / 15.4 54L79V (TTA→GTA) maeBmalic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (28/3);  new base C (0/23);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 2.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2577996‑2578291
                                                                                                                                                          |                                                                                                                                                    
gttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAAccc                                                                                                                                                   >  2:335726/1‑151 (MQ=255)
     aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGc                                                                                                                                               >  1:846111/1‑150 (MQ=255)
     aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGc                                                                                                                                               >  1:750329/1‑150 (MQ=255)
        gcggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAccg                                                                                                                                           <  2:443064/151‑1 (MQ=255)
         cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAccgc                                                                                                                                          <  1:726333/151‑1 (MQ=255)
              acacCTTAAATAAATTTGTCCGGGTCGAGTTCGTCAACTCCAATTACAAATACATCAACCCCGGCAAATTTCTTAAACAGACCCCCCTTGCCTCCCACCCCCGGTTTGCCTGCCAGCGCCCCATTGTTCCCTACCCCCAGCACCGCCGTa                                                                                                                                      <  2:153142/150‑1 (MQ=255)
              aCAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTAc                                                                                                                                     >  1:565437/1‑151 (MQ=255)
                  cTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGtt                                                                                                                                 >  2:261675/1‑151 (MQ=255)
                      aaTAAATTTGCCCGGGTCGATTTCGTCACCTTAAATGTAAAATACATCAACCCCGGCAAATTTCTAAAACAGACCCCCCTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCGCAAATGTTGCCTACCCCCAGCACCGCCGTACCGTTagag                                                                                                                             <  2:126491/151‑1 (MQ=255)
                        tAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTag                                                                                                                               >  2:656915/1‑147 (MQ=255)
                                 ccGGGTCGATTTCGTCATCTTCAATGCAAAATACATCAACCCCGCGAATTTTCTTAAAAAAACCGCCTTTGCCTTCCACCCCGGGTTTGCCGGCCAGCCCCCCAATGTGGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGcc                                                                                                                    <  1:203108/149‑1 (MQ=255)
                                      tCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGt                                                                                                             >  1:182941/1‑151 (MQ=255)
                                           ttCGTCAACTTCAATTTAAAATACATCAACCCCGGCGATTTTCTAAAACAAACCGCCTTGCCCTCCCACCCCCGGTTTGCCTGCCAGCGCGACAATGTTCCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTa                                                                                                           <  1:766575/148‑1 (MQ=255)
                                            tCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGCTACCTc                                                                                                       >  2:61988/1‑151 (MQ=255)
                                              cgtaaCTTCAAATGCAAAATACATCACCCCCGGCGAATTTCTAAAACAAACCGCCCTTCCCTCCCACCCCGGTTTTGCCTGCCAGCGCCCAAATGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTc                                                                                                       <  2:251908/147‑1 (MQ=255)
                                                cAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                                   >  1:329515/1‑151 (MQ=255)
                                                 aaCTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCg                                                                                                  >  1:527912/1‑151 (MQ=255)
                                                     taaaaTGTCAAATACACCAACCCGGGCAATTTCCTAAACCAGACCGCCCTTCCTTCCCACCCCCGGTTTGCCTCCCGGCGGCCCAGTGTTCCCACCCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCgg                                                                                                 <  2:263828/146‑1 (MQ=255)
                                                        actctaaaaTACAACAATCCGGGGGATTTTCTAAAACAGAGCCCCCTTGCCTCCCATCACTGTTGTGCCGACGAGGGCGCTAATGTTGCCTACCCCCATCTCCGCTGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                              <  2:106143/143‑1 (MQ=255)
                                                           gTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg                                                                                         >  1:471370/1‑151 (MQ=255)
                                                             cAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                                        >  1:610310/1‑150 (MQ=255)
                                                              aaaTACATCAACCCGGGCAATTTCCTAAAACAGAACCCCCTTCCCTCCCACCCCGGGTTGGCCTGCCGGCGCCCCAATTTTCCCAACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                          <  2:278147/147‑1 (MQ=255)
                                                               aaTACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAgg                                                                                     >  1:550240/1‑151 (MQ=255)
                                                                       acccccGGCAATTTTCTAAACCAGAACGCCCTTCCCTCCCACCCCCGTTTTGCCGCCCAGCGCGCCAATGTTGCCTACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                                 <  2:53558/144‑1 (MQ=255)
                                                                       aaccccGGCAAATTTCTTAAAAAAACCCCCTTTGCCTCCCACCCCCGTTTTGCCTGCCAGCGCGCAAATGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                                <  1:397668/145‑1 (MQ=255)
                                                                        aTCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                            >  1:70957/1‑151 (MQ=255)
                                                                         ccccGGCAAATTTCTAAAACAGACCCCCTTTCCTTCCCATCCCCGTTTTCCCTCCCAGCGCGCAATTTTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACgg                                                                           <  2:89602/150‑1 (MQ=255)
                                                                          cccGGCAAATTTCTAAAACAAACCGCCCTTGCCTTCCACCCCCGGTTGGCCTGCCAGCGCCCAAATGTTGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg                                                                          <  1:586182/151‑1 (MQ=255)
                                                                           ccGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGt                                                                         >  2:136254/1‑151 (MQ=255)
                                                                               ccaattttATTAAACATACCCCCCTCTCCTTCCATCACCTTTTTCCCGCCCAGCGGGCAATCGTAGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTc                                                                        <  2:352239/144‑1 (MQ=255)
                                                                                caaTTTTTTAAACAAAACCCCCCTTCCCTTCCACCCCGGTTTTCCCGCCCGGCCCGCAAATGTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt                                                                      <  1:535790/149‑1 (MQ=255)
                                                                                  attttCTTAAACAGACCGCCCTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCGCCATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                    <  2:597988/148‑1 (MQ=255)
                                                                                     tttttaaaccaaaaCGCCGTGGCCTCCCACCCCCGGTTGCCTTCCGGGGGCGCCAATTTTGCCAACCCCCAGCTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                                   <  1:649779/137‑1 (MQ=255)
                                                                                     ttCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                                 >  1:831417/1‑150 (MQ=255)
                                                                                          aaaCAGACCGCCTTTGCCTCCCACCCCCGTTTTGCCTCCCAGCGCGCAAATTTTGCCACCCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                           <  2:629442/151‑1 (MQ=255)
                                                                                            aCAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                           >  1:377796/1‑149 (MQ=255)
                                                                                               aGACCGCCCTTCCCTTCCATCACCGTTTTGCCTGCCAGCGCGCCAATGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGa                                                        <  2:344352/150‑1 (MQ=255)
                                                                                                 aacccccctgcccTTCCATCCCCGTTTTGCCTCCCAGCGCGCAAATGTTCCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                                      <  1:531329/140‑1 (MQ=255)
                                                                                                  aCGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                                      >  1:760732/1‑149 (MQ=255)
                                                                                                  aCGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                                      >  2:891161/1‑149 (MQ=255)
                                                                                                     gCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                                   >  2:501072/1‑150 (MQ=255)
                                                                                                      cccTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCCCCAATGTTGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                                  <  1:359629/150‑1 (MQ=255)
                                                                                                           gCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGc                                             >  2:324371/1‑150 (MQ=255)
                                                                                                            ccTCCCACCCCCGGTTGCCCCCCCAGCGCGCAAATGTTCCATACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGc                                             <  2:597208/149‑1 (MQ=255)
                                                                                                                  cATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCt                                      >  1:478516/1‑151 (MQ=255)
                                                                                                                       ccGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg                                  >  1:717913/1‑150 (MQ=255)
                                                                                                                        cGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg                                  >  2:914997/1‑149 (MQ=255)
                                                                                                                         ggTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                               >  2:696202/1‑151 (MQ=255)
                                                                                                                          gTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                               >  2:70602/1‑150 (MQ=255)
                                                                                                                            ttGCCTCCCAGCCCGCCATTGTTCCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAg                             <  2:283862/150‑1 (MQ=255)
                                                                                                                                     ccggcccgcaattGTTCCCAACCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcg                      <  2:410431/138‑1 (MQ=255)
                                                                                                                                        gcgcgcCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc                    >  2:431551/1‑149 (MQ=255)
                                                                                                                                               attGTTGCCTACCCCCACCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCgcgc           <  2:921807/149‑1 (MQ=255)
                                                                                                                                                       cTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGc   >  2:649723/1‑151 (MQ=255)
                                                                                                                                                         aaCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGcc  >  1:365470/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                    
GTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2577996‑2578291

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