Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC44.6% 96.7 / 13.8 56T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (26/5);  new base C (0/25);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 2.73e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTC  >  NC_000913/1376356‑1376640
                                                                                                                                                        |                                                                                                                                       
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                        >  1:578932/1‑151 (MQ=255)
  ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAaca                                                                                                                                         <  2:877320/148‑1 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                    <  2:52254/151‑1 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                     <  1:757801/149‑1 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                     <  2:373392/149‑1 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                    <  2:617578/150‑1 (MQ=255)
       gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                 <  1:290081/151‑1 (MQ=255)
       gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                 >  2:222780/1‑151 (MQ=255)
          gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                              >  1:25584/1‑151 (MQ=255)
            tGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGc                                                                                                                            >  2:753430/1‑151 (MQ=255)
                         tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                               >  2:136427/1‑151 (MQ=255)
                            ccTGGACGGCCACCGCAAAAGGGTGGGTGATCCCCTGCGCGTACTGAACCAGGTCGCCGCGCGTCCCGGAACCGTGGTGTATCTCGACCCTTTAACCCGTTACCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                            <  1:101325/151‑1 (MQ=255)
                               ggccGGCGCCCAAAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCGGGGGTGTATCTCGAACCTTTAACCCGCTATCAGAATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                            <  2:632448/145‑1 (MQ=255)
                               ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                          >  1:844391/1‑150 (MQ=255)
                                  cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                      >  1:338955/1‑151 (MQ=255)
                                   ggCGACCGCAAAAGGGTGAGTGATCCCCTGCGCGTCCGGGAACGGGTCGCCGCGCGTACCGGACCCGTGGTTTATCTCGACCCGTTAACCCGCTATCAGGATCATATGATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                      <  1:981266/150‑1 (MQ=255)
                                   ggCGACCGAAAAAGGGTGGGTGTTTCCCTGCGCGTACTGGAACAGTTCGCCGCGCGTACCGGACCCGTGGTTTATCCCGACCCTTAAACCCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                        <  2:920407/148‑1 (MQ=255)
                                         ccaaaaaTGGGTGAGTGTTCCCCTGCGCGTACGGGAACGGGGCCCCGCGCGCCCGGGCCCCGTGGTATCTCTCGACCCGTAAACCCGCTATCAGGACAATAAGCAAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                    <  1:238782/144‑1 (MQ=255)
                                          gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                               >  1:997079/1‑150 (MQ=255)
                                             aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                            >  2:420417/1‑150 (MQ=255)
                                                  ggTGAGTGTTCCCCTGCGCTTCTGGAAACGGGTCGCGGCGCGCCCCGGACCGGTGTTTTATCCCGACCCTTTAACCCGTTACCAGGACCAAATAAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                         <  1:620905/149‑1 (MQ=255)
                                                    tggggtaTTCCCGCCGCTTATTGAAACACTTCGCGGCGCGTACCGAAACTGTGTTTACTCTCAACCTGTAAACCCGTTACACGAAACTTATATTCATCCCCCTCGCTGATTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                       <  1:770568/143‑1 (MQ=255)
                                                      agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                      >  2:195952/1‑148 (MQ=255)
                                                             cccTGCGCGTACGGGAACAGGTCGCCGCGCGTCCCGGAACCGTGGTGTATCCCGACCCGTCAACCCGTTACCAGGATCATATGAACCACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                              <  1:265986/149‑1 (MQ=255)
                                                                   gcgTACGGGAACGGGTCGCCGCGCGTACCGGACCCGTGGTGTATCTCGACCCGTAAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                          <  2:426722/147‑1 (MQ=255)
                                                                      gtatggaaaCGGGTGGCCGCGCTTACGGGAACCGGGGTGTATCTCGACCCGTTAACCCGCTATAAGGATCATAAGAACAACCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                      <  1:996130/142‑1 (MQ=255)
                                                                         tGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                 >  2:292202/1‑150 (MQ=255)
                                                                           gaaaCAGGTCGCCGCGCGTCCCGAACCGGTGTTTTATCTCGACCCGTTAACCCGCTATCAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                               <  2:206982/149‑1 (MQ=255)
                                                                             aaCGGGTCGCCGCGCTTCCCGGACCCGTGTTTTATCTCAAACCGTTAACCCGCTATCGGGATCATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                              <  2:902928/150‑1 (MQ=255)
                                                                             aaCAGGTCGCTGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                             >  2:622030/1‑151 (MQ=255)
                                                                               cAGGTCGCCGCGCGTACCGGACCCGTGGTTTACCTCGACCCGTAAAACCGCTATCAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                           <  2:6174/151‑1 (MQ=255)
                                                                                 ggTCGCCGCGCGTACCGGACCCGTGGTGTATCTCGAACCGTAAACCCGCTATCAGGATCATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                           <  2:560568/149‑1 (MQ=255)
                                                                                    cgccgcGCGTAGCGGAACCGTGTTGTATCTCGAACCGTTAACCGGCTATAGGGATCATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                        <  2:553733/149‑1 (MQ=255)
                                                                                     gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                     >  2:919648/1‑151 (MQ=255)
                                                                                      ccgcGCGTCCCGGACCCGTTGTGTATCTCAACCCGTAAACCCGTTATAAGGATCATATGATCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                    <  1:332681/151‑1 (MQ=255)
                                                                                        gcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                    >  2:690502/1‑149 (MQ=255)
                                                                                        gcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                    >  2:562185/1‑149 (MQ=255)
                                                                                         cgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                    >  2:809209/1‑148 (MQ=255)
                                                                                         cgcgTACCGAACCCGTGTTTTACCCCAACCCGTTAACCCGTTATCAGAATCATATAACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                  <  1:32725/150‑1 (MQ=255)
                                                                                          gcgTACCGGACCCGTGGTGTATCTCAACCCTTTAAACCGCTATCAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                 <  2:632803/150‑1 (MQ=255)
                                                                                            gTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                 >  2:907365/1‑148 (MQ=255)
                                                                                                 gaccccGTGGTGTTTCTCAACCCTTAAACCCGTTATAAGGATAATATGATCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                          <  1:2748/146‑1 (MQ=255)
                                                                                                   aaCCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                       >  1:910029/1‑151 (MQ=255)
                                                                                                      cGTGTTGTATCCCGACCCGTAAAACCGTTATCGGGATAATATAAAAACCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                           <  1:724624/144‑1 (MQ=255)
                                                                                                         gttgtATCTCGACCCTTAAACCCGCTATAAGAATAATATGACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                   <  1:703505/147‑1 (MQ=255)
                                                                                                             tATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                             >  1:694689/1‑151 (MQ=255)
                                                                                                             tATCCCGACCCGTAAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                               <  2:201037/149‑1 (MQ=255)
                                                                                                                    acccGTTAACCCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                            <  2:524283/143‑1 (MQ=255)
                                                                                                                     aCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                     >  2:663155/1‑151 (MQ=255)
                                                                                                                        gTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                  >  1:722288/1‑151 (MQ=255)
                                                                                                                         ttAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                  <  2:143847/150‑1 (MQ=255)
                                                                                                                               ccGTTATCAGGATCATATAAACACCCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGa             <  1:940540/150‑1 (MQ=255)
                                                                                                                                 gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt          >  2:1066/1‑151 (MQ=255)
                                                                                                                                  cTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg         >  2:621563/1‑151 (MQ=255)
                                                                                                                                  cTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg         >  1:594792/1‑151 (MQ=255)
                                                                                                                                    aTCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt       >  2:618310/1‑151 (MQ=255)
                                                                                                                                        ggATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt   >  1:250896/1‑151 (MQ=255)
                                                                                                                                         gATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc  >  1:255332/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                       
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTC  >  NC_000913/1376356‑1376640

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: