Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 876,898 C→G 23.5% A225G (GCG→GGG)  cmr → multidrug efflux system protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606876,8980C→G23.5% 2.9 ‑2.3 ‑0.117674cmrmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  new base (3/1):  ref base (3/10):  total (6/11)
Fisher's exact test strand distribution p-value = 9.87e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.65e-01

GCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGG  >  REL606/876863‑876933
                                   |                                   
gCTAGTGCTTAAGAACGGCCGCTTTGTGGCGGGggc                                     <  1:1362239/36‑1 (MQ=37)
  tGGTGCTGAAGACCGGCCGCTTTGTGGCGGGggcgc                                   <  1:6434282/36‑1 (MQ=255)
  tGGTGCTGAAGAACGGCCGCTTTGTGGCGGGggcgc                                   <  1:7610876/36‑1 (MQ=255)
  gggTGCTGAAGAACGGCCGCTTTGTGGCGGGggcgc                                   <  1:2934162/35‑1 (MQ=255)
   ggTGCTGAAGAACGGCCGCTTTGTGGCGGGggcgct                                  <  1:247450/36‑1 (MQ=255)
     tGCTGAAGAACGGCCGCTTTGTGGCGGGggcgctgg                                <  1:6287095/36‑1 (MQ=255)
     tGCTGAAGAACGGCCGCTTTGTGGCGGGggcgctgg                                <  1:4130012/36‑1 (MQ=255)
      gCTGAAGAACGGCGGCTTTGTGGCGGGggcgctggc                               <  1:7609655/36‑1 (MQ=255)
         gaagaaCGGCCGCTTTGTGGCGGGggcgctggcgct                            <  1:940639/36‑1 (MQ=255)
         gaagaaCGGCCGCTTTGTGGCGGGggcgctggcgct                            <  1:7149862/36‑1 (MQ=255)
           agaaCGGCCGCTTTGTGGCGGGGGggcgggggcggg                          >  1:264217/1‑36 (MQ=19)
            gaaCGGCCGCTTTGTGGCGGGGGCGCTGGCGCTggg                         <  1:2850743/36‑1 (MQ=255)
             aaCGGCCGCTTTGTGGCGGGGGGGTTGTTGCTgggg                        >  1:6136197/1‑35 (MQ=12)
              aCGGCCGCTTTGTGGCGGGGGGGGTGGCGGTggggg                       >  1:3574705/1‑34 (MQ=21)
                ggCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTc                     <  1:5823950/36‑1 (MQ=255)
                 gCCGCTTTGTGGCGGGGGCGCTTGGGCTTGGGtttg                    >  1:1703938/1‑34 (MQ=16)
                 gCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCg                    <  1:5681225/36‑1 (MQ=255)
                  ccGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGt                   <  1:3292755/36‑1 (MQ=255)
                  ccGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGt                   <  1:1744389/36‑1 (MQ=255)
                     cTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAg                <  1:2488428/36‑1 (MQ=255)
                      tttGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGt               <  1:4586683/36‑1 (MQ=255)
                       ttGTGGCGGGGGGGGTGGCGGTGGGGTTTGTTTGtt              >  1:6439308/1‑35 (MQ=2)
                       ttGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTc              <  1:3876202/36‑1 (MQ=255)
                          tgGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGc           <  1:5463690/36‑1 (MQ=255)
                          tgGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGc           >  1:4169678/1‑36 (MQ=255)
                          gggCGGGGGGCCTTTCGCTGGGATTCGTTAGTCTGc           <  1:7683928/35‑1 (MQ=19)
                           ggCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGcc          >  1:442966/1‑36 (MQ=255)
                              gggggCGCTGGCGCTGGGATTCGTTAGTCTGCCGtt       <  1:8328807/36‑1 (MQ=255)
                                 ggcgctggcgctGGGATTCGTTAGTCTGCCGTTGCt    <  1:1579143/36‑1 (MQ=255)
                                 ggcgctggcgctGGGATTCGTTAGTCTGCCGTTGCt    <  1:4786031/36‑1 (MQ=255)
                                  gcgctggcgctgGGATTCGTTAGTCTGCCGTTGCTg   >  1:8182930/1‑36 (MQ=255)
                                   cgctggcgctggGATTCGTTAGTCTTCCGTTGCTgt  >  1:6351308/1‑35 (MQ=255)
                                   cgctggcgctggGATTCGTTAGTCTGCCGTTGCTgg  >  1:718930/1‑36 (MQ=255)
                                   cgctggcgctggGATTCGTTAGTCTGCCGTTGCTgg  <  1:7141372/36‑1 (MQ=255)
                                   |                                   
GCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGG  >  REL606/876863‑876933

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG