Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,641,460 C→G 11.6% Q375E (CAG→GAG)  PP_4951 → EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,641,4600CG11.6% 137.4 / 6.7 49Q375E (CAG→GAG) PP_4951EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (18/25);  new base G (2/4);  total (20/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

TCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTACTGC  >  NC_002947/5641314‑5641605
                                                                                                                                                   |                                                                                                                                                 
ccgtaCGGTGCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGCCATCCTGCTGGCAAGCCAGGGCCTGGTGGACGTGCACACACCAAATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAg                                                                                                                                                 <  4:350406/145‑1 (MQ=255)
  gAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGtt                                                                                                                                               <  7:117655/149‑1 (MQ=255)
  gAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGtt                                                                                                                                               >  3:103745/1‑149 (MQ=255)
      ggTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCAT‑GGAGTTCACt                                                                                                                                           <  2:126553/149‑1 (MQ=255)
      ggTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACt                                                                                                                                           <  2:170490/149‑1 (MQ=255)
      ggTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACt                                                                                                                                           <  5:218054/149‑1 (MQ=255)
        ttCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCAT‑GGAGTTCACTGc                                                                                                                                         <  6:116660/149‑1 (MQ=255)
          cGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCg                                                                                                                                        <  7:178906/148‑1 (MQ=255)
           ggCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGAt                                                                                                                                      <  3:305250/149‑1 (MQ=255)
           ggCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGAt                                                                                                                                      <  5:192812/149‑1 (MQ=255)
            gCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATg                                                                                                                                     <  4:293854/149‑1 (MQ=255)
              tATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAAATCCAT‑GCAGTTCACTGCGATGaa                                                                                                                                   >  7:218327/1‑149 (MQ=255)
                tCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAaca                                                                                                                                 >  2:118815/1‑149 (MQ=255)
                  gCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAAcaca                                                                                                                               <  2:347848/149‑1 (MQ=255)
                           tCACCGGTAATACCTTACTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCAT‑GGAGTTCACTg                                                                                                                                          >  3:16258/1‑129 (MQ=255)
                           tCACCGGTAATACCTTACTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCAT‑GGAGTTCACTg                                                                                                                                          <  4:16258/129‑1 (MQ=255)
                               cGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATCGCAGTTCACTGCGATGAACACAGTGACGCTCATc                                                                                                                   <  1:50307/149‑1 (MQ=255)
                                 gTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGAt                                                                                                                <  8:100729/149‑1 (MQ=255)
                                        cTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTcgacg                                                                                                         >  7:205205/1‑149 (MQ=255)
                                            cTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACACTCATCGATCTCGACGAcgc                                                                                                     <  4:103745/149‑1 (MQ=255)
                                                        ctgctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCcagcag                                                                                         >  7:181969/1‑149 (MQ=255)
                                                           ctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCgg                                                                                      <  8:181969/149‑1 (MQ=255)
                                                           ctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCg                                                                                       <  8:205205/148‑1 (MQ=255)
                                                             ggCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCTTCGATCTCGACGACGCCAGCGCcagcaggctgt                                                                                    >  2:18817/1‑144 (MQ=255)
                                                             ggCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATc                                                                                                                   <  7:338080/118‑1 (MQ=255)
                                                             ggCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATc                                                                                                                   >  8:338080/1‑118 (MQ=255)
                                                                aaGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTaa                                                                                   >  6:191251/1‑147 (MQ=255)
                                                                aaGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACa                                                                                 <  1:118815/149‑1 (MQ=255)
                                                                       ggccgggGGGACGAGCAGACCCCTTAATGGCGGCTGTCGGTGCAACTTGGCCTTCTGGGCGCAGGCAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGACAGCGCCAGCAGCGGTAACAGCctgct                                                                          <  8:190795/149‑1 (MQ=255)
                                                                        gcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCctgctg                                                                         >  7:323584/1‑149 (MQ=255)
                                                                            ggtggACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCgg                                                                     >  4:61991/1‑149 (MQ=255)
                                                                                gACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCggtgg                                                                  <  3:61991/148‑1 (MQ=255)
                                                                                 aCGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTcgac                                                                                                          >  8:278721/1‑107 (MQ=255)
                                                                                 aCGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTcgac                                                                                                          <  7:278721/107‑1 (MQ=255)
                                                                                         aCCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgcca                                                                                                   >  1:323311/1‑106 (MQ=255)
                                                                                         aCCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgcca                                                                                                   <  2:323311/106‑1 (MQ=255)
                                                                                         aCCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAAc                                                         <  8:242347/148‑1 (MQ=255)
                                                                                         aCCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAAc                                                         >  7:242347/1‑148 (MQ=255)
                                                                                                       gctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTgg                                          >  6:182153/1‑149 (MQ=255)
                                                                                                          gTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCg                                                                                       <  7:313094/101‑1 (MQ=255)
                                                                                                          gTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCg                                                                                       >  8:313094/1‑101 (MQ=255)
                                                                                                           ttGGTGCACATTGTCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGgtgt                                      <  8:323584/149‑1 (MQ=255)
                                                                                                              gTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGc                                   <  3:141057/149‑1 (MQ=255)
                                                                                                                           ttCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCa                      >  4:238959/1‑149 (MQ=255)
                                                                                                                                atcgcAGTGAACTGCAT‑GGAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGtt                 >  8:142525/3‑149 (MQ=255)
                                                                                                                                 gcgcAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTg                >  4:12606/1‑149 (MQ=255)
                                                                                                                                          aCTGCAT‑GGAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCtt        <  7:142525/148‑1 (MQ=255)
                                                                                                                                            tCCAT‑GCAGTTCACTGCGTTGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTAc     >  4:15897/1‑149 (MQ=255)
                                                                                                                                               aT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTACTGc  <  5:191251/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAT‑GCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTACTGC  >  NC_002947/5641314‑5641605

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