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

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT  >  NC_000913/1376355‑1376650
                                                                                                                                                          |                                                                                                                                                 
cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAaca                                                                                                                                                   <  1:1059868/151‑1 (MQ=255)
   ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                                  <  2:393594/149‑1 (MQ=255)
    tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAcc                                                                                                                                                 <  2:60230/149‑1 (MQ=255)
     aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                <  2:574436/149‑1 (MQ=255)
       cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                            <  2:458886/151‑1 (MQ=255)
        gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                           >  2:71456/1‑151 (MQ=255)
         ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                          >  1:1009979/1‑151 (MQ=255)
           gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                        >  1:299850/1‑151 (MQ=255)
                ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccg                                                                                                                                   <  1:433921/151‑1 (MQ=255)
                       gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATAATATGATCAACACCCTCTCCGATGCGCGCCGTTAc                                                                                                                            >  2:539386/1‑151 (MQ=255)
                               tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                      >  2:270199/1‑149 (MQ=255)
                                ggACGGCGCCCGAAAAAGGGTGAGTGATCCCCGGCGCGTACGGAAACAGGCCGCCGCGCGTCCGGGACCGGTGTTTTATCCCAACCTTTTAACCCCCTATAAGGATCATATAAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                   <  1:246679/151‑1 (MQ=255)
                                    ggCGACCGAAAAGGGGGGGTTGTTCCCCGGCGCGTACGGAAACGGGTCCCCGCGCGTACCGAACCTGTGGTGTATCTCGACCCGTTAACCCTCTATCAGGTTCATTTGATCAACACCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                      <  2:61414/144‑1 (MQ=255)
                                    ggCGACCAAAAAGTGGTGAGTGTTCCCGTGCGCTTACTGAAACAGGTCGCCGCGCGTACCGGACCGGTGGTGTTTCTCAACCCTTAAACCCGCTACCAGGATCATATGACCAACCCCCTCGCAGATCCCCGCCGTTACATCGTCGAAAACg                                                                                                               <  2:494400/151‑1 (MQ=255)
                                            aaaaggggTGAGTGATTCCCTGCGCTTCCGGAAACGGGCCCCCGCGCTTCCGGAACCCGTGTTGTATCCCGACCCTTAAACCGCCTATCAGAATAATATAAAAACCCCCCGGGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                        <  1:295993/144‑1 (MQ=255)
                                             aaaaTGGTGAGTGCTCCCCGGCGCGTACGGAAACGGGTCGCCGCCCGTACCGGAACGGTGGTGTTCCTCGACCCTTAAACCCGCTACCAGGATAATATGACCACCCCCCTCGCCAAGGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                         <  1:739458/148‑1 (MQ=255)
                                             aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                      >  1:441991/1‑151 (MQ=255)
                                             aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                      >  2:478845/1‑151 (MQ=255)
                                             aaaaGGGTGAGTTATCCCCGGCGCTTCTTGAACCAGGTCCCGGCGCGTCCCGAAACGGGGGTGTATCTGAAACTTTAAACCCGTTATCGGTATCATAGGATCAAACCCCCCGCGGAGCCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                        <  2:316349/149‑1 (MQ=255)
                                                 tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                  >  2:168262/1‑151 (MQ=255)
                                                  ggggagtgaTTCCCGGCGGGTCCGGGAAAAGGTCGCGGCGCGTACCGGACCCGTGTTGTATCCCGAACCTTTAAACCGTTATCAGTTTAATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                                 <  1:187148/148‑1 (MQ=255)
                                                         gaTTCCCGGCGCGTACGGAAGCAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAACCCGTTATCAGGATCATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                           <  1:462221/150‑1 (MQ=255)
                                                          aTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                         >  1:807054/1‑151 (MQ=255)
                                                             cccTGCCCGTACGGGAACGGTCCGCCGGGCGTCCCGAACCCGGGGTGTACCTCAAACCTTAAACCCCTTATAAGGACAATAGGAACACCCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                         <  2:861185/148‑1 (MQ=255)
                                                                 gcgcTTACGGGAACGGGTCGCGGCGCGTCCGGGACCGGTGGTGTATCTCGACCCTTTAACCCGTTATCAGGTTCTTATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                       <  1:521170/146‑1 (MQ=255)
                                                                   gcgTACTGGAACAGGTCGCCGCGCGTCCCGGAACCGTGGTGTATCTCGACCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                <  1:303794/151‑1 (MQ=255)
                                                                          ggAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                         >  1:94164/1‑151 (MQ=255)
                                                                             aCAGTTCCCCGCGAGTCCCGGTGCTGTGGTCTATCCCGCACCATTAAGCCGCTAACGGCATCATTTTACCAGCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                       <  1:200742/150‑1 (MQ=255)
                                                                                ggTCGCCGCGCGTACCGGAACCGTGGTGTACCCCGAACCTTTAACCCGCTACCAGTATCATATGACCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                     <  2:883575/149‑1 (MQ=255)
                                                                                      cgcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                              >  1:179052/1‑150 (MQ=255)
                                                                                         gcgTCCCGGAACCGGGGTGTATCCCGACCCGTAAACCCGCTACCAGAATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                          <  1:630470/151‑1 (MQ=255)
                                                                                           gtcccGGACCGGGGGTGTTTCCCAAACCTTTAACCCGATATAAGGATCATATAATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                           <  2:1027272/145‑1 (MQ=255)
                                                                                             aCCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                       >  2:451707/1‑150 (MQ=255)
                                                                                               cGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                    >  1:867865/1‑151 (MQ=255)
                                                                                                 ggaccggtggtgTACCCCGACCCGTAAAACCGCTATAGGGTTCATATGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                         <  1:644204/139‑1 (MQ=255)
                                                                                                  aaCCGTGGTTTATCCCGACCCGTAAACCCGCTATCAGAATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                    <  1:515169/148‑1 (MQ=255)
                                                                                                           tttttcccaaaCCGTTAACCCCTTATCAGAATAATATAACCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                           <  2:247383/140‑1 (MQ=255)
                                                                                                           tgtATCTCGACCCTTTAACCCCCTAGCAGAATAATATAAACAACACCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                              <  1:510171/146‑1 (MQ=255)
                                                                                                           tgtACCCCGAACCTTTAACCCTTTATCAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                          <  1:295879/150‑1 (MQ=255)
                                                                                                               tctcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                     >  2:175361/1‑151 (MQ=255)
                                                                                                                 cccgaCCATTAAACCCGCTATCAGAATAATATAACCAACCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                <  2:616037/135‑1 (MQ=255)
                                                                                                                  caacccGTTAACCCGCAATCAGAATCATATGATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                     <  1:451718/144‑1 (MQ=255)
                                                                                                                  cGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                  >  2:692078/1‑151 (MQ=255)
                                                                                                                   gAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                 >  1:58826/1‑151 (MQ=255)
                                                                                                                      cctttAACCCGTTACCGGGATCATATGACCACCCCCCCCGCCGAGCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                               <  1:10284/147‑1 (MQ=255)
                                                                                                                      ccGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                               >  1:264112/1‑150 (MQ=255)
                                                                                                                      ccGTAAACCCGCTACCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                              <  1:432980/151‑1 (MQ=255)
                                                                                                                          tAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                          >  2:1050152/1‑151 (MQ=255)
                                                                                                                                 cccttaccAGGATCATATAACAACCCCCCTCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                          <  1:665333/140‑1 (MQ=255)
                                                                                                                                   cTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                   >  2:887756/1‑151 (MQ=255)
                                                                                                                                     tATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                  >  1:594812/1‑151 (MQ=255)
                                                                                                                                      aTCAGGATAATATAAACAACCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                    <  2:994172/148‑1 (MQ=255)
                                                                                                                                         aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg              >  2:905641/1‑151 (MQ=255)
                                                                                                                                         aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg              >  2:942239/1‑151 (MQ=255)
                                                                                                                                          ggATAAAATGACCAACCCCCTCGCGGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                    <  2:222139/144‑1 (MQ=255)
                                                                                                                                              cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg         >  1:975080/1‑151 (MQ=255)
                                                                                                                                                  tGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc      <  1:409078/150‑1 (MQ=255)
                                                                                                                                                     tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat  >  1:549800/1‑151 (MQ=255)
                                                                                                                                                          |                                                                                                                                                 
CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT  >  NC_000913/1376355‑1376650

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: