Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,280,890 +G 5.3% coding (1011/2052 nt) yjcQ ← predicted multidrug efflux system component

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,280,8901.→G 5.3% 3.5 ‑1.8 ‑0.038coding (1011/2052 nt)yjcQpredicted multidrug efflux system component
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (14/22):  total (16/22)
Fisher's exact test strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.30e-01

TCTGGATTGGTAAAAGCATCGGCGGCCATTGATGGC.GGTTTGGCTGCGGGCGTTGGCGGCGTATTCGGGT  >  REL606/4280855‑4280924
                                    |                                  
tCTGGATTGGTAAAAGCATCGGCGGCCATTGATGGc                                     <  1:1935658/36‑1 (MQ=255)
tCTGGATTGGTAAAAGCATCGGCGGCCATTGATGGc                                     >  1:2752696/1‑36 (MQ=255)
 cTGGATTGGTAAAAGCATCGGCGGCCATTGATGGC.g                                   >  1:3618625/1‑36 (MQ=255)
 cTGGATTGGTAAAAGCATCGGCGGCCATTGATGGC.g                                   >  1:2254262/1‑36 (MQ=255)
 cTGGATTGGTAAAAGCATCGGCGGCCATTGATGGC.g                                   >  1:2697224/1‑36 (MQ=255)
  tGGATTGGTAAAAGCATCGGCGGCCATTGATGGC.gg                                  <  1:6906391/36‑1 (MQ=255)
   ggATTGGTAAAAGCATCGGCGGCCATTGATGGC.GGt                                 <  1:7147252/36‑1 (MQ=255)
   ggATTGGTAAAAGCATCGGCGGCCATTGATGGC.GGt                                 <  1:6376748/36‑1 (MQ=255)
    gATTGGTAAAAGCATCGGCGGCCATTGATGGC.gggt                                >  1:1135459/1‑34 (MQ=255)
    gATTGGTAAAAGCATCGGCGGCCATTGATGGC.GGtt                                <  1:7218097/36‑1 (MQ=255)
     aTTGGTAAAAGCATCGGCGGCCATTGCTGGC.GTttt                               >  1:5341785/1‑36 (MQ=37)
     aTTGGTAAAAGCATCGGCGGCCATTGATGGC.GGttt                               <  1:4853816/36‑1 (MQ=255)
     aTTGGTAAAAGCATCGGCGGCCATTGATGGC.GGttt                               >  1:3128552/1‑36 (MQ=255)
     aTTGGTAAAAGCATCGGCGGCCATTGATGGC.GGttt                               >  1:96827/1‑36 (MQ=255)
      ttGGTAAAAGCATCGGCGGCCATTGATGGC.GGTTTg                              <  1:5162893/36‑1 (MQ=255)
        ggTAAAAGCATCGGCGGCCCTTGATGGC.GGttttgg                            >  1:2863192/1‑33 (MQ=38)
        ggTAAAAGCATCGGCGGCCATTGATGGC.GGttttgg                            >  1:4403509/1‑33 (MQ=255)
        ggTAAAAGCATCGGCGGCCATTGATGGC.GGTTTGGc                            <  1:6240368/36‑1 (MQ=255)
        ggTAAAAGCATCGGCGGCCATTGATGGC.GGTTTGGc                            <  1:4056296/36‑1 (MQ=255)
         gTAAAAGCATCGGCGGCCCTTTATGGCGGGTTTggt                            >  1:3737242/1‑35 (MQ=16)
         gTAAAAGCATCGGCGGCCATTGATGGC.GGTTTGGCt                           <  1:3074673/36‑1 (MQ=255)
         atAAAAGCATCGGCGGCCATTGATGGC.GGTTTGGCt                           <  1:5742347/35‑1 (MQ=255)
          tAAAAGCATCGGCGGCCATTGATGGC.GGTTTGGCTg                          <  1:1097278/36‑1 (MQ=255)
           aaaaGCATCGGCGGCCATTGATGGC.GGTTTGGCTGc                         <  1:4237604/36‑1 (MQ=255)
           aaaaGCATCGGCGGCCATTGATGGC.GGGTTGGCTgg                         >  1:6700219/1‑35 (MQ=255)
            aaaGCATCGGCGGTCATTGATGGCGGGTTTGGTTgg                         >  1:4302239/1‑35 (MQ=16)
            aaaGCATCGGCGGCCATTGATGGC.GGGTTTGCTggg                        >  1:1923971/1‑34 (MQ=25)
             aaGCATCGGCGGCCATTGATGGC.GGTTTGGCTGCgg                       <  1:7432097/36‑1 (MQ=255)
              aGCATCGGCGGCCATTGATGGC.GGTTTTGGTGCggg                      >  1:6756413/1‑36 (MQ=37)
              aGCATCGGCGGCCATTGATGGC.GGTTTGGCTGCggg                      <  1:3039610/36‑1 (MQ=255)
               gCATCGGCGGCCATTGATGGC.GGTTTGGCTGCGGGc                     <  1:651665/36‑1 (MQ=255)
                cATCGGCGGCCATTGATGGC.GGTTTGGCTGCGGGCg                    <  1:538691/36‑1 (MQ=255)
                 aTCGGCGGCCATTGATGTC.GGTTTGGCTGCGGGCGt                   <  1:2209530/36‑1 (MQ=255)
                   cggcggCCATTGATGGC.GGTTTTGCTGGGGGCGTTg                 >  1:5952108/1‑36 (MQ=25)
                     gcggcCATTGATGGC.GGTTTTGCTGCGGGCGtttgg               >  1:2047455/1‑33 (MQ=37)
                       ggcCATTGATGGC.GGTTTGGCTGCGGGGGTTggcgg             >  1:5002300/1‑36 (MQ=255)
                           aTTGATGGC.GGTTTGGCTGCGGGCGTTTGGGGCGTa         >  1:40446/1‑36 (MQ=37)
                           aTTGATGGC.GGTTTGGCTGCGGGCGTTGGCGGCGTa         <  1:3389124/36‑1 (MQ=255)
                             tGATGGC.GGTTTGGCTGCGGGCGTTGGCGGCGTAtt       <  1:2919174/36‑1 (MQ=255)
                              gATGGC.GGTTTGTCTGCGGGCGTTGGCGGCGTATTc      <  1:1643044/36‑1 (MQ=255)
                              gATGGC.GGTTTGGCTGCGGGCGTTGGCGGCGTATTc      <  1:170864/36‑1 (MQ=255)
                               agggC.GGTTTGGCTGCGGGCGTTGGCGGCGTATTCg     <  1:3449127/34‑1 (MQ=255)
                                tGGC.GGTTTGGCTGCGGGCGTTGGCGGCGTATTCgg    <  1:3319111/36‑1 (MQ=255)
                                 ggC.GGTTTGGCTGCGGGCGTTGGCGGCGGATTCggg   >  1:3822538/1‑36 (MQ=255)
                                  gC.GGTTTGGCTGCGGGCGTTGGCGGCGTATTCgggg  >  1:2526040/1‑35 (MQ=255)
                                    |                                  
TCTGGATTGGTAAAAGCATCGGCGGCCATTGATGGC.GGTTTGGCTGCGGGCGTTGGCGGCGTATTCGGGT  >  REL606/4280855‑4280924

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG