Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,584,499 A→T 12.5% L172Q (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,4990AT12.5% 104.7 / 5.6 40L172Q (CTG→CAG) PP_1388EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base A (19/16);  new base T (3/2);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGTCGCCAAT  >  NC_002947/1584351‑1584644
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                                    <  4:41877/103‑1 (MQ=255)
                                                                             ccTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAg                                                                                                                    >  3:41877/1‑103 (MQ=255)
                                                                             ccTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGcc                                                                      >  4:291560/1‑149 (MQ=255)
                                                                               tCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAgccgcc                                                                                               <  6:219176/122‑1 (MQ=255)
                                                                               tCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAgccgcc                                                                                               >  5:219176/1‑122 (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:42212/1‑149 (MQ=255)
                                                                                                                    gAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCcg                               >  8:48100/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)
                                                                                                                                           cccAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTa                                               >  2:233867/1‑110 (MQ=255)
                                                                                                                                                 ggCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGTCGCCAAt  <  8:184370/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGTCGCCAAT  >  NC_002947/1584351‑1584644

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