Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,648,882 +T intergenic (‑119/+177) oprB‑II ← / ← PP_1446 carbohydrate‑selective porin/TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,648,8811.T98.0% 158.7 / ‑4.2 50intergenic (‑118/+178)oprB‑II/PP_1446carbohydrate‑selective porin/TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (25/24);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTAGCAGGTTTTAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGG  >  NC_002947/1648734‑1649020
                                                                                                                                                     |                                                                                                                                           
gATGTAGCAGGTTTTAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGatt                                                                                                                                           >  2:346238/1‑149 (MQ=255)
    tAGCAGGTTTTAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGattatag                                                                                                                                       <  2:764078/150‑6 (MQ=255)
     aGCAGGTTTTAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAgg                                                                                                                                       >  2:392717/1‑150 (MQ=255)
       cAGGTTTTAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAgggg                                                                                                                                     >  1:1036420/1‑150 (MQ=255)
       cAGGTTTTAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAgggg                                                                                                                                     >  2:1440916/1‑150 (MQ=255)
           ttttAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAATACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGGTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGt                                                                                                                                 >  1:1463756/1‑150 (MQ=255)
           ttttAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGatta                                                                                                                                          >  2:904652/1‑138 (MQ=255)
           ttttAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGatta                                                                                                                                          <  1:904652/140‑3 (MQ=255)
             ttAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTta                                                                                                                               <  1:392717/150‑1 (MQ=255)
              tAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTtat                                                                                                                              <  2:928182/150‑1 (MQ=255)
                     ttGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTAAAGGGGAAGTTATAAGaaaa                                                                                                                       >  2:71161/1‑150 (MQ=255)
                            aTGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACaa                                                                                                                <  1:128663/150‑1 (MQ=255)
                            aTGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACaa                                                                                                                >  2:1450950/1‑150 (MQ=255)
                                 aTAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGt                                                                                                           <  2:1017353/150‑1 (MQ=255)
                                 aTAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGt                                                                                                           >  1:1017353/1‑150 (MQ=255)
                                     cgcTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCg                                                                                                         >  2:218000/1‑148 (MQ=255)
                                     cgcTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCg                                                                                                         <  1:218000/148‑1 (MQ=255)
                                            ttttGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAgg                                                                                                                                       <  2:651672/111‑1 (MQ=255)
                                            ttttGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAgg                                                                                                                                       >  1:651672/1‑111 (MQ=255)
                                                gATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTc                                                                                                          >  1:233072/1‑136 (MQ=255)
                                                gATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTc                                                                                                          <  2:233072/136‑1 (MQ=255)
                                                 attattGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTAAAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGAttgttgt                                                                                           <  1:71161/150‑1 (MQ=255)
                                                  ttattaAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGAttgttgtt                                                                                          <  2:1079328/150‑1 (MQ=255)
                                                  ttattGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGAttgttgtt                                                                                          >  2:1352999/1‑150 (MQ=255)
                                                     ttGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGaa                                                                                       >  2:811615/1‑150 (MQ=255)
                                                         aaaaaaCGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACggg                                                                                   <  2:1569774/150‑1 (MQ=255)
                                                                 gCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAgtc                                                                           >  1:563215/1‑150 (MQ=255)
                                                                   aaGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGAttgtt                                                                                             >  2:1320309/1‑130 (MQ=255)
                                                                   aaGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGAttgtt                                                                                             <  1:1320309/130‑1 (MQ=255)
                                                                     gcgcATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCa                                                                       >  1:422623/1‑150 (MQ=255)
                                                                     gcgcATCGAAGGGCGGGCCCCTCAAGGGATATGGCGACAAGATTGTTCGATGAAGCGTTGTAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCa                                                                       <  1:1197667/150‑1 (MQ=255)
                                                                      cgcATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGGAGATGTGAGTATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCAt                                                                      <  1:1450950/150‑1 (MQ=255)
                                                                             aaGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGa                                                               <  2:422623/150‑1 (MQ=255)
                                                                             aaGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGa                                                               >  1:209124/1‑150 (MQ=255)
                                                                              aGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAg                                                              <  2:209124/150‑1 (MQ=255)
                                                                               gggcgggcCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAgg                                                             <  1:1352999/150‑1 (MQ=255)
                                                                                  cgggcCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAg                                                          <  2:993099/150‑1 (MQ=255)
                                                                                    ggcCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAgat                                                        >  2:1421256/1‑148 (MQ=255)
                                                                                     gcCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACaa                                                                                                                <  2:575292/93‑1 (MQ=255)
                                                                                     gcCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACaa                                                                                                                >  1:575292/1‑93 (MQ=255)
                                                                                         cTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGaaaaaa                                                                                                                     <  1:1239580/84‑1 (MQ=255)
                                                                                         cTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGaaaaaa                                                                                                                     >  2:1239580/1‑84 (MQ=255)
                                                                                         cTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAgg                                                   >  2:1172625/1‑150 (MQ=255)
                                                                                          tCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGggagga                                                                                                                                              >  2:533414/1‑58 (MQ=255)
                                                                                          tCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGggagga                                                                                                                                              <  1:533414/58‑1 (MQ=255)
                                                                                             aTGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCg                                               >  1:5593/1‑150 (MQ=255)
                                                                                               ggATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCggg                                             <  1:1421256/150‑1 (MQ=255)
                                                                                               ggATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAAGGGAAGTAGTCAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCggg                                             <  1:811615/150‑1 (MQ=255)
                                                                                                 atatGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGtt                                           >  2:459468/1‑150 (MQ=255)
                                                                                                           aaGCTTGTTCGATGAAGCGTTTCAGGTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGaa                                 <  2:1269249/150‑1 (MQ=255)
                                                                                                                     gATGAAGCGTTTCAGCTTGCGGCGGGGAGAAATATAGGGGAAGTTATAAGAAAAAAGACAAATTTTCGCAGAGGATTGTAGTTGAACGGGCAACACTCCTCATTAATCGAGGCAAGTGTAGAAGAGGGTTAACCAGCAAAAGCGcccatt                       >  2:148894/1‑146 (MQ=255)
                                                                                                                        gAAGCGTTTCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGcc                    >  2:1425276/1‑150 (MQ=255)
                                                                                                                               ttCAGCTTGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAt             >  2:682402/1‑150 (MQ=255)
                                                                                                                                     ttGCGGCGGGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACAATCGTAtt       >  2:1318820/1‑150 (MQ=255)
                                                                                                                                         ggcggGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTcgcg   <  2:774721/150‑1 (MQ=255)
                                                                                                                                          gcggGGAGGATTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg  <  1:1172625/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GATGTAGCAGGTTTTAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGG  >  NC_002947/1648734‑1649020

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