Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,6610TC51.7% 62.4 / 86.9 58I44V (ATC→GTC) dppBdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (22/6);  new base C (0/30);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 6.70e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-11

AACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTA  >  NC_000913/3705516‑3705805
                                                                                                                                                 |                                                                                                                                                
aaCTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCAc                                                                                                                                             >  2:532104/1‑151 (MQ=255)
 aCTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACg                                                                                                                                            >  1:297580/1‑151 (MQ=255)
   tctTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGtt                                                                                                                                          >  2:140662/1‑151 (MQ=255)
      tCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTc                                                                                                                                         <  1:444830/149‑1 (MQ=255)
              cGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCAAAATGTAATGGGGATACGGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCAGGTTcgcccgcc                                                                                                                                  <  2:365925/148‑1 (MQ=255)
                ggATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGat                                                                                                                             >  1:342838/1‑151 (MQ=255)
                ggATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGat                                                                                                                             >  2:129649/1‑151 (MQ=255)
                     cGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGTCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGatcatc                                                                                                                         >  2:933691/1‑150 (MQ=255)
                         ttttCATTAAAATGCCTAGCCCGCCATGCATACCCCCCAAAATGTAAGGGAGATCCTGCTCCCACATCGTTTTATCTAGGCCGAGTTCACCCCGCCGCTGCGCGTGCCGCTCGGGGGAGCCCCCACGTTCGCCCGCCATGATCATCAc                                                                                                                       <  2:158276/148‑1 (MQ=255)
                          ttttataaaaaTGCCAAGACCCCCATGCATACCGCCCAAAATGGAAGGGAGATACTGCTGCCACAGCGGTTTATCAAACCCGAGTTCACCCAGCAGCTGCGCTTGCCGCTCGGGGGAGCTCCCACGTTCGCCCGCCATGAtcatca                                                                                                                        <  2:730989/138‑1 (MQ=255)
                             cATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCg                                                                                                                >  1:34502/1‑151 (MQ=255)
                                tGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGccc                                                                                                             >  1:630716/1‑151 (MQ=255)
                                  aaaTGCCTAGATCGCCATGCATACCCCCCAAAATGTAATGGAGATACTGCTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGTGCCGCTGGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGccc                                                                                                             <  2:658974/149‑1 (MQ=255)
                                   aaTGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                          >  2:469487/1‑151 (MQ=255)
                                       ccTAGTTCGACATCAAAACCGCCCAAAAAGTAAGGGCGATAGTGCTCCCACTTGGTTTTACCGACGCCGTGCCCAACGACCCGCTGGGCGGGACGCAGGGGGGTGTCCCCGTGTTCGTCCGCCATGATCGTCACCGGATCGCCCGGGATc                                                                                                       <  2:548639/150‑1 (MQ=255)
                                       ccTAGACCGCCAGCCATACCCCCCAAAATGTAAGGGAGATACTGCTCCCACAGCGTTTTATCTAAGCCGATTCCAGCAAGCAGCGGCGCGTGCCGCTGTGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGAt                                                                                                        <  2:109740/149‑1 (MQ=255)
                                          aGATCCCCATGAATACCCCCCAAAATGAAAGGGAGATCCTGCGCCCACATCGTTTTATCAAAGCCGATTTCAGCCAGCGGCTGCGCGTGCCGCTCGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGAt                                                                                                        <  1:446838/146‑1 (MQ=255)
                                          aGACCGCCATGCAAACCCCCCAAAATGTAAGGGAGATCCTGCCGCCACATCGTTTTATCAAAGCCGAGTTCACCCAGCAGCTGCGCGTGCCGCTGGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgt                                                                                                   <  2:15531/151‑1 (MQ=255)
                                              cGCCAGGCATACCCCCCAAAATGTAAGGGAGATCCTGCTCCCACATCGTTTTATCTAAGCCGATTGCAGCCCCCAGCTGCGCGGGCCGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgtg                                                                                                  <  2:305256/148‑1 (MQ=255)
                                                ccATGCATAACGCCCAAAATGTATGGGAGATCCTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGGGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaa                                                                                             <  1:305167/151‑1 (MQ=255)
                                                  aTGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaa                                                                                             >  2:521469/1‑149 (MQ=255)
                                                  aTGCATAACCCCCAAAATGTAGGGGAGATCCTGCTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCCGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAg                                                                                           <  2:510243/151‑1 (MQ=255)
                                                  aTGAATACCGCCCAAAATGTAAGGGGGATCCTGCGCCCACCTGGTTTTATCAAAGCCGAGTCCGGCCAGCAGCTGGGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaaa                                                                                            <  1:944739/150‑1 (MQ=255)
                                                    gCATAACCCCCAAAATGTAATGGAGATACTGCTGCCACATGGTTTAATCTAAGCCAAGTCCACCCGGCGGCGGCCCGTGCCGCTCTGGGGAGCCCCCACGTCGGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaa                                                                                             <  2:548975/147‑1 (MQ=255)
                                                        accccccaAAATTTAAGGGAGATCCTGCTGCCACATGGTTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaaa                                                                                      <  1:647264/146‑1 (MQ=255)
                                                        aaCGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAt                                                                                     >  1:35200/1‑151 (MQ=255)
                                                          ccccccAATGGTAGTGGAGATATTGCCGCCACGGCTTTTTATAAACGCGTATTCCGCCCGGCAGCGGCGAGTGCCGCTTTGGGGATTTCCAGCGTGGGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaaa                                                                                      <  1:467516/146‑1 (MQ=255)
                                                          cGCCCCAAATGTAATGGAGATCCTGCTGCCACAGCGTTTTATCTAAGCCGAGTCCAGCCAGCAGCTGCGGGTGACGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtg                                                                                    <  1:259750/150‑1 (MQ=255)
                                                               cAAATGTAAGGGAGATCCGGCGCCCACTTGTTTTTATCAAAGCCGAGTTCAGCCCCCAGCTGGGCGGGACGCTGGGGGGAGATCCCAGGTGCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                                               <  2:65768/150‑1 (MQ=255)
                                                                aaaTGTAAGGGAGATCCTGCGGCCACATCGGTTTATCTAAGCCGAGTTCCGCCAGCCGCTGCGCGTGACGCTTGGGGGAGACCCCGCGTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAg                                                                             <  2:524322/151‑1 (MQ=255)
                                                                      aagggAGATCCTGCTCCCACTTCGGTTTATCAAAGCCGATTTCAGCCAGCACCGGCGCGTGCCGCTTGGGGGAGACCCCACGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                                          <  2:766858/145‑1 (MQ=255)
                                                                      aaTGGAGATTCTGTTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                           <  1:774613/147‑1 (MQ=255)
                                                                             aTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGAt                                                                 >  2:79701/1‑150 (MQ=255)
                                                                                 tgcgcccACATCGTTTTTTCTAGGCCGATTCCAGCCAGCAGCGGCGCGTGCCGCTGTGGGGAGACCCCCCTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGa                                                                  <  1:684809/140‑1 (MQ=255)
                                                                                    tgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgg                                                         >  2:528806/1‑151 (MQ=255)
                                                                                    tgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGCAATACCGATCAACGTgg                                                         >  2:513560/1‑151 (MQ=255)
                                                                                      ccacagggttttATCAAAGCCGAGTTCGCCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg                                                          <  2:654016/139‑1 (MQ=255)
                                                                                      ccACATCGTTTTATCTAAGCCGAGTCCAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGt                                                           <  1:25075/147‑1 (MQ=255)
                                                                                      ccACATCGGTTTATCTAGGCGGATTTCAGCCAGCAGCGGCGGGTGCCGCTCGGGGGAGATCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg                                                          <  1:663696/148‑1 (MQ=255)
                                                                                             gttttATCAAAGCCGATTCCAGCCAGCAGCTGCGCGTGCCGCTGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGAc                                                  <  1:260746/147‑1 (MQ=255)
                                                                                                     tAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTcc                                            >  2:343415/1‑147 (MQ=255)
                                                                                                      aggcggAGTTCACCCAGCAGCTGCGCCTGCCGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                       <  1:88728/146‑1 (MQ=255)
                                                                                                      aggCCGGGTTCAGCCGGCCGCGGCGCGGGCCGCCTTGGGGGCCCCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTcc                                            <  1:11053/144‑1 (MQ=255)
                                                                                                         ccGAGTCCAGCCAGCAGCTGCGCGGGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCa                                          <  2:672916/145‑1 (MQ=255)
                                                                                                            aGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                  >  2:934924/1‑150 (MQ=255)
                                                                                                                         gtggcgctggACGCTGTGGGGCGCCCCCAGGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTg                       <  1:848142/139‑1 (MQ=255)
                                                                                                                           tGCGCGTGACACTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAAc                   >  2:533303/1‑150 (MQ=255)
                                                                                                                              gcgTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATa                >  1:739441/1‑150 (MQ=255)
                                                                                                                                gggccgctggggggAGACCCCACTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAc              <  1:361501/140‑1 (MQ=255)
                                                                                                                                gTGCCGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCaa                    <  1:538118/144‑1 (MQ=255)
                                                                                                                                 tGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAccc            >  2:70839/1‑151 (MQ=255)
                                                                                                                                   aCGCTGTGGGGGGGCCCCCTGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGTGTCCCAAACGTCGGAGAATAAACTGCAACAt                 <  2:216329/144‑1 (MQ=255)
                                                                                                                                    cgctggggggAGCCCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAccc            <  2:375015/142‑1 (MQ=255)
                                                                                                                                    cGCTCTGGGGAGCCCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCgg          <  1:159683/150‑1 (MQ=255)
                                                                                                                                       tctGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtc      >  2:906155/1‑151 (MQ=255)
                                                                                                                                         tGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtctc    >  2:914680/1‑151 (MQ=255)
                                                                                                                                          ggggAGACCCCACGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAt        <  2:566514/146‑1 (MQ=255)
                                                                                                                                              agaTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGGGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCta  >  2:490482/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTA  >  NC_000913/3705516‑3705805

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