Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,584,487 A→T 11.6% L176Q (CTG→CAG)  PP_1388 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,584,4870AT11.6% 109.1 / 6.0 43L176Q (CTG→CAG) PP_1388EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base A (19/19);  new base T (3/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

GCGTGGCCCTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCAC  >  NC_002947/1584340‑1584632
                                                                                                                                                   |                                                                                                                                                 
gCGTGGCCCTGGCCCACCAAGGTGATGCCCAGAAGCAGGATGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGACCACCTGTGCCTCACGACGCTGAACTGGCATGTCAGCGAGCGCACGGCGACTGGCCCACAAAGCAACCAg                                                                                                                                                  >  8:279925/1‑149 (MQ=255)
 cGTGGCCCTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTTCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGc                                                                                                                                                 <  1:272547/149‑1 (MQ=255)
   tggccctggcccTCCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGccc                                                                                                                                               <  4:36736/149‑1 (MQ=255)
       cctggcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGc                                                                                                                                           <  8:213268/149‑1 (MQ=255)
           gcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCAACTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCa                                                                                                                                       <  6:318787/149‑1 (MQ=255)
             ccaccaTGGTGATGCCCAGCAGCAGGCTGCCATGGCGGTGATTCAACAGTTTGTCGCCAAGGTAGTCCACCAGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAgg                                                                                                                                     <  1:119168/149‑1 (MQ=255)
               accaAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGtt                                                                                                                                   <  1:83979/149‑1 (MQ=255)
                 caAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGa                                                                                                                                 >  6:195788/1‑149 (MQ=255)
                 caAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCAAGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGa                                                                                                                                 <  6:213486/149‑1 (MQ=255)
                          gCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAaca                                                                                                                        >  6:102086/1‑149 (MQ=255)
                           cccAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCATGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCTGCCGCTGGGGCTGGTTgctttgtg                                                                                                                           >  1:264414/1‑138 (MQ=255)
                           cccAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCTGCCGCTGGGGCTGGTTgctttgtg                                                                                                                           <  2:264414/145‑8 (MQ=255)
                                    agGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAgcgc                                                                                                                >  4:336944/1‑147 (MQ=255)
                                    agGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAgcgc                                                                                                                <  3:336944/147‑1 (MQ=255)
                                                  gagaATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCgg                                                                                                >  2:47292/1‑149 (MQ=255)
                                                                  gcgcCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAg                                                                                <  5:102086/149‑1 (MQ=255)
                                                                    gcCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAgg                                                                                                       >  6:261438/1‑124 (MQ=255)
                                                                    gcCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAgg                                                                                                       <  5:261438/124‑1 (MQ=255)
                                                                                      tgCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGC‑GTGGCCGGGCCagtgctgg                                                           <  5:195788/149‑8 (MQ=255)
                                                                                      tgCCTCACGACGCTGCACTGCCATGTCAGCGAGCGCTCGGCGAATGGCCCACAATGCAACCAGGCCCAGCGGCAGGTTGATGCAGAACACCCAGCGCCACGACAGGTACTCGGGCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTg                                                             <  7:279925/148‑1 (MQ=255)
                                                                                        ccTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAg                                                                                                        >  3:41877/1‑103 (MQ=255)
                                                                                        ccTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAg                                                                                                        <  4:41877/103‑1 (MQ=255)
                                                                                        ccTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGcc                                                          >  4:291560/1‑149 (MQ=255)
                                                                                          tCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAgccgcc                                                                                   >  5:219176/1‑122 (MQ=255)
                                                                                          tCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAgccgcc                                                                                   <  6:219176/122‑1 (MQ=255)
                                                                                            acgacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCAc                                                                 <  6:262181/138‑1 (MQ=255)
                                                                                            acgacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCAc                                                                 >  5:262181/1‑138 (MQ=255)
                                                                                                            aTGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGctgctg                                        >  4:357785/1‑147 (MQ=255)
                                                                                                                            ggCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAGGTAGCCCTGGTAGCgg                      >  6:259401/1‑149 (MQ=255)
                                                                                                                               gAATGGCCCACAAAGCAACCAGCCCCAGCTGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACTCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTGGCGGCcg                   <  1:218578/149‑1 (MQ=255)
                                                                                                                               gAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCcg                   >  8:48100/1‑149 (MQ=255)
                                                                                                                               gAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCcg                   >  8:42212/1‑149 (MQ=255)
                                                                                                                                aaTGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCcgc                  <  3:291560/149‑1 (MQ=255)
                                                                                                                                       cacaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcgcg           >  5:125959/1‑149 (MQ=255)
                                                                                                                                        acaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCTGCCGCGTTcgcgc          >  8:160125/1‑149 (MQ=255)
                                                                                                                                         caAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcgcgcg         >  5:35055/1‑149 (MQ=255)
                                                                                                                                           gatcAACCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGc                                                                    >  2:10114/4‑88 (MQ=37)
                                                                                                                                           gatcAACCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGc                                                                    <  1:10114/85‑1 (MQ=37)
                                                                                                                                            aGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCgggg      <  1:47292/149‑1 (MQ=255)
                                                                                                                                             gCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCggggg     <  6:35055/149‑1 (MQ=255)
                                                                                                                                              cAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCggggg     >  7:184370/1‑148 (MQ=255)
                                                                                                                                               aaCCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGGGTACATGCGGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCggggggg   >  4:306368/1‑147 (MQ=255)
                                                                                                                                                 ccAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCAc  >  2:192597/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCGTGGCCCTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCAC  >  NC_002947/1584340‑1584632

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