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

TTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTG  >  NC_000913/2577997‑2578289
                                                                                                                                                         |                                                                                                                                                   
ttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAAcccc                                                                                                                                                >  2:378765/1‑151 (MQ=255)
  cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCa                                                                                                                                               <  1:691542/150‑1 (MQ=255)
   gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGc                                                                                                                                             <  2:721852/151‑1 (MQ=255)
   gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGc                                                                                                                                             <  2:297466/151‑1 (MQ=255)
        cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAccgc                                                                                                                                        >  1:903849/1‑151 (MQ=255)
        cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAccgc                                                                                                                                        <  2:859124/151‑1 (MQ=255)
          gcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCATGTTCAAATACATCAAGCCCGGCGAATTTCTTAAACAGACCGCCCTTCCCTTCCATCACGGCTTGGCCTGCCAGGGCGCCAATGTAGCCTAACCCCTGCCccgccg                                                                                                                                      <  1:268987/151‑1 (MQ=255)
           cgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGt                                                                                                                                     >  2:58620/1‑151 (MQ=255)
            gACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGt                                                                                                                                     >  2:709781/1‑150 (MQ=255)
                     aaTAAATTTGCCGGGGTCGAGTCCGTCAACTTAAATGTAAAATACACCAACCCCGGCAATTTTCTAAAACAGACCCCCCTTCCCTCCCACCCCCGTTTTGCCTCCCAGCGCGCAAATGTGGCCTACCCCCAGCACCGCCGTACCGTTa                                                                                                                              <  2:510795/148‑1 (MQ=255)
                      aTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTag                                                                                                                             >  1:878895/1‑148 (MQ=255)
                        aaaTTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTag                                                                                                                             >  2:21830/1‑146 (MQ=255)
                              gTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGcc                                                                                                                  >  1:92609/1‑151 (MQ=255)
                                  gggTGGTGTTCTCCAACTTGAATGCAAAATACGCAAATCCGGGAGAATTTTTAAAACAAAACCCCTTTGCTTCCCACCCCTGGTTTGCCTCCCGCCGCCCAACGTTGCCCAACCCCCGGCTCCGCCGTACCGTTAGAGATCACCGccacc                                                                                                               <  1:468924/150‑1 (MQ=255)
                                     tCGAGTTCGTCACCTTAAATTTAAAATCCAACACCCCGGGCGAATTTCTAAAACAGACCGCCCTTCCCTCCCACCACGGTTTTGCCTGCCAGCGCCCAAATTTTCCCTACCCCAAGCACCGCCGTACCGTTAGAGATCACCGccacca                                                                                                              <  1:1000105/148‑1 (MQ=255)
                                       gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                          >  1:206691/1‑150 (MQ=255)
                                         tttCGTCAATTTCAATGTAAAATACACCCCTCCCGGCGAATTTCTAAAACAACCCCCCCTTCCTTTCCACCCCCGTTTTGCCGGCCCGCGCGCATATGTTGCCTACCCCCGGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTa                                                                                                         <  1:134065/148‑1 (MQ=255)
                                                accTTCAATGCAAAAAACATCAACCCCGGCGAATTTCTAAAACAGACCGCCCTTCCCTTCCACCCCCGTTTTGCCTGCCGCCGCGCCAATTTTGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCggg                                                                                                  <  2:196247/147‑1 (MQ=255)
                                                  cTTAAATGTCAAATACATCACCCCGGGCAAATTTTTAAACAAAAACCCCCTTCCCTCCCACCCCGGTTTGGCCGCCCGGCCCCCCATTTTTCCAAACCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCg                                                                                                <  2:938859/149‑1 (MQ=255)
                                                   ttCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                              >  2:512162/1‑150 (MQ=255)
                                                      aaTGTCAAATACACCAACCCCGGCGAATTTCTAAAACAGACCGCCCTTCCCTTCCACCCCCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTACCCCAAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                            <  2:973122/149‑1 (MQ=255)
                                                        tGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                        >  2:917692/1‑151 (MQ=255)
                                                          tCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                                      >  2:555838/1‑151 (MQ=255)
                                                            aaaTACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                                      >  2:149942/1‑149 (MQ=255)
                                                             aaatacaccaccccGGGCAATTTCCTAAAAAAAACCCCCCTTCCCTCCCACCCCCGTTTTGCCTCCCGGCGCCCATATTTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                        <  1:530265/136‑1 (MQ=255)
                                                                       ccccGGCAATTTTCTAAAACAGACCCCCTTCCCCTCCCACCCCCGTTTGCCCCGCCCGCGCCCCATGTTCCCTACCCCCCCGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                          <  1:605280/149‑1 (MQ=255)
                                                                        cccGGAGATTTTTTAAAACAAACCCCCCTTCCCTCCCACCCCGGGTTGGCCGCCCAGACCCACAATGTTGCATACCCCAGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc                                                                           <  1:379411/148‑1 (MQ=255)
                                                                         ccGGGGAATTTCTTAAAAAAGACCCCCTTTCCCCCCCATCCCCGTTTTCCCGCCCGCCGCCCAAATTTGCCCTCCCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                               <  1:645972/144‑1 (MQ=255)
                                                                         ccGGCTAATTTATAAAACAAAACGCCTTTCCCTCCATCCCCCGTTTGCCCTCCCAGCGCGCAAATTTTGCCTCCCCCAAGACCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                              <  2:630061/144‑1 (MQ=255)
                                                                              gAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                  >  2:1012915/1‑151 (MQ=255)
                                                                                 tttCTAAAACAGACCCCCTTTCCCTCCCCTCCCCGGTTTCCCTCCCAGCGCCCAAATTTGCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                <  1:123899/150‑1 (MQ=255)
                                                                                 tttCTAAAAAAGACCCCCTTTCCCTTCCTCCCCGGTTTTCCCCCCCAGCCCCCCAATGTTGCCTACCCCCGCCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                <  1:108728/150‑1 (MQ=255)
                                                                                  ttCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                               >  1:541440/1‑150 (MQ=255)
                                                                                  ttCTAAAACAGACCCCCCTTGCCTTCAACCCCGGTTTTGCCTGCAAGCGCGCCATTGTTCAAAACCCCCAGCTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                               <  2:21537/150‑1 (MQ=255)
                                                                                   tCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                             >  1:435385/1‑151 (MQ=255)
                                                                                       aaaaCAGAACCCCCTTGCTTCCCACCCCGGGTTTGCCTGCCGCCCCCCTATTGTGGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                          <  1:784297/150‑1 (MQ=255)
                                                                                          aCAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                         >  2:480820/1‑149 (MQ=255)
                                                                                          aCAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                         >  1:25728/1‑149 (MQ=255)
                                                                                             aGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                     >  2:142614/1‑151 (MQ=255)
                                                                                                  cGCCCTTCCCTTCCACCCCGGGTTGGCCCGCAAGCCCGCAAATTTTCCCAAACCCAAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                                 <  2:891095/151‑1 (MQ=255)
                                                                                                   gCCCTTGCCTTCCGCAACAGCATTATCTGCTATCACGCCAATCTGGCGGATCCCCAGCACCGCCGTACCGTGAGAGACCTCCGCCACCAGGTTACCTCGGGCCGTATATCTGAAGGCTTTTAACGGGTCTTTTTCGACTTCAAGACAAgg                                                 >  1:492966/1‑150 (MQ=255)
                                                                                                   gCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                                >  2:144350/1‑151 (MQ=255)
                                                                                                    cctttcccTCCCACCCCCGTTTTGCCTGCCGGCGCGCAATTTTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGa                                                      <  1:336034/138‑1 (MQ=255)
                                                                                                       ttcccTCCCACCCCCGTTTTGCCTGCCGGCGCCCCATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGc                                              <  1:225129/146‑1 (MQ=255)
                                                                                                         tccTTCCCACCCCGGTTTTCCCGCCCGGCGCCCAAATGTCCCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                            <  2:981811/149‑1 (MQ=255)
                                                                                                              ccccccccGGTTGCCCTCCCGCCGCGCAATTTTTGCCAACCCCCCGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                        <  1:822777/142‑1 (MQ=255)
                                                                                                                  ctccGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACACCTgg                                  >  1:164755/3‑150 (MQ=255)
                                                                                                                   ccccgTTTGCCCTCCCGGCCCGCCAATGTCGCATACCCCCATCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACg                                       <  1:189997/141‑1 (MQ=255)
                                                                                                                      ggTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAg                              >  1:753034/1‑150 (MQ=255)
                                                                                                                      ggTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                             >  2:541738/1‑151 (MQ=255)
                                                                                                                                            cAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg          >  2:992638/1‑151 (MQ=255)
                                                                                                                                             attGTTGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAg                <  2:564744/142‑1 (MQ=255)
                                                                                                                                             attGTTCCCAACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcg            <  2:656555/146‑1 (MQ=255)
                                                                                                                                                 ttGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCt        <  2:773341/148‑1 (MQ=255)
                                                                                                                                                    ccTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTg  <  1:199827/151‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
TTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTG  >  NC_000913/2577997‑2578289

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