Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG21.2% 193.2 / 53.9 80intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base G (16/1);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 2.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATC  >  NC_000913/3705887‑3706078
                                                                                                 |                                                                                                   
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                   >  1:411160/1‑100 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTCTcc                                                                                                >  2:987141/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                             <  1:1010741/100‑1 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                             <  1:593242/100‑1 (MQ=37)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                            <  1:106762/100‑1 (MQ=255)
          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                         >  2:1439212/1‑100 (MQ=255)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                      <  1:1211545/100‑1 (MQ=255)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                                    <  1:267284/100‑1 (MQ=255)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                   <  1:478482/100‑1 (MQ=255)
                 cccGCAGAGGCAAAGCGCCACCGGGTATAAAGATTGTTGCCACAGTGAGATTTGTTTATGCCGGATGCGGCGTTAACGACTTATCAGGCCTACAACATCg                                                                                 <  1:1101272/98‑1 (MQ=18)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                 <  1:321228/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                 <  2:371093/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                               <  2:887462/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                               <  1:324218/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                              <  2:604494/100‑1 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                             >  2:232380/1‑100 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                             >  2:1527325/1‑100 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                        <  1:202366/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                       >  2:690466/1‑100 (MQ=255)
                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGATGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                      <  1:550814/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCTACAAAATCGTGCAAATTCAAt                                                                     >  1:726380/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                     <  1:1279620/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                     >  2:976319/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                     <  2:838362/100‑1 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                >  2:1148658/1‑100 (MQ=255)
                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                               <  1:226481/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                           >  2:1485127/1‑100 (MQ=255)
                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑Attt                                                   >  1:1099855/1‑100 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑Atttt                                                  >  1:164209/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCg                                                <  2:849645/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCg                                                >  1:789106/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCg                                                <  1:976327/100‑1 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGt                                               >  1:1486182/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTa                                              >  1:251155/1‑100 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAgg                                            <  2:1187436/100‑1 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAgg                                            <  2:569334/100‑1 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAgg                                            >  1:847838/1‑100 (MQ=255)
                                                     ttaaCAAAGTGCGCTTTGTTTATGCAGGATTCGTCGTTAACGCCTGATCCGGCCTACAAAATCGTTCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGcc                                          <  2:153546/96‑1 (MQ=17)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTg                                        >  2:1325421/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGa                                       >  2:1526449/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGAt                                      >  2:1348016/1‑100 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATa                                     <  1:294522/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATa                                     >  2:131025/1‑100 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATaa                                    >  2:987632/1‑100 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGc                                  <  1:273675/100‑1 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCg                                 >  2:139166/1‑100 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCgg                                <  2:563755/100‑1 (MQ=255)
                                                                  cTTTGTTTATGGCGGATGCGGCGTTAAAGCCTTATCCGGCCTAGAAAAGCGTGCAAATTCAATACATTGCAGA‑‑ATGTTCGTAGGCCTGATAAGCGGAgcg                            <  2:1417513/100‑1 (MQ=17)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAgcg                            >  1:44386/1‑100 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCa                          <  1:103631/100‑1 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATc                        >  1:1572164/7‑100 (MQ=18)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATc                        >  2:485241/1‑100 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCa                       <  1:232381/100‑1 (MQ=255)
                                                                       tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGA‑‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAg                      >  1:260919/6‑100 (MQ=18)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg                      >  1:358572/1‑100 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGG‑‑AGTTATGTAGGCCTGATAAGCGTAGCGCATCAgg                     >  1:1299694/5‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGG‑‑ACTATGTAGGCCTGATAAGCATAGCGCATCAGGc                    >  2:649396/3‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                    >  1:1146109/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATG‑‑AGCCACGTAGGCCTGATAAGCGTAGCGCATCAGGc                    >  1:249526/2‑99 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAG‑‑ATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                   >  2:1130087/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGG‑‑GACTATGTAGGCCTGATAAGCATAGCGCATCAGGCa                   >  2:461992/2‑100 (MQ=18)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                  >  1:1317673/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCG‑‑ATCCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                  >  1:993691/1‑99 (MQ=21)
                                                                             ccggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGA‑‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                >  1:1230649/9‑100 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATG‑‑AGCCACGTAGGCCTGATAAGCGTAGCGCATCATGCagt                 >  1:947137/9‑98 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑‑ACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt               >  1:585039/7‑100 (MQ=18)
                                                                               ggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt               >  1:399837/7‑100 (MQ=18)
                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt               <  1:134947/100‑1 (MQ=255)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑‑ACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt              >  1:831241/6‑100 (MQ=18)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGG‑‑AAATTC‑‑AGGCCTGATAAGCGTAGCGCATCAGGCttttcac             >  1:320365/6‑97 (MQ=37)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg             <  1:345541/100‑1 (MQ=255)
                                                                                  taaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGA‑‑ACGATGTAGGCCTGATAAGCGTAGCGCATCAGGCaactacgt            >  1:1515347/4‑93 (MQ=25)
                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc            >  2:1003362/1‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt         <  2:505645/100‑1 (MQ=255)
                                                                                      gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATCAGCGGAGCGCATCAGGCCTTttttccttt        >  2:131786/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg       <  1:690469/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg       <  1:987149/100‑1 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc     >  2:391425/1‑100 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa    <  2:1300810/100‑1 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc  <  1:592866/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATC  >  NC_000913/3705887‑3706078

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