Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,1760TC46.6% 54.1 / 32.8 58E304G (GAA→GGA) ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (24/7);  new base C (0/27);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 2.48e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGT  >  NC_000913/3601028‑3601316
                                                                                                                                                     |                                                                                                                                            
ttattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTc                                                                                                                                             <  2:560218/150‑1 (MQ=255)
  attCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGc                                                                                                                                          <  2:429896/151‑1 (MQ=255)
  attCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGc                                                                                                                                          >  1:33579/1‑151 (MQ=255)
       ggCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAAcc                                                                                                                                      <  2:1036408/150‑1 (MQ=255)
       ggCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAAcc                                                                                                                                      >  1:584479/1‑150 (MQ=255)
              aGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATg                                                                                                                               >  1:343292/1‑150 (MQ=255)
                  gCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAAt                                                                                                                          >  1:840113/1‑151 (MQ=255)
                      aaaTGTTGTCCCGGGGCAACCCCCGCTCCCCCCCAGCCAATCATCGAGCATACCAGCAGAAATAGCAGCCATTCAACAAATGATAACCCATTAATATAAACCTTCGTCCCAAAAC‑CCTGTCCCCCTCCCGCACCCGCCGATGACAATCgc                                                                                                                       <  1:250380/150‑1 (MQ=255)
                       aaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCgc                                                                                                                       >  1:233939/1‑149 (MQ=255)
                           ttggccggTGGCAACCCCCCCGCCCCCCCAACCAATCATCGACCACCCCAGCAGCAATACCAGGCATTCATCGAATAAAACGCAATTAATAAAAACCTTCGTACCAAAAA‑CCGGTCCCCCTCCCGAAACCGCCGATGACAATCgcagca                                                                                                                   <  2:655928/142‑1 (MQ=255)
                              tACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCAcc                                                                                                                 >  2:645955/1‑148 (MQ=255)
                              tACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCAcc                                                                                                                 >  2:226167/1‑148 (MQ=255)
                                  gTGGAAACCCACGCGCCCCCCCACCAAATCACCGAGCATCCCAGAAGAAATAGCAGGCATTCATCAAATGAAAAGCCATTGATAAAAAACTTCGTCCCAAAAC‑CCTGTGCCACTCCCACAACCGCCGATGACAATCGCAGCACCAg                                                                                                               <  1:674552/146‑1 (MQ=255)
                                       aaGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGAc                                                                                                     >  1:1007179/1‑151 (MQ=255)
                                          ccACGCTCCCCCCCAGCAAATCACCGACCATCCCAGAAAAAAAAGCAGGCATTCATCAAATGAAAAGCCATTAATATAAACCTTCGTCCCAAAAC‑CCTGTCCCACTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaa                                                                                                   <  1:901563/150‑1 (MQ=255)
                                               cccccccccacccaATCACCGACCAACCCAGCAGCAATAGCAGGCATTAACCAAATGATAAGCAATTGATATCAACCTCCTTCCCAAACC‑CCTGTCCCCTTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAt                                                                                                <  1:487611/142‑1 (MQ=255)
                                               cTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTa                                                                                              >  1:68932/1‑150 (MQ=255)
                                                  ccccccAGCCAATCATCGACCATCCCAGCAACAAAGGCAGGCATTCACCGAATAAAAAGCCATTAAAATCAACTTTCTTCCCAAAAA‑CCTGTGCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAAt                                                                                            <  1:249740/146‑1 (MQ=255)
                                                     cccAGCAAATCAACGAGCATACCAGAAAGAAGAACGGGCATTCATAGAAGGATAAGCATTTGATACAAATCTCCGTCCCAAAAG‑CCGGGCCCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGaca                                                                                        <  1:613177/150‑1 (MQ=255)
                                                     cccAACAAATCACCGAGCATCCCAGAAGAAATAGCAGGCTTTCATCGAATGATAAGCATTTAATATAAAACTTCGTTCCAAAAC‑CCGGTGCCCTTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGac                                                                                         <  1:305910/149‑1 (MQ=255)
                                                     ccaaacaaaaaatcGAGCATACCAGCAACAAGCGCGGGCATTCATCAAATAAAAAGCATTTGATATAAAACTTCGTCCCAAAAC‑CCTGTCCCTCCTCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaa                                                                                       <  2:313526/140‑1 (MQ=255)
                                                            aaTCATCAACCATCCCAGCAGAAATAGCGGGCATTCATCAAATAATAAGCTATTAATAAAAATCTTCGTCCCAAAAC‑CCTGTGCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATg                                                                                  <  1:183656/149‑1 (MQ=255)
                                                               catcGAGCATACCAGCAGAAAAAGCAGGCATTCATAAAATAAAAAGCCATTAATAAAAACCTCCGTCCCAAAAC‑CCTGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgcg                                                                                <  2:879077/148‑1 (MQ=255)
                                                               caccGACCACACCAGCAACAATAGCAGGTTTTCATCAAATGATAAGCCTTTAATATAAATCTTCGTCCCAAAAC‑CCTCGCCCCCTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATg                                                                                  <  1:849206/143‑1 (MQ=255)
                                                                       atcccAGCAGAAATAGCAGGCATTCATCGAATAAAAAGCTATTGATATCAACCTTCGTCCCAAAAC‑CCTGTGCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgc                                                                                 <  1:45768/136‑1 (MQ=255)
                                                                          ccAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCTTTGATATAAACCTTCGTCCCAAAAC‑CCTGTCCCTCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATccc                                                                    <  1:588243/149‑1 (MQ=255)
                                                                            agcagaAATAGCAGGCTTTCATGAAATAAAAACCAATTTATATAAACCTTCTTCCCAAAAC‑CCTGTGCCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGaaaa                                                                        <  1:734121/143‑1 (MQ=255)
                                                                                  aaaagcagcCATTCAAAGAATAATAAGCTATTAATATCAACCTTC‑TTCCAAAAACCCTGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgc                                                            <  1:895661/146‑1 (MQ=255)
                                                                                  aaTAGCGGGCATTCATCAAATAATAAGCCATTAATATAAACCTTCGTCCCAAAAC‑CCTGTCCCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTg                                                             <  1:51117/148‑1 (MQ=255)
                                                                                    tAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGcc                                                         >  1:690981/1‑150 (MQ=255)
                                                                                      gCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGAAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGccacc                                                      >  2:633354/1‑151 (MQ=255)
                                                                                         ggaaTTCATCAAATAAAAACCCATTGAAATCAACCTTCTTTCAAAAAC‑CCTGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGcca                                                        <  2:829643/143‑1 (MQ=255)
                                                                                            aTTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATAc                                                  >  1:791923/1‑149 (MQ=255)
                                                                                             ttCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAgg                                               >  1:391741/1‑151 (MQ=255)
                                                                                                aTCAAATGAAAACCCTTTGATATAAACCTTCTTCCCAAAAC‑CCTGTCCCACTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGa                                              <  1:824871/149‑1 (MQ=255)
                                                                                                        aTAACCCATTTATATAAACCTTCGTCCCAAAAC‑CCTGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAg                                    <  1:500607/151‑1 (MQ=255)
                                                                                                          aaGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAgg                                   >  1:800312/1‑150 (MQ=255)
                                                                                                            gCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATg                                >  2:1066848/1‑151 (MQ=255)
                                                                                                             ctaTTAATATCAAACTTCGTCCCGAAAC‑CCGGCCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt                                 <  1:513101/147‑1 (MQ=255)
                                                                                                             ccATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATg                                >  1:909162/1‑150 (MQ=255)
                                                                                                                ttGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATc                            >  2:499557/1‑151 (MQ=255)
                                                                                                                 taatatAAAACTTC‑TTCCAAAAACCCTGTCCCCCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATc                            <  1:429904/148‑1 (MQ=255)
                                                                                                                 tGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATcc                           >  1:288160/1‑151 (MQ=255)
                                                                                                                  gATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATcc                           >  1:972970/1‑150 (MQ=255)
                                                                                                                   atatCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCa                          >  1:601232/1‑150 (MQ=255)
                                                                                                                    tatCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGA‑AATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATc                        >  2:315017/1‑150 (MQ=255)
                                                                                                                      acAACCTTCGTACCAAAAC‑CCTGCCCCCCTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCt                       <  1:858177/149‑1 (MQ=255)
                                                                                                                       cAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGt                     >  2:332493/1‑151 (MQ=255)
                                                                                                                        aaaCTGCGTCCCAAAAC‑CCTGTCCCCTTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTc                    <  2:825544/151‑1 (MQ=255)
                                                                                                                         acctccgtccaaaaac‑CCTGTCCCCCTCCCGAAGCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCt                       <  2:584347/135‑1 (MQ=255)
                                                                                                                         accTTCGTCCCAAAAC‑CCTGTGCCACTCCCCCAACCGCCGATGACAATCGCAGCACCTGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCg                   <  1:676544/149‑1 (MQ=255)
                                                                                                                                ttcaaaaac‑CCTGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCaa             <  1:998692/145‑1 (MQ=255)
                                                                                                                                ttCCAAAAC‑CCTGTGCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCa                 <  1:11515/146‑1 (MQ=255)
                                                                                                                                  ccGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATc           >  2:582461/1‑150 (MQ=255)
                                                                                                                                     aaag‑CCTGTGCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCAc                <  1:671560/142‑1 (MQ=255)
                                                                                                                                     aaac‑CCTGTCCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCa          <  1:468366/148‑1 (MQ=255)
                                                                                                                                      aaacctgaccccACTTCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCa          <  2:354680/138‑1 (MQ=255)
                                                                                                                                           cTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGt  >  1:43449/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
TTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGT  >  NC_000913/3601028‑3601316

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