Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 482,230 G→A 100% 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,2300GA95.3% 145.1 / ‑1.8 51P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (25/22);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_pae_tpiA/482084‑482375
                                                                                                                                                   |                                                                                                                                                 
tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                                 <  1:589374/149‑1 (MQ=255)
  acatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGt                                                                                                                                               >  2:630614/1‑149 (MQ=255)
      aaaCTTTCGTCCACGGACCGCGGTCGATAACGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTTTTAGAAACACCCAGCGCCCGCGCTTTTTCTTGGTCGATATCAATCTTAAACTGCGGGGTATCt                                                                                                                                           <  2:98932/149‑1 (MQ=255)
      aaaCGTTCTTCACACGTCCGCCGTCCATCAACTCGTTCACATAGCCGCCGCCCCATGCAGCGCCCAGAGTGTTGTTAATGTCGTTGATAGAAACACCCAGCTCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACCGCGGGGTATCt                                                                                                                                           <  2:134541/149‑1 (MQ=255)
              ttcACACGACGGCGGTCGATAATGTTGTTCTCATGGCTGCTGCCCCAGGCAGCGCCCAGTGTGGTCTTAATGTCCTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGAcc                                                                                                                                   <  2:46/149‑1 (MQ=255)
              ttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTCTCTTGCCGAcc                                                                                                                                   >  1:520473/1‑149 (MQ=255)
                 acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCCTGCGGCGCCCAGAGTGGTGTTAATGTCGTGGATAGAAACACCCAGCGCCTGCGCGTTTTCCTGATCGCTATCAACATTAACATGCTGGGCATCTTCCACGCCGtt                                                                                                                                >  2:150568/1‑149 (MQ=255)
                 acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCCGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCGTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTGTCTTCCAGCCCGtt                                                                                                                                >  1:597782/1‑149 (MQ=255)
                      gACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                            <  1:36909/149‑1 (MQ=255)
                       aCCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGAc                                                                                                                           <  1:122024/149‑1 (MQ=255)
                           cgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTAc                                                                                                                       <  1:522202/149‑1 (MQ=255)
                                gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGATATCTTCCAGACCGTTTGGGCGGACGCTgg                                                                                                                  >  1:104127/1‑149 (MQ=255)
                                  tAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTc                                                                                                                >  1:26945/1‑149 (MQ=255)
                                       tCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCCGCGCTTTTTCCTGGTCGATATCAATCTTACACTGCAAGATACCTTCCCGCCCGTTGGGACGTCCGCTGGTCAACat                                                                                                           >  2:369373/1‑149 (MQ=255)
                                         gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATATTAAACTCCAGGGTATGTTCCAGCCGGTTTGGACGTACGCTGGTCAActagt                                                                                                         >  1:456825/1‑145 (MQ=255)
                                             acaTGGTTGCCGCGCCACGCAGCGCCCAGAGTGGTGTAAAGGTCGGTGAGAGAAACGCCCAGTGCCTGCGCTTTTTACCGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGCCCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     <  2:364187/149‑1 (MQ=255)
                                             acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATCTCGTTGATAGCTACACTCAGCGCCTGCGCTCTTTCCTGCTCGATATCAATCTTAAACTGCAGGGTATCTTCCATACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     <  2:364500/149‑1 (MQ=255)
                                              gaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGAGATCAATCTTAAACTGCAGGGTATCTTCAGGACCGTTTGGACGTACGCTGGTCAACATATCAggg                                                                                                    <  1:356859/148‑1 (MQ=255)
                                                      ccgtcCCACGCAGCGCCGAGAGTGTTGTTAAGCTCTTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGAGATTAATCTTCAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgc                                                                                            <  2:128370/149‑1 (MQ=255)
                                                      ccgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTccc                                                                                            >  1:148610/1‑147 (MQ=255)
                                                             aTGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCt                                                                                     >  1:292798/1‑149 (MQ=255)
                                                              tGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATGTCAATCTTAAACTGCAGGGTATCTTCCAGACAGTTTGGAAGTACGCTGGTCAACATATCACGGTGCTTCGCTGGTTCTg                                                                                    >  1:558338/1‑149 (MQ=255)
                                                              tGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTg                                                                                    <  1:685902/149‑1 (MQ=255)
                                                                 aGCGCCCATAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCCATCTTACACTCCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                                 <  1:102131/149‑1 (MQ=255)
                                                                 aGCGCCCAGCTTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                                 <  1:400857/149‑1 (MQ=255)
                                                                      ccggagTGGGGTTAATGTAGTGGAGAGAAACACCCAGGGGCTGCGTTTTTTTCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                            <  2:583028/146‑1 (MQ=255)
                                                                        agagTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTg                                                                          >  1:144228/1‑149 (MQ=255)
                                                                               tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                   >  1:473895/1‑149 (MQ=255)
                                                                               tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                   <  2:406421/149‑1 (MQ=255)
                                                                               tgtTAATGTCGTTGATAGAAACACCCAGCGAATCCGCAGGTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                   >  1:729141/1‑149 (MQ=255)
                                                                                   aaTGTCGTTGATAGGAACACCCAGCGGCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                               <  2:518414/149‑1 (MQ=255)
                                                                                    aTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGATCAACATATCAGGGTGCTTCGCTGGTTCTGCAAGCAACTGGTTACGccccgg                                                              >  1:189842/1‑147 (MQ=255)
                                                                                       tCGTTGATAGAAACACCCAGCGCATGCAATCGGTCTGATTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGCGt                                                           >  1:174936/1‑149 (MQ=255)
                                                                                         gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCCACTGAGTCa                                                         >  2:435674/1‑149 (MQ=255)
                                                                                          ttGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAg                                                        <  1:35716/149‑1 (MQ=255)
                                                                                                 aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAAGCGCCTGAGTCGGTTTGTCg                                                 >  1:148449/1‑149 (MQ=255)
                                                                                                  aaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGTTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                <  1:518602/149‑1 (MQ=255)
                                                                                                      accAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACa                                                                                                            >  1:419712/2‑85 (MQ=255)
                                                                                                      accAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACa                                                                                                            <  2:419712/84‑1 (MQ=255)
                                                                                                             ccTGCGCTTGTTCCTGGTCGATATCAATCTTACACTGCAGGGTATCCTCCAGTCCGTTTGGACGTACGCTGGTCAACCTATCAGCGTGCTCCTCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                     <  2:593772/149‑1 (MQ=255)
                                                                                                                  gcTCTTTCCTGGTCGATATCAATCCTAAACTGCAGTGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                <  1:70597/149‑1 (MQ=255)
                                                                                                                              tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtc                    >  1:226847/1‑149 (MQ=255)
                                                                                                                                 atatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa                 >  1:471529/1‑149 (MQ=255)
                                                                                                                                  tatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaag                >  1:229239/1‑149 (MQ=255)
                                                                                                                                        tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCACTCAGCtcaaagtcaaag          >  1:370864/1‑149 (MQ=255)
                                                                                                                                        tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCCAAtcaaagtcaaag          >  2:362924/1‑149 (MQ=255)
                                                                                                                                        tCTTAAACTGCAGGGTATCTTCCAGACCGGTTGGACGTACGCTGCTCAACCTATGAGGGTGGTTCGCTGCGTCTGCAAGCAACTGGCTACGCGCCTGAGTCAGATTTTCGTGACCAAGGCCAGCCAGGTCAATCAGCtcaaagtgaaag          >  1:621380/1‑149 (MQ=255)
                                                                                                                                           tAAGCTGCAGGGTATTGTCCAGACCGTTTGGACGTATGTTGGTCAACATATCGGGGTGCTTCGCTGCTTCTGCAAGCGACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCGAGCCTGGTCAATCAGCTCAAAGTCAAAGCCg       <  2:38347/149‑1 (MQ=255)
                                                                                                                                            aaaCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGCTGCTTCGCTGCTTCTGCAAGCAACTGTTTAGGCGCCTGAGTCAGTTGTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAACCCgg      >  1:211256/1‑149 (MQ=255)
                                                                                                                                              aCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGtt    <  2:144228/149‑1 (MQ=255)
                                                                                                                                              aCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGtt    <  1:303268/149‑1 (MQ=255)
                                                                                                                                                tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGATACGCGCCTGCGTCAGTTTTTCGTGCCCAAGGCCAGCCTAGCCAATCAGCTCAAAGTCAACGCCgactgc  >  2:505966/1‑144 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_pae_tpiA/482084‑482375

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