Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC45.3% 82.2 / 11.7 54T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (21/8);  new base C (0/24);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 1.45e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376356‑1376648
                                                                                                                                                           |                                                                                                                                               
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                                >  2:849631/1‑151 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                            <  2:590818/150‑1 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCg                                                                                                                                           <  2:981509/151‑1 (MQ=255)
      cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                          >  2:772580/1‑151 (MQ=255)
       gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                         <  1:1078152/151‑1 (MQ=255)
          gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACCCCCTCGCCGAt                                                                                                                                      <  2:1031219/151‑1 (MQ=255)
          gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                      <  2:525769/151‑1 (MQ=255)
           aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                     <  1:681460/151‑1 (MQ=255)
             gACCTCGCCGCGTAGCCTGGACGGCGACCGTAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGc                                                                                                                                    <  2:929339/150‑1 (MQ=255)
               ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccg                                                                                                                                 <  2:980912/151‑1 (MQ=255)
                      gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                          >  2:437992/1‑151 (MQ=255)
                            ccgggACGGCGCCCGAAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTGTATCTCGACCCTTTAACCCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACAtc                                                                                                                       <  1:540788/145‑1 (MQ=255)
                                   ggCGACCGAAAAATGGTGAGTGATCCCCGGCGCGTACTGGAACAGGTCCCCGCGCGTACCGAACCCGTGTTGTATCTCGACCCTTTAACCCGCTATCAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGaaaa                                                                                                               <  1:572991/149‑1 (MQ=255)
                                       aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                         >  1:21432/1‑151 (MQ=255)
                                         cGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGGGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                         >  1:573522/1‑149 (MQ=255)
                                          gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                       >  1:975607/1‑150 (MQ=255)
                                                 ggTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                >  1:1076250/1‑150 (MQ=255)
                                                  gtgggtgaTCCTCGGCGCGTATTGGAACAGTTCCCCGCGGGTCCGGAACCGGGGGGTTATCGCGACCCGTAAACCCGCTATCGGAATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                 <  2:201612/148‑1 (MQ=255)
                                                        gaTTCCCTGGGCTTACGGGAACAGTTCGCCGCGCGTCCCGGACCGGGGGTGTTTCTCGACCCTTAAACCCGCTATCAGAATCATATAATAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                           <  1:763405/148‑1 (MQ=255)
                                                        gaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                        <  2:1009534/151‑1 (MQ=255)
                                                             ccTGCGCGTACTGGAACGGGTCGCGGCGCGACCCGGACCCGTGGTGTTTCTGAACCCTTAAAACCGTTATCAGGATCAAATGACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                   <  1:57152/151‑1 (MQ=255)
                                                              cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  >  1:1047986/1‑151 (MQ=255)
                                                                  gcgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCACGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                              >  2:979371/1‑151 (MQ=255)
                                                                       atGGAACCAGGTCGCCGCGCGACCCGAACCTGTGGTTTCTCCGGACCGTTAAACCCGCTATCAGGTTTATATAACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                            <  2:33269/148‑1 (MQ=255)
                                                                        tGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         >  2:441219/1‑150 (MQ=255)
                                                                          gaaagaggcgcccgcGCGTGCCGGAACCGTGGTGTATCCCAAACCTTAAACCCGTTATCAGGATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                        <  2:530735/139‑1 (MQ=255)
                                                                           aaCGGGTCGCCGCGCGTCCCGGACCGGTGGTGTATCTCAACCCTTAAACCCGCTATAAGAATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  2:509027/147‑1 (MQ=255)
                                                                           aaCAGGCCGCCGCGCGTCCCGGAACGGGGGTTTATCTCGAACCTTTAACCCGCTATCAGGATCATATAATAACCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                     <  1:941725/151‑1 (MQ=255)
                                                                            aacgggcccccgcGCTTACCGGACCCGTGGTGTATCCAAAACCTTAAACCCCCTATAAGGATCATATGATAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                      <  2:687108/141‑1 (MQ=255)
                                                                                ggTCGCCGCGCTTCCCGAAACGGTGTTGTCTCCCAACCCTTTAACCCGCTATCAGTATAATATAACCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                   <  2:43370/149‑1 (MQ=255)
                                                                                ggTCGCCGCGCGTCCCGGACCGGTGGTGTCCCTCAACCCTTAAACCCGCTATCAGAATCATATAAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                 <  2:335120/151‑1 (MQ=255)
                                                                                   cgccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              >  1:190494/1‑151 (MQ=255)
                                                                                    gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAACCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                             >  1:912177/1‑151 (MQ=255)
                                                                                     ccggGCGTACAGAAACGGTGGTTTACCTCCACCAGTTAACCCGCTATCAGGACCATAAGAAAACCCCCCACGCCGTTCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                               <  1:325434/148‑1 (MQ=255)
                                                                                     ccgcGCTTCCCGGACCGGTGGTGTATCTCGACCCTTTAAACCGCTATCAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                             <  1:26817/150‑1 (MQ=255)
                                                                                       gcgcgTCCCGGACCGGTGGTTTATCTCAACCCGTAAACCCCCTATCAGAATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                          <  1:429437/151‑1 (MQ=255)
                                                                                            tACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  2:285324/1‑147 (MQ=255)
                                                                                                   accgggggtgTATCCCAACCCTTAAACCGCCAACCAGAATAATAGAACCAACTCCCCCGCGGTTTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                <  2:735718/144‑1 (MQ=255)
                                                                                                    ccGTGGTATATCTCGAACCGTTAACCCCCTATCGGGATAATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                             <  1:359549/151‑1 (MQ=255)
                                                                                                    ccGGGGTGTATCTCGACCCTTTAACCCGCTATAAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                              <  1:891280/150‑1 (MQ=255)
                                                                                                      cggtgttgtcccccGACCCGTAAACCCGCTATCGGGATAAAATAACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                               <  2:84654/136‑1 (MQ=255)
                                                                                                        tgttgtATCTCGACCTTTTAACCCGCTATCAGAATCATATCATCACCCCCCCCGCTAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                              <  2:766061/144‑1 (MQ=255)
                                                                                                         ggttTATCTCGACCCGTTAACCCGCTATCAGGATAATATGATCACCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                         <  1:306344/151‑1 (MQ=255)
                                                                                                         ggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           >  1:729156/1‑149 (MQ=255)
                                                                                                              tATCTCGACCCGTAAACCCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                      <  1:1003955/150‑1 (MQ=255)
                                                                                                              tATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                                     >  2:1074025/1‑151 (MQ=255)
                                                                                                                     aaCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                               >  2:456177/1‑150 (MQ=255)
                                                                                                                      aCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                             >  1:225441/1‑151 (MQ=255)
                                                                                                                              aaCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                       >  2:743446/1‑150 (MQ=255)
                                                                                                                                 ccGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                    >  1:369539/1‑151 (MQ=255)
                                                                                                                                      tATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                >  1:56485/1‑151 (MQ=255)
                                                                                                                                          aGGATCATATAATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg            <  2:768988/151‑1 (MQ=255)
                                                                                                                                            gATAATATGACCAACCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACGTGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg             <  2:345408/148‑1 (MQ=255)
                                                                                                                                                   tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc    >  1:339427/1‑150 (MQ=255)
                                                                                                                                                     aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  1:343829/1‑150 (MQ=255)
                                                                                                                                                           |                                                                                                                                               
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376356‑1376648

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: