Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,647,927 G→A Y265Y (TAC→TAT)  oprB‑II ← carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,647,9270GA94.6% 163.9 / ‑0.0 56Y265Y (TAC→TAT) oprB‑IIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (29/24);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

GCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGAAG  >  NC_002947/1647790‑1648071
                                                                                                                                           |                                                                                                                                                 
gCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCgg                                                                                                                                         >  2:97808/1‑148 (MQ=255)
  cAGGGAGGTGAGCTGCTGCTGGGCGCCAATCCACAAGCGGTGCTTGCTCGTACTGCTGCGGTAGGGGGCGCCGCGGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGATCATCT‑‑TGTGCGTTGGCATTACTGTAATAGTAGCCGGCAc                                                                                                                                      <  2:316811/149‑1 (MQ=255)
  cAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGTCCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCAc                                                                                                                                      <  3:100337/149‑1 (MQ=255)
  cAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCATGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCAc                                                                                                                                      <  1:279918/149‑1 (MQ=255)
           gACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTccc                                                                                                                             >  2:118672/1‑149 (MQ=255)
              ctgctgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCt                                                                                                                           <  6:212994/148‑1 (MQ=255)
                 ctgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGAGCGTTGCTGTCCTTGAGGACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCttt                                                                                                                         >  8:212548/1‑147 (MQ=255)
                     tgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGc                                                                                                                    >  5:340284/1‑148 (MQ=255)
                        gcgcCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCtt                                                                                                                 <  6:340284/148‑1 (MQ=255)
                         cgcCAATCCACCAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGTCATTACTGTAATAGTAGCGGTCACGATATTCCCCTTTCAAGCCTTgg                                                                                                               <  7:212548/149‑1 (MQ=255)
                          gcCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGa                                                                                                              <  2:391972/149‑1 (MQ=255)
                          gcCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGa                                                                                                              >  8:149189/1‑149 (MQ=255)
                                  acaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTCCTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTg                                                                                                        <  1:118672/149‑1 (MQ=255)
                                  acaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                      >  2:294167/1‑149 (MQ=255)
                                  acaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTgg                                                                                                       >  5:283210/1‑148 (MQ=255)
                                       ccGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCa                                                                                                 <  7:149189/149‑1 (MQ=255)
                                         gTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATa                                                                                               <  1:294167/149‑1 (MQ=255)
                                              tgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGCCCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAg                                                                                           >  6:35512/1‑148 (MQ=255)
                                                   gAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGaa                                                                                     >  3:283699/1‑149 (MQ=255)
                                                   gAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCAGTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGaa                                                                                     >  7:251363/1‑149 (MQ=255)
                                                    aaCTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAAc                                                                                    <  4:283699/149‑1 (MQ=255)
                                                       tgctgcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGATCATCT‑‑TGTGCCTGGGCATTACTGTAAGTGTAGCCGGCACGATATTACCCTTTCAAGCCTTGGATACGTGGGGTCCATGCCAGTTCGAACGGc                                                                                 >  3:358071/1‑149 (MQ=255)
                                                             ggTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCAtt                                                                              <  6:283210/146‑1 (MQ=255)
                                                                 ggcggcGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCgacc                                                                       >  2:196898/1‑146 (MQ=255)
                                                                    ggcgCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCt                                                                    >  2:347815/1‑149 (MQ=255)
                                                                    ggcgCCGCTGAGGGCGGCCGGCTGACCGGTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTTGTGTCCATAc                                                                                              >  1:262654/1‑123 (MQ=255)
                                                                                     ccggcTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTgccgctgccgct                                                   >  2:396502/1‑148 (MQ=255)
                                                                                        gcTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGc                                               >  7:287769/1‑149 (MQ=255)
                                                                                                   cTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAAcc                                       <  2:148543/146‑1 (MQ=255)
                                                                                                   cTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAAcc                                       >  1:148543/1‑146 (MQ=255)
                                                                                                     gTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGaaa                                         <  4:28072/142‑1 (MQ=255)
                                                                                                     gTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGaaa                                         >  3:28072/1‑142 (MQ=255)
                                                                                                        cTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCg                                >  8:115924/1‑148 (MQ=255)
                                                                                                        cTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACAGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGcag                                >  4:56804/1‑146 (MQ=255)
                                                                                                         ttGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGa                              >  3:59373/1‑149 (MQ=255)
                                                                                                          tGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGa                                                                                                              >  5:52611/1‑69 (MQ=255)
                                                                                                          tGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGAt                             >  8:91531/1‑149 (MQ=255)
                                                                                                           gagaACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGAt                                                                                                             <  6:52611/69‑1 (MQ=255)
                                                                                                            agaACATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGAtt                            <  2:391136/148‑1 (MQ=255)
                                                                                                                cATCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAgg                       >  1:281618/1‑149 (MQ=255)
                                                                                                                  tCT‑‑TGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGtt                     >  3:72769/1‑149 (MQ=255)
                                                                                                                        gtgCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTgg                   >  7:317445/1‑147 (MQ=255)
                                                                                                                        gtgCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGtt                     <  2:281618/145‑1 (MQ=255)
                                                                                                                          gCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGaa                 <  6:34882/147‑1 (MQ=255)
                                                                                                                           ccTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGTGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAg                <  1:97808/148‑1 (MQ=255)
                                                                                                                               ggCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCTTACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGtt            <  1:196898/147‑1 (MQ=255)
                                                                                                                               ggCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGtt            <  7:115924/147‑1 (MQ=255)
                                                                                                                                  aTTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCt       >  4:159476/1‑149 (MQ=255)
                                                                                                                                   ttACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTc      <  4:72769/149‑1 (MQ=255)
                                                                                                                                    tACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTc      >  8:301381/1‑148 (MQ=255)
                                                                                                                                      cTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTc      >  5:13746/1‑146 (MQ=255)
                                                                                                                                      cTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCg     >  4:75614/1‑147 (MQ=255)
                                                                                                                                      cTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGaa   <  8:251363/149‑1 (MQ=255)
                                                                                                                                       tGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTTAAAGCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGAAg  <  3:56804/149‑1 (MQ=255)
                                                                                                                                       tGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCg     >  2:259103/1‑146 (MQ=255)
                                                                                                                                       tGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCg     >  1:835/1‑146 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCT‑‑TGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTG‑GGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGAAG  >  NC_002947/1647790‑1648071

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: ‑