Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 987,461 G→T 30.6% L327L (CTG→CTT ycaQ → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606987,4610G→T30.6% 7.3 ‑2.8 0.038981ycaQhypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/1):  ref base (15/10):  total (27/11)
Fisher's exact test strand distribution p-value = 5.97e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCCGAAACGCCAGTATGGCTATTTTGTTCTGCCGTTATTACATCGTGGGCAATTAGTTGGGCGAAT  >  REL606/987431‑987496
                              |                                   
gcCGAAACGGCAGTATGGCNATTTTTGTCTTCCggt                                >  1:2161717/1‑34 (MQ=12)
gcCGAAACGCCAGTATGGCTTTTTTGTTCTGCCGtt                                <  1:7305731/36‑1 (MQ=255)
 ccGAAACGCCAGTATGGCTATTTTTTTTTTCCTtta                               >  1:2791993/1‑36 (MQ=16)
 ccGAAACGCCAGTATGGCTATTTTTTTCTTCCGtta                               >  1:4498557/1‑36 (MQ=37)
 ccGAAACGCCAGTATGGCTATTTTGTTCTGCCGtta                               <  1:5388104/36‑1 (MQ=255)
  cGAAACGCCAGTATGGCTATTTTGTTCTGCCGtttt                              >  1:1741199/1‑34 (MQ=255)
    aaaCGCCAGTATGGCTATTTTTTTCTTCCGtttttt                            >  1:3046655/1‑35 (MQ=21)
    aaaCGCCAGTATGGCTATTTTGTTCTTCCGttatta                            >  1:1250938/1‑36 (MQ=255)
    aaaCGCCAGTATGGCTATTTTGTTCTTCCGttatta                            >  1:6830399/1‑36 (MQ=255)
    aaaCGCCAGTATGGCTATGTTGTTGTGCCGttgtta                            >  1:5723134/1‑36 (MQ=21)
     aaCGCCAGTATGGCTATTTTTTTCTTCCGTTTTTAc                           >  1:3117652/1‑36 (MQ=21)
     aaCGCCAGTATGGCTATTTTTTTCTTCCGTTATTAc                           >  1:5649252/1‑36 (MQ=37)
     aaCGCCAGTATGGCTATTTTTTTCTTCCGTTATTAc                           >  1:4925116/1‑36 (MQ=37)
      aCGCCAGTATGGCTATTTTTTTCTTCCGttattttc                          >  1:1759794/1‑33 (MQ=21)
      aCGCCAGTATGGCTATTTTTTTCTGCCGttattata                          >  1:1038006/1‑34 (MQ=37)
       cGCCAGTATGGCTATTTTGTTCTGCCGTTATTACAt                         <  1:3502621/36‑1 (MQ=255)
         ccAGTATGGCTATTTTTTTCTTCCGTTATTACATCg                       >  1:5948342/1‑36 (MQ=37)
         ccAGTATGGCTATTTTGTTCTGCCGTTATTACATCg                       <  1:1381277/36‑1 (MQ=255)
          cAGTATGGCTATTTTGTTCTGCCGTTATTACATCGt                      <  1:2850584/36‑1 (MQ=255)
          cAGGATGGCTATTTTGTCCTTCCGTTATTACATCGt                      <  1:830915/36‑1 (MQ=37)
           aGTATGGCTATTTTGTTCTGCCGTTATTACATCGTg                     <  1:1592038/36‑1 (MQ=255)
             tATGGCTATTTTGTTCTGCCGTTATTACATCGTggg                   >  1:6490146/1‑36 (MQ=255)
             tATGGCTATTTTGTTCTGCCGTTATTACATCGTggg                   <  1:2241650/36‑1 (MQ=255)
              aTGGCTATTTTGTTTTTCCGTTATTACATCGTGGGc                  >  1:7459390/1‑36 (MQ=37)
                ggCTATTTTGTTCTGCCGTTATTACATCGTGGGCaa                >  1:5682376/1‑36 (MQ=255)
                 gCTCTTTTGTTCTGCCGTTATTACATCGTGGGCAAt               <  1:2274132/36‑1 (MQ=255)
                 gCTATTTTGTTCTGCCGTTATTACATCGTGGGCAAt               >  1:3135034/1‑36 (MQ=255)
                   tATTTTGTTCTGCCGTTATTACATCGTGGGCAATTa             >  1:5938201/1‑36 (MQ=255)
                   tATTTTGTTCTGCCGTTATTACATCGTGGGCAATTa             >  1:4247696/1‑36 (MQ=255)
                    aTTTTGTTCTGCCGTTATTACATCGTGGGCAATTAg            >  1:778279/1‑36 (MQ=255)
                     ttttGTTCTGCCGTTATTACATCGTGGGCAATTAGt           <  1:6999754/36‑1 (MQ=255)
                      tttGTTCTGCCGTTATTACATCGTGGGCAATTAGtt          <  1:2857208/36‑1 (MQ=255)
                       ttGTTCTGCCGTTATTACATCGTGGGCAGTTAGttt         >  1:4347796/1‑35 (MQ=255)
                        tGTTCTGCCGTTATTACATCGTGGGCAATTAGTTgg        >  1:7039285/1‑36 (MQ=255)
                        tGTTCTGCCGTTATTACATAGTGGGCAATTAGTTgg        >  1:6337823/1‑36 (MQ=255)
                          ttCTGCCGTTATTACATCGTGGGCAATTAGTTGGGc      >  1:7339579/1‑36 (MQ=255)
                           tCTGCCGTTATTACATCGTGGGCAATTAGTTGGGCg     >  1:3390682/1‑36 (MQ=255)
                             tGCCGTTATTACATCGTGGGCAATTAGTTGGGCGaa   >  1:7506000/1‑36 (MQ=255)
                              gCCGTTATTACATCGTGGGCAATTAGTTGGGCGAAt  <  1:900725/36‑1 (MQ=255)
                              gCCGTTATTACATCGTGGGCAATTAGTTGGGCGAAt  >  1:98435/1‑36 (MQ=255)
                              |                                   
GCCGAAACGCCAGTATGGCTATTTTGTTCTGCCGTTATTACATCGTGGGCAATTAGTTGGGCGAAT  >  REL606/987431‑987496

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