Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,641,455 C→G 9.8% S373C (TCC→TGC)  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,4550CG 9.8% 120.6 / 4.9 41S373C (TCC→TGC) PP_4951EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (20/17);  new base G (2/2);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

CCGCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTA  >  NC_002947/5641307‑5641601
                                                                                                                                                    |                                                                                                                                                  
ccgcgtaTCGAACGGTTCGGCTATCGGATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTc                                                                                                                                                    <  6:119241/143‑1 (MQ=255)
  gCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCa                                                                                                                                                  >  6:17000/1‑149 (MQ=255)
      aTCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATg                                                                                                                                                <  5:139124/147‑1 (MQ=255)
      aTCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATg                                                                                                                                                <  8:340488/147‑1 (MQ=255)
      aTCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGCGCGCAGTGAACTCCATg                                                                                                                                                >  1:150734/1‑147 (MQ=255)
          aaCGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGGCTTCTGGGCGCAGTGACCTCCATGCAGtt                                                                                                                                           >  3:305624/1‑148 (MQ=255)
             ggTTCGGCTATCGGATCATCCTCACCGGTAATACCTTGCTGTTGGGCTTCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGTTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCATGGAGTTCACt                                                                                                                                       <  7:121719/149‑1 (MQ=255)
              gTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCATGGAGTTCACt                                                                                                                                       <  3:231512/148‑1 (MQ=255)
              gTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCATGGAGTTCACt                                                                                                                                       >  4:231512/1‑148 (MQ=255)
                       tCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGaa                                                                                                                               <  2:36133/147‑1 (MQ=255)
                                  tCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAt                                                                                                                                                 <  8:281318/118‑1 (MQ=255)
                                  tCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAt                                                                                                                                                 >  7:281318/1‑118 (MQ=255)
                                   cACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCa                                                                                                                 >  1:352189/1‑149 (MQ=255)
                                           aTACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGAtctc                                                                                                         <  3:170680/149‑1 (MQ=255)
                                                                  ctgGCAAGCCTGGGCCTGGTGGACGAGCGGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGaa                                                                                                                               <  2:226745/104‑1 (MQ=255)
                                                                  ctgGCAAGCCTGGGCCTGGTGGACGAGCGGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGaa                                                                                                                               >  1:226745/1‑104 (MQ=255)
                                                                  ctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCgg                                                                                  >  8:11943/1‑149 (MQ=255)
                                                                  ctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCACCGAcagcaggct                                                                                  >  1:147877/1‑146 (MQ=255)
                                                                   tgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCgg                                                                                  >  2:233377/1‑148 (MQ=255)
                                                                     gCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGt                                                                                 <  2:352189/147‑1 (MQ=255)
                                                                       aaGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTaa                                                                               >  4:317405/1‑147 (MQ=255)
                                                                               gcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCctgctg                                                                     >  5:218509/1‑149 (MQ=255)
                                                                                  ggttcgACGAGCAGAACCCTTATTGGCTGCTGTTGGGGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACCCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCg                                                                  <  4:305624/144‑1 (MQ=255)
                                                                                          gAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTcgc                                                          >  6:83706/1‑149 (MQ=255)
                                                                                              aGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGt                                                                                 <  5:208526/122‑1 (MQ=255)
                                                                                                  cccTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGt                                                  >  3:70940/1‑149 (MQ=255)
                                                                                                   ccTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTc                                                 <  7:11943/149‑1 (MQ=255)
                                                                                                        ttGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGa                                            <  5:17000/149‑1 (MQ=255)
                                                                                                          ggCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCgcgc                                                        >  5:103158/1‑135 (MQ=255)
                                                                                                          ggCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCgcgc                                                        <  6:103158/135‑1 (MQ=255)
                                                                                                           gctgctGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGcc                                         >  8:338535/1‑149 (MQ=255)
                                                                                                              gctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTgg                                      <  4:70940/149‑1 (MQ=255)
                                                                                                               ctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTgg                                      <  7:338535/148‑1 (MQ=255)
                                                                                                                     gTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACCGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGc                               >  8:202977/1‑149 (MQ=255)
                                                                                                                                  ttCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGt                                                                                 >  4:120047/1‑86 (MQ=255)
                                                                                                                                  ttCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGt                                                                                 <  3:120047/86‑1 (MQ=255)
                                                                                                                                       ggCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGtt             >  3:149179/1‑149 (MQ=255)
                                                                                                                                         cgcAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTgg           >  8:116989/1‑149 (MQ=255)
                                                                                                                                            aGTGAACTGCATGGAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGggcg        >  5:306567/1‑149 (MQ=255)
                                                                                                                                            aGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGggcg        >  5:303010/1‑149 (MQ=255)
                                                                                                                                              tGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGggcggc      <  5:83706/149‑1 (MQ=255)
                                                                                                                                                  cTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTa  <  4:149179/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTA  >  NC_002947/5641307‑5641601

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