Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,704,256 +G 5.3% coding (1089/1185 nt) ygaX → predicted transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,704,2561.→G 5.3% 3.4 ‑3.2 ‑0.138coding (1089/1185 nt)ygaXpredicted transporter
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (6/30):  total (8/30)
Fisher's exact test strand distribution p-value = 3.98e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.46e-01

GGTTCGCTAATTTCAGCCTCAGCCTGGCAACATGGC.GGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGAT  >  REL606/2704221‑2704291
                                    |                                   
ggTTCGCTAATTTCAGCCTCAGCCTGGCAACATGGc                                      <  1:2757580/36‑1 (MQ=255)
 gTTCGCTAATTTCAGCCTCAGCCTGGCAACATGGC.g                                    <  1:5084044/36‑1 (MQ=255)
 gTTCGCTAATTTCAGCCTCAGCCTGGCAACATGGC.g                                    <  1:2991518/36‑1 (MQ=255)
   tCGCTAATTTCAGCCTCAGCCTGGCAACATGGC.GGt                                  <  1:1307544/36‑1 (MQ=255)
    cGCTAATTTCAGCCTCAGCCTGGCAACATGGC.GGtt                                 <  1:7727023/36‑1 (MQ=255)
    cGCTAATTTCAGCCTCAGCCTGGCAACATGGC.GGtt                                 <  1:161570/36‑1 (MQ=255)
    cGCTAATTTCAGCCTCAGCCTGGCAACATGGC.GGtt                                 <  1:7722634/36‑1 (MQ=255)
    cGCTAATTTCAGCCTCAGCCTGGCAACATGGC.GGtt                                 <  1:5710020/36‑1 (MQ=255)
    cGCTAATTTCAGCCTCAGCCTGGCAACATGGC.GGtg                                 >  1:1735845/1‑35 (MQ=255)
      cTAATTTCAGCCTCAGCCTGGCAACATGGC.GGTTgg                               <  1:5192455/36‑1 (MQ=255)
      cTAATTTCAGCCTCAGCCTGGCAACATGGC.GGGTgg                               >  1:4421048/1‑36 (MQ=255)
        aaTTTCAGCCTCAGCCTGGCAACATTGC.GGTTGGGc                             <  1:1845321/36‑1 (MQ=255)
        aaTTTCAGCCTCAGCCTGGCAACATGGC.GGTTGGGc                             <  1:5319751/36‑1 (MQ=255)
        aaTTTCAGCCTCAGCCTGGCAACATGGC.GGTTGGGc                             <  1:6884284/36‑1 (MQ=255)
         aTTTCAGCCTCAGCCTGGCAACATGGCGGGTTgggg                             >  1:2882683/1‑35 (MQ=37)
         aTTTCAGCCTCAGCCTGGCAACATGGC.GGGTggggt                            >  1:2846366/1‑34 (MQ=38)
          tttCAGCCTCAGCCTGGCAACATGGC.GGTTtggggg                           >  1:5414796/1‑33 (MQ=37)
          tttCAGCCTCAGCCTGGCAACATGGC.GGTTGGGCTg                           <  1:4602689/36‑1 (MQ=255)
          tttCAGCCTCAGCCTGGCAACATGGC.GGTTGGGCTg                           <  1:1027244/36‑1 (MQ=255)
           ttCAGCCTCAGCCTGGCAACATGGCGGGTTgggggg                           >  1:5389044/1‑33 (MQ=25)
           ttCAGCCTCAGCCTGGCAACATGGC.GGTTGGGCTgg                          <  1:4095932/36‑1 (MQ=255)
           ttCAGCCTCAGCCTGGCAACATGGC.GGTTGGGCTgg                          <  1:7794317/36‑1 (MQ=255)
           ttCAGCCTCAGCCTGGCAACATGGC.GGGTGGGGTgg                          >  1:3507948/1‑36 (MQ=37)
            tcagcctcagccTGGCAACATGGC.GGTTGGGCTGGc                         <  1:3991203/36‑1 (MQ=255)
            tcagcctcagccTGGCAACATGGC.GGTTGGGCTGGc                         <  1:3342051/36‑1 (MQ=255)
              agcctcagcctGGCAACATGGC.GGTTGGGCTGGCGt                       <  1:7809321/36‑1 (MQ=255)
                cctcagcctGGCAACATGGC.GGTTGGGCTGGCGttt                     <  1:5692745/36‑1 (MQ=255)
                  tcagcctGGCAACATGGC.GGGTGGGGGGGGGttttt                   >  1:789861/1‑34 (MQ=16)
                   cagcctGGCAACATGGC.GGTTGGGCTGGCGTTTGTc                  <  1:4814708/36‑1 (MQ=255)
                    agcctGGCAACATGGC.GGTTGGGCTGGCGTTTGTct                 <  1:6543763/36‑1 (MQ=255)
                    agcctGGCAACATGGC.GGTTGGGCTGGCGTTTGTct                 <  1:20168/36‑1 (MQ=255)
                     gcctGGCAACATGGC.GGTTGGGCTGGCGTTTGTctg                <  1:64509/36‑1 (MQ=255)
                       ctGGCAACATGGC.GGTTGGGCTGGCGTTTGTctggc              <  1:3348340/36‑1 (MQ=255)
                        tGGCAACATGGC.GGTTGGGCTGGCGTTTGTctggct             <  1:2761747/36‑1 (MQ=255)
                         ggCAACATGGC.GGTTGGGCTGGCGTTTGTctggctg            <  1:281043/36‑1 (MQ=255)
                          gCAACATGGC.GGTTGGGCTGGCGTTTGTctggctgg           <  1:7004337/36‑1 (MQ=255)
                           cAACATGGC.GGTTGGGCTGGCGTTTGTCtggttggc          <  1:1721840/36‑1 (MQ=255)
                           cAACATGGC.GGTTGGGCTGGCGTTTGTCtggctggc          <  1:7124666/36‑1 (MQ=255)
                             aCATGGC.GGTTGGGCTGGCGTTTGTCTGGCTGgcgc        <  1:4234615/36‑1 (MQ=255)
                               aTGGC.GGTTAGGCTGGCGTTTGTCTGGCTGGCGcga      <  1:6336005/36‑1 (MQ=255)
                                tGGC.GGTTGGGCTGGCGTTTGTCTGGCTGGCGcgac     <  1:6488293/36‑1 (MQ=255)
                                   c.GGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGAt  <  1:4969124/36‑1 (MQ=255)
                                   c.GGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGAt  >  1:3697175/1‑36 (MQ=255)
                                    |                                   
GGTTCGCTAATTTCAGCCTCAGCCTGGCAACATGGC.GGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGAT  >  REL606/2704221‑2704291

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG