Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,648,966 A→G 13.1% intergenic (‑245/+51) 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,9660AG13.1% 128.8 / 9.3 53intergenic (‑245/+51)oprB‑II/PP_1446carbohydrate‑selective porin/TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base G (2/5);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCA  >  NC_002947/1648821‑1649114
                                                                                                                                                 |                                                                                                                                                    
cAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATc                                                                                                                                                    <  3:144325/148‑1 (MQ=255)
         gcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCgg                                                                                                                                          <  3:216997/149‑1 (MQ=25)
         gcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCgg                                                                                                                                          <  2:283003/149‑1 (MQ=25)
         gcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAAGAGCGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCgg                                                                                                                                          <  1:268999/149‑1 (MQ=25)
                 gaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAAccc                                                                                                                                  >  2:152113/1‑149 (MQ=255)
                    ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAAc                                                                                                                                    >  6:70766/1‑144 (MQ=255)
                    ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAAc                                                                                                                                    <  5:70766/144‑1 (MQ=255)
                        gggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTAtt                                                                                                                                               <  3:14050/129‑1 (MQ=255)
                        gggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTAtt                                                                                                                                               >  4:14050/1‑129 (MQ=255)
                          gAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTaca                                                                                                                         >  3:297773/1‑149 (MQ=255)
                          gAAGTCATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTaca                                                                                                                         <  1:91818/149‑1 (MQ=255)
                                ataAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGAt                                                                                                                    <  4:297773/148‑1 (MQ=255)
                                    gAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtcg                                                                                                               >  6:93115/1‑149 (MQ=255)
                                     aaaaaaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtcgt                                                                                                              >  4:156507/1‑149 (MQ=255)
                                       aaaaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtcgtcg                                                                                                            >  6:255818/1‑149 (MQ=255)
                                          aGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                                         >  7:185265/1‑149 (MQ=255)
                                            aCAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTc                                                                                                                              >  8:102801/1‑126 (MQ=17)
                                            aCAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTc                                                                                                                              <  7:102801/126‑1 (MQ=17)
                                              aaaTTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAACGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTgg                                                                                                     >  6:138023/1‑149 (MQ=255)
                                              aaaTTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCAACCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTgg                                                                                                     >  7:336066/1‑149 (MQ=255)
                                                   gTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCa                                                                                                 <  8:185265/148‑1 (MQ=255)
                                                       cAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                                         >  6:148380/1‑136 (MQ=255)
                                                       cAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                                         <  5:148380/136‑1 (MQ=255)
                                                            gATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                          >  7:150915/1‑146 (MQ=255)
                                                            gATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                          <  8:150915/146‑1 (MQ=255)
                                                                 ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACAcc                                                                                                                       >  1:101156/1‑112 (MQ=255)
                                                                 ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACAcc                                                                                                                       <  2:101156/112‑1 (MQ=255)
                                                                               aCAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGc                                                                    <  3:156507/149‑1 (MQ=255)
                                                                                 aGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTa                                                                  >  5:282046/1‑149 (MQ=255)
                                                                                   gcgtcTTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCAGCGAATGCGGCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTAc                                                                 <  8:336066/147‑1 (MQ=255)
                                                                                     gtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGcc                                                              >  4:147468/1‑149 (MQ=255)
                                                                                      tcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTcc                                                                                                                             <  4:35532/85‑1 (MQ=255)
                                                                                      tcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTcc                                                                                                                             >  3:35532/1‑85 (MQ=255)
                                                                                         ttAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCgg                                                          >  3:113417/1‑149 (MQ=255)
                                                                                            aTCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAACGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGGGGTACGCAACCCGACTTCCACGCTACGCCCCGgcgc                                                       <  5:138023/149‑1 (MQ=255)
                                                                                                   cAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCg                                                <  1:152113/149‑1 (MQ=255)
                                                                                                    aGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGa                                               <  3:147468/149‑1 (MQ=255)
                                                                                                      gtgtAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGata                                             >  5:177182/1‑149 (MQ=255)
                                                                                                            gaacGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAGCCCGACATCCACGCTACGCCCCGGCGCCGGTGCGATATcgcg                                        >  4:93344/4‑148 (MQ=18)
                                                                                                            gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCgcgc                                       <  6:282046/149‑1 (MQ=255)
                                                                                                             aGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                      <  5:255818/149‑1 (MQ=255)
                                                                                                              gCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAg                                     >  4:104532/1‑149 (MQ=255)
                                                                                                               cGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACt                                                                           <  1:101452/110‑1 (MQ=255)
                                                                                                               cGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACt                                                                           >  2:101452/1‑110 (MQ=255)
                                                                                                                  gTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATc                                 >  5:38662/1‑149 (MQ=255)
                                                                                                                  gTTCACCCGCGAATACGATGGTAGGGCATCCGGCGCATTCGCGGGTGAACCCGTTCCTACTCCGATCTTCGTCGCTTGTGGTTCAGAAGGGGGTAGGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATc                                 <  3:93344/149‑1 (MQ=21)
                                                                                                                    caaCCGGCGCATTCGCGGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCg                                >  5:277243/3‑148 (MQ=255)
                                                                                                                                  gcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCa                 >  4:236072/1‑149 (MQ=255)
                                                                                                                                   ccGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCAc                >  7:138964/1‑149 (MQ=255)
                                                                                                                                      gTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGTCATAGCGCACGGt             >  6:365907/1‑149 (MQ=255)
                                                                                                                                          cccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATc                                 >  6:164058/1‑125 (MQ=255)
                                                                                                                                          cccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATc                                 <  5:164058/125‑1 (MQ=255)
                                                                                                                                               ccATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcggt    <  3:104532/149‑1 (MQ=255)
                                                                                                                                                 aTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcggt    >  6:311478/1‑147 (MQ=255)
                                                                                                                                                 aTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCa  <  6:38662/149‑1 (MQ=255)
                                                                                                                                                 aTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCa  <  5:93115/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCA  >  NC_002947/1648821‑1649114

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