Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 483,580 T→C Y275C (TAC→TGC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913483,5800TC94.9% 131.3 / ‑0.8 39Y275C (TAC→TGC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (13/24);  total (14/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 = 8.37e-01

CAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGC  >  NC_000913/483441‑483724
                                                                                                                                           |                                                                                                                                                
cAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTcacca                                                                                                                                       <  1:139399/151‑1 (MQ=255)
     ccGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCAccc                                                                                                                                    >  1:454479/1‑149 (MQ=255)
       gACGGAAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGCCCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCt                                                                                                                                <  2:26745/151‑1 (MQ=255)
               gaaCGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTc                                                                                                                        >  1:92542/1‑151 (MQ=255)
                        cATCTTCGCCAGTTCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCGCCACCCAGCTCAATCTTCGCGACGTc                                                                                                                <  2:567999/150‑1 (MQ=255)
                        cATCTTCGCCAGTTCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCa                                                                                                               <  2:527999/151‑1 (MQ=255)
                             tCGCCAGTTCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACg                                                                                                             <  2:269318/148‑1 (MQ=255)
                                    ttCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCAcg                                                                                                   >  1:467238/1‑151 (MQ=255)
                                     tCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCAcg                                                                                                   <  1:80380/150‑1 (MQ=255)
                                      cAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCAcgc                                                                                                  <  2:116028/150‑1 (MQ=255)
                                       aGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCgg                                                                                                <  1:63577/151‑1 (MQ=255)
                                              aTTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCAt                                                                                         >  1:618566/1‑151 (MQ=255)
                                                 gCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCt                                                                                      <  2:457575/151‑1 (MQ=255)
                                                 gCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCGTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCt                                                                                      <  2:212223/151‑1 (MQ=255)
                                                  ccGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCt                                                                                      <  1:322881/150‑1 (MQ=255)
                                                  ccGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTg                                                                                     <  2:44616/151‑1 (MQ=255)
                                                  ccGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTg                                                                                     >  1:178275/1‑151 (MQ=255)
                                                      aGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTc                                                                                 <  1:579488/151‑1 (MQ=255)
                                                               cAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTcagc                                                                        >  2:146415/1‑151 (MQ=255)
                                                                   gcgTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAgg                                                                     >  2:269040/1‑150 (MQ=255)
                                                                     gTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGAtt                                                                  <  1:611120/151‑1 (MQ=255)
                                                                          cACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGcc                                                             <  1:622041/151‑1 (MQ=255)
                                                                           aCCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCg                                                            <  2:242432/151‑1 (MQ=255)
                                                                                     aGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTc                                                    >  1:326365/1‑149 (MQ=255)
                                                                                              cccAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCa                                         <  2:436138/151‑1 (MQ=255)
                                                                                               ccAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAg                                        <  2:278016/151‑1 (MQ=255)
                                                                                                    ccGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGAcgcg                                   >  1:536401/1‑151 (MQ=255)
                                                                                                      ggAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTc                                 <  1:146418/151‑1 (MQ=255)
                                                                                                       gAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCt                                >  1:416205/1‑151 (MQ=255)
                                                                                                         aGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGa                              <  1:326824/151‑1 (MQ=255)
                                                                                                         aGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGa                              <  1:512249/151‑1 (MQ=255)
                                                                                                                 ggCCGTTAAACTCTGCGATGATGTCGTAGTTCTCTCCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCaataat                      <  1:252956/151‑1 (MQ=255)
                                                                                                                   ccGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCaataat                      >  1:142674/1‑149 (MQ=255)
                                                                                                                   ccGTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATag                    <  1:546196/151‑1 (MQ=255)
                                                                                                                     gTTAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATagag                  <  2:546664/151‑1 (MQ=255)
                                                                                                                       tAAACTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAgg                 <  2:199777/150‑1 (MQ=255)
                                                                                                                        aaaCTCTGCGATGATGTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCg               >  2:60582/1‑151 (MQ=255)
                                                                                                                                    gatgTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTgg   >  2:309074/1‑151 (MQ=255)
                                                                                                                                     atgTCGCAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGc  >  2:9124/1‑151 (MQ=255)
                                                                                                                                           |                                                                                                                                                
CAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCGCAGCGGTATCCAGCGCGTTTGCACCGGTCGCCAGCTTGATCCCCAGACCGGAAGCCGGTTGGCCGTTAAACTCTGCGATGATGTCGTAGTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGC  >  NC_000913/483441‑483724

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑