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

CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376357‑1376652
                                                                                                                                                        |                                                                                                                                                   
cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAcc                                                                                                                                                   >  2:188136/1‑151 (MQ=255)
     cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                              <  2:1227086/151‑1 (MQ=255)
            gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGc                                                                                                                                        <  1:952777/150‑1 (MQ=255)
               cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCACCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                    >  2:287403/1‑151 (MQ=255)
                        tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                                           >  2:399644/1‑151 (MQ=255)
                              ggAGGGCGCCCGAAAAGGGGGGATTGTTCCCCGGCGCGTACGGGACCAGGTCCCCGCGCGTCCCGGACCCGTGGTGTACCCCGAACCGTAAAACCGTTATCAGGTCAATATGACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGa                                                                                                                      <  2:964136/150‑1 (MQ=255)
                                 cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                  >  2:641444/1‑151 (MQ=255)
                                       ccGAAAAAGGGTGAGTGTTCCCCTGCGCTTACTGGAACAGGCCGCCCCGCGTCCGGAACCCGGGGTGTATCCCGACCCTTAAACCCGTTATCAGGATAATAGGATCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                              <  2:291339/149‑1 (MQ=255)
                                       ccGAAAAAGGGTGAGTGATTCCCTGCGCGTCCGGGAACGGGCCGCCGCGCGTCCCGGACGCGTGGTGTTTCTCGACCCGTTAACCCGCTACCAGGATAATATGATCACCCCCCTCCCCGAGGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                            <  2:1174174/151‑1 (MQ=255)
                                                       gaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                             >  1:403749/1‑150 (MQ=255)
                                                       gaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                             >  1:990632/1‑150 (MQ=255)
                                                        aTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                            >  2:1092658/1‑150 (MQ=255)
                                                         ttCCCTGCGCGTACTGGAACAGGTCGACGCGCGTACCGGAACCGTGGTGTATCTCGAAGCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                           >  2:263068/1‑150 (MQ=255)
                                                           cccTGCGCGTACTGGAACGGGTCGCCGCGCGTCCCGAACCGGTGGTGTATCCCGACCCTTAAAACCGCTACCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                               <  1:1099159/144‑1 (MQ=255)
                                                                 gcgTACTGGAACAGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCTCGACCCTTTAACCCGCTATCAGGATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                      <  1:1011685/147‑1 (MQ=255)
                                                                  cgtgctgaactaGGTCGCCGGGCGTACCGGAACGGTGGTGTACCTCTAACCCTTAAACGCCTTTAAGAATCATAAAATCAACACCCCCGCCGCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                 <  1:1069759/140‑1 (MQ=255)
                                                                    tACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                >  2:116533/1‑150 (MQ=255)
                                                                     acgggAAAAGGTCGCGGCGCGTACCGGAACCGTGGTGTATCTCAACCCGTTAACCCGTTAACAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                   <  1:1188211/143‑1 (MQ=255)
                                                                      cggaaacgggcccccgcGCGTCCGGGACCCGGGGTGTACCTCAACCCGTAAACCCGCTATCAGAATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                              <  2:705345/137‑1 (MQ=255)
                                                                      cTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  1:764668/1‑151 (MQ=255)
                                                                       tgaaaCGGGTCCCCGCGCTTCCCGGAACCGGGTTTTATCTCCACCCTTTAAACCGCTATCAGAATCATATGACCAACACCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                              <  1:1156175/146‑1 (MQ=255)
                                                                        gaaaCAGGTCGCCGCGCGTCCGGGACCCGTGGTGTATCTCAACCGGTAAACCCGCTACCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                  <  1:834833/142‑1 (MQ=255)
                                                                         taaCAGGTCCCGCGGCGTACCGGACCCGGGTTGTACCTCAACCCTTAAACCCGCTATAAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             <  2:78151/147‑1 (MQ=255)
                                                                           acgggTCGCGGCGGGTCCCGAACCGGGGGTTTCTCTCAAACCGTAAACCCGCTATAAGAATAATATAAAAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                  <  2:869157/138‑1 (MQ=255)
                                                                           aagagcccgcggggcgTCCGGGAACGGTGTGTTATCGCGACCCGATAAACCGCTATCAGGATCATAGTAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                           <  1:1199596/136‑1 (MQ=255)
                                                                           aCAGGTCGCCGCGCGTACCGAACCCGTGGTGTATCTCAACCCGTAAACCCGCTATCAGGATAATATAATCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAGAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                        <  2:285464/151‑1 (MQ=255)
                                                                              ggTCGCCGCGCGTCCCGGACCCGTGGTTTATCTCAACCCGTAAAACCGATATCAGAATCATATGACCAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                      <  1:216829/150‑1 (MQ=255)
                                                                                 cgccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                                  >  2:1143794/1‑151 (MQ=255)
                                                                                     gcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                >  1:83909/1‑149 (MQ=255)
                                                                                      cgcgTCCCGGACCGGGGGTTTATCTCGAACCGTAAACCCGCTATCAGGATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                              <  2:1136583/150‑1 (MQ=255)
                                                                                          gTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  1:470648/1‑148 (MQ=255)
                                                                                               ggACCCGGGGTGTATCTCGACCCGTAAACCCGCTATCAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                      <  1:1189998/150‑1 (MQ=255)
                                                                                                  accggtggttTATCTCAACCCATAAACCGGTTATAAAGATCTAATGAACAACCCCCGCGCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                      <  2:1125435/144‑1 (MQ=255)
                                                                                                  accgggggttTATCTCAAACCGTAAACCCGCTACCAGAATAATATGATCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                   <  1:160267/145‑1 (MQ=255)
                                                                                                  aCCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                  >  2:1136628/1‑151 (MQ=255)
                                                                                                    ccGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 >  1:725320/1‑151 (MQ=255)
                                                                                                         gtgtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                            >  2:1210588/1‑151 (MQ=255)
                                                                                                          tttttctcGACCCGTAAACCCGTTATAAGAATAATATCATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                            <  2:989447/146‑1 (MQ=255)
                                                                                                           gtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                           >  2:416117/1‑150 (MQ=255)
                                                                                                             aTCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                        >  1:994411/1‑151 (MQ=255)
                                                                                                             aTCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                        >  1:993384/1‑151 (MQ=255)
                                                                                                             aTCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                        >  2:91480/1‑151 (MQ=255)
                                                                                                              tctcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                       >  1:563648/1‑151 (MQ=255)
                                                                                                                   aaCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                   >  1:722561/1‑150 (MQ=255)
                                                                                                                     ccGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                 >  1:1108971/1‑150 (MQ=255)
                                                                                                                           aacccGCTATCAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                            <  1:545151/147‑1 (MQ=255)
                                                                                                                           aacccGCAATAAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACaa                               <  1:896185/144‑1 (MQ=255)
                                                                                                                               cGCTATCAGGATAAAATAAACCACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       <  2:764655/151‑1 (MQ=255)
                                                                                                                                        ggATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt               >  2:756388/1‑151 (MQ=255)
                                                                                                                                        ggATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                 <  1:643898/149‑1 (MQ=255)
                                                                                                                                        ggATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                  <  1:479515/148‑1 (MQ=255)
                                                                                                                                                      acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg  >  2:523414/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376357‑1376652

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