Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 482,230 G→A P725L (CCG→CTG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA482,2300GA96.2% 157.4 / NA 53P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (26/24);  total (27/26)
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

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGG  >  NC_000913_3_pae_tpiA/482084‑482371
                                                                                                                                                  |                                                                                                                                             
tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                             <  2:212637/149‑1 (MQ=255)
tgacagAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                             <  1:274456/149‑1 (MQ=255)
    atAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTAt                                                                                                                                         >  1:190273/1‑149 (MQ=255)
     tAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGTGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATc                                                                                                                                        >  1:468704/1‑149 (MQ=255)
          tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                   >  1:306712/1‑149 (MQ=255)
             cttcCCCCGCCCGCGGGCGATAAAGCCGTTCACATAGCGGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACCCCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAACCTGCAGGGTATCTTCCAGAc                                                                                                                                <  1:209923/149‑1 (MQ=255)
              ttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCGATCTTAAACTGCAGGGTATCTTCCAGAcc                                                                                                                               >  1:403092/1‑149 (MQ=255)
               tcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCg                                                                                                                              >  2:619409/1‑149 (MQ=255)
                  cacGACCGCGGTCGCTAAAGTCCTTCCCATAGCTGCCGCCCCCTGCAGCGCCCAGAGTGGTTTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                           <  1:66320/149‑1 (MQ=255)
                   acGACCGCGGTCGATAAAGTCGTGCACATAGGTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACCCCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTg                                                                                                                          <  2:747037/149‑1 (MQ=255)
                               gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCtag                                                                                                              >  2:795843/1‑147 (MQ=255)
                               gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTgg                                                                                                              >  1:184789/1‑149 (MQ=255)
                               gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCAAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGCAATCTTCCAGACCGATTGGACGCACCCTgg                                                                                                              >  1:295605/1‑149 (MQ=255)
                               gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGGGCTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGCTATCAATCTTCACCTGCACGCGACCTTCCGCACCGTTTGCGCGTagctcgc                                                                                                              >  2:682413/1‑143 (MQ=255)
                                aTAAAGTCGTTCACATAGCGGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGt                                                                                                             <  1:403222/149‑1 (MQ=255)
                                     gTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCGGCTCGGTTTCCTGTTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACa                                                                                                        <  2:340704/149‑1 (MQ=255)
                                        gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACatat                                                                                                     >  1:352721/1‑149 (MQ=255)
                                            acaTAGCTGGCGGCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                 <  2:219163/149‑1 (MQ=255)
                                            acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTCGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                 <  2:341678/149‑1 (MQ=255)
                                                aGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATTAGGGTGc                                                                                             >  1:567845/1‑149 (MQ=255)
                                                 gCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCt                                                                                            <  2:184789/149‑1 (MQ=255)
                                                  cTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATAAGGGTGCtt                                                                                           >  1:821474/1‑149 (MQ=255)
                                                        ccccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                     >  1:632367/1‑101 (MQ=255)
                                                        ccccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                     <  2:632367/101‑1 (MQ=255)
                                                                aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                             <  1:480112/149‑1 (MQ=255)
                                                                     cccgagTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTGGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                        <  1:888094/146‑1 (MQ=255)
                                                                       agggtggtgTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTg                                                                      <  2:196956/146‑1 (MQ=255)
                                                                           tggtgTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCACTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTa                                                                  <  1:689984/149‑1 (MQ=255)
                                                                                   aTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTg                                                          >  1:712572/1‑149 (MQ=255)
                                                                                    tGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTTCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGa                                                         <  2:177066/149‑1 (MQ=255)
                                                                                         ttGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGGTATCAGTCTTAAACTGCAGGTTATCTTCCAGACCGTTTGGACGTCCGCTGGTCAACAACTCCGGCGGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAg                                                    >  2:573/1‑149 (MQ=255)
                                                                                         ttGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAg                                                    >  1:658437/1‑149 (MQ=255)
                                                                                          tGATAGAAAAACCCAGCGCCTGCGCCGTTTCCTGGTCGATATCAATCTAAAACTGCCGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATAACAGGGTGCTTCTCTGCTTCTGCAAGCAAATCGTTACGCCCCTGAGTCAGt                                                   <  2:904411/149‑1 (MQ=255)
                                                                                            aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                 >  1:271659/1‑149 (MQ=255)
                                                                                               aaaaCACCCACCTCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                              <  1:85528/148‑1 (MQ=255)
                                                                                                aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCACTCTTAAGCTGCAGGGTATCTTCCAGACCGTTTGGCCGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGTCTGAGTAAGTTTTTCg                                             >  1:462258/1‑149 (MQ=255)
                                                                                                    acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                         <  2:587294/149‑1 (MQ=255)
                                                                                                      ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                       >  1:839012/1‑149 (MQ=255)
                                                                                                       cAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCaa                                      >  1:265990/1‑149 (MQ=255)
                                                                                                                 gcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                            <  2:354122/149‑1 (MQ=255)
                                                                                                                         cTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCCGCCTGGTCAATCa                    >  2:489193/1‑149 (MQ=255)
                                                                                                                          tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCATTCAg                   >  2:833983/1‑149 (MQ=255)
                                                                                                                          tGGTCGATATCAATCTTAAACTGCAGCGAATTTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACGGGTTACACGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                   >  2:675269/1‑149 (MQ=255)
                                                                                                                            gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                 >  1:128340/1‑149 (MQ=255)
                                                                                                                                gtatCAATCTTAAACTGCAGGGTATTTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaat             <  2:734491/148‑2 (MQ=255)
                                                                                                                                atatCACACTTACACTGCAGGGTATCTTCCCGACCGTTTGCACCGCCGCTGGCTGACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa             <  1:27041/149‑1 (MQ=255)
                                                                                                                                atatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa             <  1:801454/149‑1 (MQ=255)
                                                                                                                                ataaCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa             <  2:190273/149‑1 (MQ=255)
                                                                                                                                   tCAAGCTTAAACTTCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATGAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACAAAGGCCAGCCTGGTCAACCACCtcaaagtc          >  2:836749/1‑149 (MQ=255)
                                                                                                                                     aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCCAGGCCAGCCTGGTCAATCAGCgcaaagtcaa        >  2:547981/1‑149 (MQ=255)
                                                                                                                                       tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATGAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag      <  2:712572/149‑1 (MQ=255)
                                                                                                                                          tAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCATCATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTCGTCAATCAGCTCAAAGTCAAAGCCg   >  2:627465/1‑149 (MQ=255)
                                                                                                                                          tAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCg   <  1:852974/149‑1 (MQ=255)
                                                                                                                                           aaaCTGCGGGTGAGCTTCCAGACCGTTTGGACGTACGCTGGTCAACATTTCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATTAGCTCAAAGTCAAAGCCgg  <  2:461092/149‑1 (MQ=255)
                                                                                                                                           aaaCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg  <  2:461573/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGG  >  NC_000913_3_pae_tpiA/482084‑482371

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