Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,921,1770AC41.8% 93.6 / 12.3 55G308G (GGT→GGGgudPputative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base A (26/6);  new base C (0/23);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 4.88e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCG  >  NC_000913/2921033‑2921324
                                                                                                                                                |                                                                                                                                                   
cATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCTCCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGc                                                                                                                                               <  1:739199/151‑1 (MQ=255)
 aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGc                                                                                                                                               >  2:693491/1‑150 (MQ=255)
     aGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAcc                                                                                                                                          >  2:186436/1‑151 (MQ=255)
        cGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgcc                                                                                                                                       >  1:782160/1‑151 (MQ=255)
         gatgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgccg                                                                                                                                      <  1:750190/151‑1 (MQ=255)
          atgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGa                                                                                                                                     <  2:688203/151‑1 (MQ=255)
              tcatcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATaa                                                                                                                                  >  2:1342752/1‑150 (MQ=255)
                  ccACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATaaaa                                                                                                                                <  1:448212/148‑1 (MQ=255)
                  ccACTCAACGTTGACGTAGTGGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCCCCCAGCACACCGCCGATAaaacc                                                                                                                              <  1:665954/150‑1 (MQ=255)
                   cACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgc                                                                                                                            >  2:998178/1‑151 (MQ=255)
                    aCTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgc                                                                                                                            >  2:1186672/1‑150 (MQ=255)
                    aCTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgc                                                                                                                            >  1:788379/1‑150 (MQ=255)
                        aaCGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaaacc                                                                                                                       >  1:820267/1‑151 (MQ=255)
                            ttGACGTAGTTGCAGAACCCCATCCCCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCACCCAGCACACCGCCGATAaaaccgcaaaccgc                                                                                                                     <  1:419635/149‑1 (MQ=255)
                             tGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcgg                                                                                                                  >  1:332450/1‑151 (MQ=255)
                                 gTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                              >  1:1285183/1‑151 (MQ=255)
                                 gTAGTTGCAGAACACCACCCCCACCGACAGAAACATGCCCATTACGATCGGTTTTTTACGCGAAATGTTCAGCGACCCCGGGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaa                                                                                                                <  2:473134/149‑1 (MQ=255)
                                 gTAGTTGAAGAACCCCACCCCCATCGACAGCAACATGCCCATTACGATGGGTGTTTTACGCGCAATGTCCAGCGACCCCGTGCGGCGCATCAGCCAGTCGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                              <  1:776559/151‑1 (MQ=255)
                                    gTTGCAGACCCCCACCCCCACCGACAGCAACATCCCCATTACGATCGGTGTTTTACGCGAAATGTCCAGCGACCCCGTGCGGCGCATCAGCCAGCCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacggaa                                                                                                           <  1:227670/151‑1 (MQ=255)
                                         aGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCca                                                                                                       >  2:1183818/1‑150 (MQ=255)
                                                     ccatcGACAGCAACATCCCCATTCCGACCGGTTTTTTACGCGAATTGTCCAGCGACCCCGTGCGGCGCATCAGCCAGCCGAAAATACCCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGccc                                                                                               <  1:810042/146‑1 (MQ=255)
                                                            cAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                                   >  1:787855/1‑151 (MQ=255)
                                                                acaatccccTTTCCGATGGGTGTTTTCCGCGCAATGTCCAGCGATCCCGGGCGGCCCACCAGCCAGCCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACa                                                                               <  2:866431/145‑1 (MQ=255)
                                                                aaCATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGACCCCGGGCGGCCCATCAGCCAGTGGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                                   <  1:1235999/147‑1 (MQ=255)
                                                                   aTGCCATTTCCGATGGGTTTTTTACGCGAAATGTCCAGCGACCCCGTGCGGCCCATCAGCCAGTCGAAAAAAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACAt                                                                              <  2:135330/149‑1 (MQ=255)
                                                                     gCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg                                                                          >  1:349627/1‑151 (MQ=255)
                                                                      cccATACCGATCGGTGTTTTACGCGCAATGTCCAGCGACCCCGTGCGGCCCATCAGCCAGCGGAAAAAAACCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                         <  2:1047466/151‑1 (MQ=255)
                                                                       ccATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                        >  2:1110524/1‑151 (MQ=255)
                                                                       ccATTACGATCGGTGTTTTACGCAAAATGTTCAGCGATCCCGGGCGGCGCATCAGCAAGTCGGAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                        <  1:214161/151‑1 (MQ=255)
                                                                        cATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                        >  1:439620/1‑150 (MQ=255)
                                                                            aCGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt                                                                   >  2:494334/1‑151 (MQ=255)
                                                                                tCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAc                                                               >  2:217529/1‑151 (MQ=255)
                                                                                      ttttACGCGAAATGTCCAGCGATCCCGTGCGGCGCATCACCCAGTCGAAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAg                                                            <  1:779960/148‑1 (MQ=255)
                                                                                       tttACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATa                                                         >  1:1012713/1‑150 (MQ=255)
                                                                                             cgcAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCggg                                                  >  2:875963/1‑151 (MQ=255)
                                                                                              acattcttCAGGCATCCCGCGCGGCGTCTCACCGAGTCGAAAATAACCCCACCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAcc                                                     <  2:228382/141‑1 (MQ=255)
                                                                                                attttccAGCGACCCCGGGCGGCGCACCAGCCGGTCGAAAAAAACCCCCCCAACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAc                                               <  1:857413/145‑1 (MQ=255)
                                                                                                agtGTCCAGCAACCCCGTGCGGCGCACCGGCCAGTGGAAAAAAACCCCCCCCAACACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa                                                 <  2:1328694/147‑1 (MQ=255)
                                                                                                aatttccAGCGACCCCGGGCGGCGCACCAGCCAGTCGAAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa                                                 <  2:1367440/143‑1 (MQ=255)
                                                                                                aaTGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAc                                               >  2:161542/1‑151 (MQ=255)
                                                                                                   gTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAcc                                              >  1:791232/1‑149 (MQ=255)
                                                                                                    ttCAGCGATCCGGGGCGGCGCACCACCCAGGCGAAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAg                                            <  1:251077/150‑1 (MQ=255)
                                                                                                         cGACCCCGGGCGGCGCATCAGCCAGCCGGAAATAACCCCCCCAACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGta                                         <  2:981828/148‑1 (MQ=255)
                                                                                                          gcccccGGGCGGCGCAGCAGCCGGTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAg                                            <  2:1285155/141‑1 (MQ=255)
                                                                                                           aTCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaataa                                     >  2:378835/1‑150 (MQ=255)
                                                                                                              ccGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaaga                                  >  2:452208/1‑150 (MQ=255)
                                                                                                                gTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagt                               >  1:1033047/1‑151 (MQ=255)
                                                                                                                    ggcgCATCAGCCAGTCGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaag                                <  1:1321847/146‑1 (MQ=255)
                                                                                                                      cgcATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTc                          >  2:1377409/1‑150 (MQ=255)
                                                                                                                              ccAGTGGGAAATAATCCCCCCCATCACACCGCCGAGAAAACCGCAAACCGCCGGAACGCAAGCCACAAAGCGAGCTGTCAGAAGCGACATCCCGCGTGCCTGCACCAGAGAAACCGGGAAACAGGTAATAAAGAAGGAAGGAAgggagggg                 <  1:809208/151‑6 (MQ=255)
                                                                                                                              ccAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg                    <  1:386354/148‑1 (MQ=255)
                                                                                                                              ccAGTCGAAAATAACCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAgg                       <  1:1385191/145‑1 (MQ=255)
                                                                                                                                    gaaaaTAATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAg           <  1:998199/149‑1 (MQ=255)
                                                                                                                                     aaaaTAACCCCCCCCAGCACGCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACa            <  1:699904/148‑1 (MQ=255)
                                                                                                                                       aataatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTAc        >  1:31914/1‑151 (MQ=255)
                                                                                                                                             cccacccaGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCg  >  2:105143/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCG  >  NC_000913/2921033‑2921324

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: