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,2300GA98.0% 156.9 / NA 51P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/25);  total (26/25)
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

ATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_pae_tpiA/482088‑482375
                                                                                                                                              |                                                                                                                                                 
atAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAAAGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTACTGGTCGATATCAATCTTAAACTGCAGGGTAt                                                                                                                                             >  1:129536/1‑149 (MQ=255)
   aaCCTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCtt                                                                                                                                          >  1:508661/1‑149 (MQ=255)
    agtttCTTCACACGACCGCGGTCGATAAAGTCGTTCAAATCGCTGCCGCTCCATGCAGCGCCCAGAGTGGGGCTACAGTCGGTCATACAGGCAGGCAGGGCCGGAAGATTTTCCTGGTCGATATCAAATTTAGACTGCAGTGTATCTTc                                                                                                                                         <  2:362185/147‑1 (MQ=255)
     cTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTcc                                                                                                                                        >  2:398738/1‑149 (MQ=255)
      gttcttcACTCTTCTGCGGTCGATAAAGTCGTTTACATAGCTGCCGCTCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                       <  2:161165/148‑1 (MQ=255)
             acacGACGGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCTTTGATAGAAACACGCAGCGCCTGCGCTTTTTCTTGGGCGCTATCAATCTTAGAATGCGGGGTATCTGCCGGATCGtt                                                                                                                                >  2:237903/1‑149 (MQ=255)
                  aCCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCGATCTTAAACTGCAGGGTATCTTCCGGGCCGTTTGGAc                                                                                                                           >  2:69761/1‑149 (MQ=255)
                  aCCGCGGTCGATAAAGTCGTTCACAGAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCTTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGTAc                                                                                                                           <  2:621957/149‑1 (MQ=255)
                        gTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                               >  2:478175/1‑139 (MQ=255)
                        gTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                               <  1:478175/139‑1 (MQ=255)
                          ctaTAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGGCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTg                                                                                                                   <  2:128907/147‑1 (MQ=255)
                              aaaGTCGTTCACATAGCTGCCGCCCCATGCAGCGACCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCCGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTCTCTTCCAGACCGTTTGGACGTCCCCTCGTCa                                                                                                               >  2:651602/1‑149 (MQ=255)
                               aaGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCaa                                                                                                              <  1:687374/149‑1 (MQ=255)
                                aGTCGTTCCCATAGCTGCAGCCCCATGCAGCGCCCAGAGCGGTGTTATTGTCGTTGATAGAAACACCCAGCGCCTGCGCATTTTCCTGGTCGATATCAATCTTAATCTGCAGGGTATCTTCCATACCGTTTGGACGTACGCTGGTCAAc                                                                                                             <  1:594280/149‑1 (MQ=255)
                                aGTCGTTCACATAGCTGCTGCCCCATGCAGCGCCCAGAGTGGTGTTATTGTCGTTGATAGAAACACCCAGCGCCTGCGCATTTTCCTGGTCGATATCAATCTTAATCTGCAGGGTATCTTCCATACCGTTTGGACGTACGCTGGTCAAc                                                                                                             <  1:594319/149‑1 (MQ=255)
                                        acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     <  1:615767/149‑1 (MQ=255)
                                                          gCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgc                                                                                   <  2:514902/149‑1 (MQ=255)
                                                                cccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaagcaa                                                                             <  1:35639/149‑1 (MQ=255)
                                                                    gagTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTgg                                                                         <  1:398738/149‑1 (MQ=255)
                                                                        ggtaTTAATGTCTTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATAGCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAc                                                                     <  2:123835/149‑1 (MQ=255)
                                                                          tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                   >  2:758106/1‑149 (MQ=255)
                                                                                  tCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGt                                                           >  1:486465/1‑149 (MQ=255)
                                                                                  tCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCCCCTGAGt                                                           >  2:176936/1‑149 (MQ=255)
                                                                                    gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCACTCTTAAACTGCAGGGTATCTTCCCGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCACGCAACTGGTTACGCGCCTGataca                                                         >  2:502192/1‑145 (MQ=255)
                                                                                    gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                         <  1:255904/149‑1 (MQ=255)
                                                                                     ttGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTGGATATCAATGTTAAACGGCAGGGTATCTTCCTGACAGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTTCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAg                                                        >  2:641796/1‑149 (MQ=255)
                                                                                     ttGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAg                                                        >  2:527507/1‑149 (MQ=255)
                                                                                     atcaTAGAAGCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGGAACTGGTTACGCGACTGAGTCAg                                                        >  2:133968/4‑149 (MQ=255)
                                                                                       gATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtt                                                      <  1:463138/149‑1 (MQ=255)
                                                                                        aTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                     >  1:771627/1‑149 (MQ=255)
                                                                                           gAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                  <  1:194438/149‑1 (MQ=255)
                                                                                             aaCACCCAGCGCCCGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                <  2:676359/149‑1 (MQ=255)
                                                                                                acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                             >  2:327189/1‑149 (MQ=255)
                                                                                                  ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCACCTGGTTCCGCGCCTGGGTCAGCTTTTCGTGACCa                                           >  2:178952/1‑149 (MQ=255)
                                                                                                        ccTGCGGTTTTTCCTGGTCGATATGAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                     <  2:115577/149‑1 (MQ=255)
                                                                                                            tgcTTTTTCCTGCTCGATATCAATCTTAAACGGCAGGCTACCTTCCAGACCGGTTGGACGTACGCTTGTCAACATATCAGGGTGCTTCGCTGCTTCTGTAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCt                                 <  1:258767/148‑1 (MQ=255)
                                                                                                            cgcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCt                                 >  1:41414/1‑149 (MQ=255)
                                                                                                               tatttCCTGGTCGATAAAACACTCAAACTGCAGAACATCTTCCAGACCGTTTGGACGTACGCTGGCCAACTTCTCAGACGGCTTCGCTGCTTCTGCAAGCAACCGGTTACGCCCCTACGTCAATTCTTCGACACCAAGGCCCGCCTgct                              >  2:133176/3‑147 (MQ=255)
                                                                                                                 tttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                >  1:15012/1‑145 (MQ=255)
                                                                                                                 tttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                <  2:15012/145‑1 (MQ=255)
                                                                                                                   tCCTGGTCGATATCAAGCTTAATCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                          <  2:557517/149‑1 (MQ=255)
                                                                                                                    ccTGGTCTTTATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATc                         <  1:337517/149‑1 (MQ=255)
                                                                                                                      tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                       >  1:530192/1‑149 (MQ=255)
                                                                                                                         tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTACTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCCGCtc                    >  2:659621/1‑149 (MQ=255)
                                                                                                                          cGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                   >  1:464113/1‑149 (MQ=255)
                                                                                                                            ttagtcaTCTTAAACTGCAGGGTATCGTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTTAATCAGCtcaaa                 <  1:530945/143‑1 (MQ=255)
                                                                                                                              atCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTACTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt               >  2:272651/1‑149 (MQ=255)
                                                                                                                                  aTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaa           <  1:327189/149‑1 (MQ=255)
                                                                                                                                   tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaacg          >  1:322187/1‑147 (MQ=255)
                                                                                                                                    cTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc         <  1:556765/149‑1 (MQ=255)
                                                                                                                                           tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGc  >  2:661468/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
ATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_pae_tpiA/482088‑482375

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