Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 751,788 G→T 21.3% intergenic (+295/‑552) ybgG → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606751,7880G→T21.3% 2.0 ‑2.2 0.031intergenic (+295/‑552)ybgG/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  new base (6/2):  ref base (9/14):  total (15/16)
Fisher's exact test strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATT  >  REL606/751757‑751823
                               |                                   
cTAAATGTTAATATTGGCGGGGTGGATTTATTCCtt                                 <  1:6367914/36‑1 (MQ=255)
 tAAATGTTAATATTGGCGGGGTGGATTTATGCCttt                                <  1:7015988/36‑1 (MQ=255)
  aaaTGTTAATATTGGCGGGGTGGGTTTTTTCCtttt                               >  1:2642377/1‑35 (MQ=21)
  aaaTGTTAATATTGGCGGGGTGGATTTATTCCTtta                               <  1:6240503/36‑1 (MQ=255)
   aaTGTTAATATTGGCGGGGTGGGTTTTTTCCttttt                              >  1:5638999/1‑34 (MQ=25)
   aaTGTTAATATTGGCGGGGTGGGTTTTTTCCtttct                              >  1:4904949/1‑34 (MQ=25)
        tAATATTGGCGGGGTGGATTTATGCCTTTATTATTa                         >  1:1313261/1‑36 (MQ=255)
        tAATATTGGCGGGGTGGAGTTATGCCTTTATTAGTa                         <  1:4765246/36‑1 (MQ=255)
        tAATATTGGCGGGGAGGCTTTATTTGTttattattt                         >  1:1573236/1‑33 (MQ=9)
         aaTATTGGCGGGGTGGATTTATTCCTTTTTTATTaa                        >  1:7387/1‑36 (MQ=21)
         aaTATTGGCGGGGTGGATTTATGCCTTTATTAGTaa                        <  1:6813817/36‑1 (MQ=255)
          atatTGGCGGGGTGGATTTATGCCTTTATTAGTAAt                       <  1:2928739/36‑1 (MQ=255)
          atatTGGCGGGGGGGGTTTTTTCCTTTTTTTTTAAt                       >  1:7270706/1‑36 (MQ=2)
           tatTGGCGGGGTGGATTTATGCCTTTATTAGTAATc                      <  1:1912715/36‑1 (MQ=255)
           tatTGGCGGGGTGGATTTATGCCTTTATTAGTAATc                      <  1:6633205/36‑1 (MQ=255)
             ttGGCGGGGTGGATTTATGCCTTTATTAGTAATCCt                    >  1:4086382/1‑36 (MQ=255)
             ttGGCGGGGTGGATTTATGCCTTTATTAGTAATCCt                    >  1:4677319/1‑36 (MQ=255)
               ggCGGGGTGGATTTATGCCTTTATTAGTAATCCTGa                  <  1:1798152/36‑1 (MQ=255)
                gCGGGGTGGATTTATGCCTTTATTAGTAATCCTGaa                 >  1:4153834/1‑36 (MQ=255)
                   gggTGGAATTATGCCTTTATTAGTAATCCTGAAAct              <  1:6120842/36‑1 (MQ=255)
                    ggTGGATTTATGCCTTTATTAGTAATCCTGAAActc             >  1:7210700/1‑36 (MQ=255)
                     gTGGATTTATGCCTTTATTAGTAATCCTGAAActct            <  1:6534644/36‑1 (MQ=255)
                       ggATTTATGCCTTTATTAGTAATCCTGAAACTCTGc          <  1:1998097/36‑1 (MQ=255)
                        gATTTATGCCTTTATTAGTAATCCTGAAACTCTGcg         >  1:378416/1‑36 (MQ=255)
                          tttATGCCTTTATTAGTAATCCTGAAACTCTGcgtc       <  1:165108/36‑1 (MQ=255)
                           ttATGCCTTTATTAGTAATCCTGAAACTCTGcgtcg      <  1:7240872/36‑1 (MQ=255)
                             aTGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTa    >  1:5148183/1‑36 (MQ=255)
                             aTGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTa    >  1:1819286/1‑36 (MQ=255)
                             aTGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTa    <  1:1080432/36‑1 (MQ=255)
                              tGCTTTTATTAGTAATCCTGAAACTCTGCGTCGTAt   <  1:7486542/36‑1 (MQ=255)
                              tGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTAt   >  1:2270539/1‑36 (MQ=255)
                               gCCTTTATTAGTAATCCTGAAACTCTGCGTCGTAtt  <  1:4543802/36‑1 (MQ=255)
                               |                                   
CTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATT  >  REL606/751757‑751823

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG