Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,753,755 A→G 6.9% T51A (ACG→GCG)  bamE → outer membrane protein assembly factor BamE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,753,7550AG 6.9% 147.1 / 3.0 60T51A (ACG→GCG) bamEouter membrane protein assembly factor BamE
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base G (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-02

TAAAACGCTGACTGCTGCAGCAGCAGTACTATTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAAC  >  NC_000913/2753613‑2753892
                                                                                                                                              |                                                                                                                                         
tAAAGCGCTGGCTGCTGCAGCAGCAGTACTGTTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGa                                                                                                                                           <  2:866852/143‑1 (MQ=255)
     cGCAGACGGCTGCAGCAGCAGTACTATTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGcaa                                                                                                                                      <  1:395098/143‑1 (MQ=255)
      gCTGACTGCTGCAGCAGCAGTACTATTTATGTGGAGCACAGGCTGTGCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGGGTTGGCATGACCcaac                                                                                                                                     >  1:495098/1‑143 (MQ=255)
        tGACTACTGCAGCAGCAGTACTATTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCAGTACGcaacaa                                                                                                                                   >  2:66748/1‑143 (MQ=255)
         gtcggctgcAGCGGCAGTACGATTCAGGTTGACCGCAGGCTGCTCCACTCTGGAGCGAGTGGTGTACCGTCCTGCCATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGcaacaac                                                                                                                                  <  2:669995/139‑1 (MQ=255)
         gACTGCTGCCGCAGCCGGACTATTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGcaacaac                                                                                                                                  <  2:671198/143‑1 (MQ=255)
          aCTGCTGCAGCAGCAGTACTATTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGcaacacca                                                                                                                                 >  1:384492/1‑143 (MQ=255)
           cggctgcAGCAGCAGTACTATTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGcaacaacaa                                                                                                                                <  2:1040333/141‑1 (MQ=255)
            tgctgcAGCAGCAGTACTATTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAg                                                                                                                               <  1:984164/143‑1 (MQ=255)
            tgctgcAGCAGCAGTACTATTGATGTTGACCGCAGACTGTTCCACTCTGGAGCGAGAGGGTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTACAGCCGCATCCAAAACACGTGCTCGCATGACGCAAAAACAAg                                                                                                                               >  2:383862/1‑143 (MQ=255)
              ctgcAGCAGCAGTACTATTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTTGCATGACGCAACAACAAGtt                                                                                                                             >  2:21624/1‑143 (MQ=255)
                      gcggTACTATTGATGTTGCCCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGc                                                                                                                     <  2:430386/140‑1 (MQ=255)
                          gaCTATTGCTGCTGACCGCATGCTGTTCCACGCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGCGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTg                                                                                                                 <  1:90722/142‑1 (MQ=255)
                              agtGATGTTCACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTa                                                                                                             <  1:707262/141‑1 (MQ=255)
                                tGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACGGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATGGGGTaca                                                                                                           >  1:270139/1‑143 (MQ=255)
                                   tGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCg                                                                                                        >  1:960539/1‑143 (MQ=255)
                                    gTTGACCGCAGGCTGTTCCACTCTCGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGc                                                                                                       <  1:910426/143‑1 (MQ=255)
                                     ttGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCt                                                                                                      >  1:340663/1‑143 (MQ=255)
                                       gACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCAGAGGGTACCCCGCTGa                                                                                                    >  2:391508/1‑143 (MQ=255)
                                          cGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTCACGACGCATCCAAAAGACGTGTTCGCATGACGCAACAACAAGTGCCGAACACAATGTGTACACCGcggccgt                                                                                                 >  2:4299/1‑137 (MQ=255)
                                              gggtGTTCCACTCGGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGa                                                                                             <  2:336969/140‑1 (MQ=255)
                                              ggCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGa                                                                                             <  1:825855/143‑1 (MQ=255)
                                              ggCGGTCCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGa                                                                                             <  1:558643/143‑1 (MQ=255)
                                                     ccACTCTGGAGCGAGTGGCTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCAttt                                                                                      <  1:257536/143‑1 (MQ=255)
                                                           tagAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACg                                                                                <  1:246756/141‑1 (MQ=255)
                                                                cGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATAc                                                                           >  2:927708/1‑143 (MQ=255)
                                                                   gTGGTTTACCGTCCTCACATCAACCAGGGGAACTATCTGACCGCTAACGACGTCTCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCACTACATTTAGTACGAATACCTg                                                                        >  2:888260/1‑143 (MQ=255)
                                                                    tGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTgc                                                                       >  1:451184/1‑142 (MQ=255)
                                                                    tGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGATAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTACGAACAAATTGGGTCCCCAGCCGCTTGCCGAGCGATCTGGGACGAATAccgga                                                                       >  2:506682/1‑140 (MQ=255)
                                                                      gTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGGCGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGtt                                                                     >  2:145405/1‑143 (MQ=255)
                                                                      gTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGGCGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGtt                                                                     >  2:148892/1‑143 (MQ=255)
                                                                          gctgTCCTGTCATCAACCAGGGGAACTATCTGACCGCTACCGACGTATCCACAATACGTTTTGGCATGACGCAACAACAAGTTGCGTACGCATTGCGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTAt                                                                 <  2:706834/140‑1 (MQ=255)
                                                                            cGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCACAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGAACCATTTGGTACGAATAACTAGTGCTATGt                                                               >  2:309684/1‑143 (MQ=255)
                                                                            cGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCACAACACGTGTTGGCATGCCGCAACCACACGTTGCGTACGCATTGCGTACACCCGTGATGTCCGACCCATTTGGTACGAATACCCGGTTCTACGt                                                               >  2:415237/1‑143 (MQ=255)
                                                                            cGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGt                                                               >  2:310860/1‑143 (MQ=255)
                                                                                cTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACCACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGACCCATTGCGTACGAATACCTGGTTCCATGGCTTc                                                           >  2:866130/1‑143 (MQ=255)
                                                                                cTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCCTTGGGTACACCGCTGATGACCGATCCATTTGGTACGAATACCTGGTTCTATgcgttc                                                           >  2:873488/1‑138 (MQ=255)
                                                                                 tGACATCAACCAGGGGAACTATCGGACCGCTACCGACGTATCCAAAGTACGTGTTGGCATGCCGCAACAACAAGTTGCGTACGCATTGGGTACACCGGTGATGTCCGATCCATGTGGTACGAATACCTGGCTCCATGTcaccc                                                          >  2:179538/1‑139 (MQ=255)
                                                                                   aCATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGc                                                        <  1:837747/143‑1 (MQ=255)
                                                                                    gaTCAACCAGGGGAACTAGCTGACCGCTAACGACGGATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGcc                                                       <  2:724795/142‑1 (MQ=255)
                                                                                             ggggAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAg                                              >  1:106548/1‑143 (MQ=255)
                                                                                             gggcAACTATCTTACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACCCCGCTGACGTCCGATCCATTTGATACCAATACTTGGTTCTATGTCTTCCGCCAGCAACCAg                                              <  2:270015/143‑1 (MQ=255)
                                                                                               ggAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGACTACCTGGTTCTATGTCTTCCGCCAGCAACCAGGt                                            >  1:406353/1‑143 (MQ=255)
                                                                                                      tCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAg                                     >  1:277396/1‑143 (MQ=255)
                                                                                                      tCTGACCGCTAACGACGTAGCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCAGTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAg                                     <  2:944127/143‑1 (MQ=255)
                                                                                                            cGCTAACGACGTCTCCTAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTaa                               <  2:680965/143‑1 (MQ=255)
                                                                                                            cGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGCCATGAAGGTGTaa                               >  2:730100/1‑143 (MQ=255)
                                                                                                              cTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACt                             >  2:969653/1‑143 (MQ=255)
                                                                                                                  cgacgTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGc                         <  1:927708/143‑1 (MQ=255)
                                                                                                                        aTCCAAAATACGTGTTGGCATGGCGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAacgc                   <  1:145405/143‑1 (MQ=255)
                                                                                                                        aTCCAAAATACGTGTTGGCATGGCGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATGTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAacgc                   <  1:148892/143‑1 (MQ=255)
                                                                                                                        aTCCAAAATACGTGTTGGCATGACGCAACAGCAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAacgc                   <  1:471048/143‑1 (MQ=255)
                                                                                                                        aTCCAAAAGACGTGTGGGCTTGACGCAACAACAAGTTGCCTACGCATTGTGTACACCGCTGATGTCCGAGCCACTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAacgc                   <  2:720150/143‑1 (MQ=255)
                                                                                                                         gccAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCACTGGGTACACCGCTGATGGCCGATCCATTTGGGACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAacgct                  <  1:320617/142‑1 (MQ=255)
                                                                                                                            aaaaTACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCCGCCACCAGGTCATGAAGGTGTAACTCAGCAAaccctgac               >  2:686855/1‑143 (MQ=255)
                                                                                                                             aaaTACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCACTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATCAAGGTGTAACTCAGCAAacgctgacg              >  2:267694/1‑143 (MQ=255)
                                                                                                                                   gtgtTGGCATGCCGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCACCCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGAcc        <  2:932053/143‑1 (MQ=255)
                                                                                                                                    tgtTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCt       <  1:21624/143‑1 (MQ=255)
                                                                                                                                        ggCATGACGCAACAACAAGTTGCGTACACATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCGTTaa   >  2:657529/1‑143 (MQ=255)
                                                                                                                                        ggCATGACGCAACAACAAGTTGCGTACACATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAACAAACACTGACCCTGACCtttca   >  2:657992/1‑141 (MQ=255)
                                                                                                                                         gCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCCAACGCTGACGCTGACCTTTaaa  >  1:331277/1‑142 (MQ=255)
                                                                                                                                              |                                                                                                                                         
TAAAACGCTGACTGCTGCAGCAGCAGTACTATTGATGTTGACCGCAGGCTGTTCCACTCTGGAGCGAGTGGTTTACCGTCCTGACATCAACCAGGGGAACTATCTGACCGCTAACGACGTATCCAAAATACGTGTTGGCATGACGCAACAACAAGTTGCGTACGCATTGGGTACACCGCTGATGTCCGATCCATTTGGTACGAATACCTGGTTCTATGTCTTCCGCCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAAC  >  NC_000913/2753613‑2753892

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