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,2300GA92.2% 139.1 / 0.1 51P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (20/27);  total (22/29)
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

ACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTG  >  NC_000913_3_pae_tpiA/482092‑482374
                                                                                                                                          |                                                                                                                                                
aCTTGCTTCACCCGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                        <  2:431800/149‑1 (MQ=255)
  tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAACATTAAAATGCAGCGTATCTTTCa                                                                                                                                      <  2:425682/149‑1 (MQ=255)
         acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                        >  2:405801/1‑140 (MQ=255)
         acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                                        <  1:405801/140‑1 (MQ=255)
                    gTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGc                                                                                                                    <  1:254766/149‑1 (MQ=255)
                      cGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTg                                                                                                                  <  1:16241/149‑1 (MQ=255)
                       gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCTCCCAGAGTGGTGTTAATGTCGTTGATATAAACACCCAGCGCCTGCGCTTTTACCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGCCGTACGCTgg                                                                                                                 >  1:239827/1‑149 (MQ=255)
                       gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTgg                                                                                                                 >  1:430031/1‑149 (MQ=255)
                       gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTgg                                                                                                                 >  2:110748/1‑149 (MQ=255)
                         tAAAGTCGTTCACACAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTc                                                                                                               <  1:87484/149‑1 (MQ=255)
                                  tCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCa                                                                                                      <  1:67773/149‑1 (MQ=255)
                                    acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                    >  1:629969/1‑149 (MQ=255)
                                      aTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGt                                                                                                  >  1:80741/1‑149 (MQ=255)
                                                ccccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgctgc                                                                                        >  2:77335/1‑149 (MQ=255)
                                                ccccATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgctgc                                                                                        <  1:538001/149‑1 (MQ=255)
                                                     tGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTg                                                                                   <  1:665258/149‑1 (MQ=255)
                                                      gCAGCGCCCAGAGGGGTGGTAATGGCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgc                                                                                  <  2:691472/149‑1 (MQ=255)
                                                       cAGCGCCCAGAGTTGTGTTAATGTAGTTGATAGAAACACCTAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGCAAGTAGGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgca                                                                                 <  2:304152/149‑1 (MQ=255)
                                                        aGCGCCCAGGGTTGTGTTAATGGCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                                <  2:585174/149‑1 (MQ=255)
                                                        aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                                <  2:80741/149‑1 (MQ=255)
                                                             ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                           <  2:629969/149‑1 (MQ=255)
                                                             ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTACCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                           >  2:630611/1‑149 (MQ=255)
                                                                  ggggtgTTAATGTCGTTGATAGAACCTCCCAGCGCCTGCTCTTTTTCCTGTTCGATATCAATCTTAAACTGCAGGGTATCTTCCATACCGTTTGGACGGACGCTGGTCACCATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGtt                                                                      <  2:345286/147‑1 (MQ=255)
                                                                           aTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAAACTTACACTGCAGGGTATCTTCCAGACCGTTTGACCGTACGCTGGTCAACATATCAGAGTGCTTCGCTGCTCCTCCACACCACTGGTTCCGCCCCTg                                                             >  2:162514/1‑149 (MQ=255)
                                                                                gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGATTCTGCAAGCAACTCGTTACGCGCCTCGGTCa                                                        >  2:637544/1‑149 (MQ=255)
                                                                                       gAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                 <  1:413835/149‑1 (MQ=255)
                                                                                            acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAACCTATGGCGGATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                            >  1:635118/1‑149 (MQ=255)
                                                                                            acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCGGTTTTTCGTGAc                                            >  1:521869/1‑149 (MQ=255)
                                                                                              ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                          >  1:704567/1‑149 (MQ=255)
                                                                                                 gcgcCTGCTCTTTTTCCTGGTCNATANCAGTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGAACGTACGCTGGTCAACATCTCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAgg                                       <  1:633572/149‑1 (MQ=255)
                                                                                                 gcgcCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAgg                                       <  1:110748/149‑1 (MQ=255)
                                                                                                      tGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGc                                  <  1:393367/149‑1 (MQ=255)
                                                                                                         gcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAATCAGTTTTTCGTGACCCAGGCCAGCCTg                               >  1:689656/1‑149 (MQ=255)
                                                                                                          cTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTgg                              >  1:401466/1‑149 (MQ=255)
                                                                                                               tCCTGGTCGATATCAATCTTAAACTTCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                         <  2:133603/149‑1 (MQ=255)
                                                                                                                ccTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTGGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATc                        <  1:197470/149‑1 (MQ=255)
                                                                                                                 cTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCa                       >  1:552817/1‑149 (MQ=255)
                                                                                                                 cTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCa                       <  2:704567/149‑1 (MQ=255)
                                                                                                                     tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGt                                                                                                  <  2:584229/70‑1 (MQ=255)
                                                                                                                     tCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGt                                                                                                  >  1:584229/1‑70 (MQ=255)
                                                                                                                        atatCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa                <  2:528608/149‑1 (MQ=255)
                                                                                                                        atatCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGCCCAGCCTGGTCAATCAGCtcaaa                >  1:528608/1‑149 (MQ=255)
                                                                                                                        agatCAATCTTGAGCTGCAGGGTATCTTCGAGACCGTGTGGACGTACGCTGGTCAACATATCAGGGTGCGTCGCTGCTTCTGCAAGCAGCTGGTTACGCGCTTGAGTCAGTTGTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtaaaa                <  2:635118/147‑1 (MQ=255)
                                                                                                                             aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                         <  2:407811/135‑1 (MQ=255)
                                                                                                                             aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                         >  1:407811/1‑135 (MQ=255)
                                                                                                                                 ttAAACTACAGGGTATCTTCCAGACCGGTTGGACTTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAACGTCAAAGcc       >  2:702056/1‑149 (MQ=255)
                                                                                                                                   aaaCTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg     <  1:271804/149‑1 (MQ=255)
                                                                                                                                      cTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg  >  2:252308/1‑149 (MQ=255)
                                                                                                                                      cTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg  <  1:252308/149‑1 (MQ=255)
                                                                                                                                       tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                          <  1:622095/124‑1 (MQ=255)
                                                                                                                                       tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                          >  2:622095/1‑124 (MQ=255)
                                                                                                                                          |                                                                                                                                                
ACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTG  >  NC_000913_3_pae_tpiA/482092‑482374

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