Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 715,754 A→T 19.0% R121* (AGA→TGA)  ybfC → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606715,7540A→T19.0% 5.8 ‑2.6 0.030361ybfChypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (9/15):  total (14/16)
Fisher's exact test strand distribution p-value = 7.25e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAGCAGAAAAAACGGCGTTTCTCGGCGATCGGTTTAGAAGTTATATAAGAGCGGGTATCTATCAAAAAAAG  >  REL606/715719‑715789
                                   |                                   
cagcagAAAAAACGGCGTTTCTCGGCGATCGGTTTa                                     <  1:982405/36‑1 (MQ=255)
cagcagAAAAAACGGCGTTTCTCGGCGATCGGTTTa                                     <  1:2271209/36‑1 (MQ=255)
 agcagAAAAAACGGCGTTTCTCGGCGATCGGTTTAg                                    >  1:1700855/1‑36 (MQ=255)
 agcagAAAAAACGGCGTTTCTCGGCGATCGGGTTAg                                    >  1:677239/1‑36 (MQ=255)
 agcagAAAAAACGGCGTTTCTCGCCGATCGGTTTAg                                    >  1:922056/1‑36 (MQ=255)
  gcagAAAAAACGGCGTTTCTCGGCGATCGGGTTTGa                                   >  1:2684689/1‑36 (MQ=37)
   cagAAAAAACGGCGTTTCTCGGCGATCGGTTTAGaa                                  <  1:1769746/36‑1 (MQ=255)
   cagAAAAAACGGCGTTTCTCGGCGATCGGTTTAGaa                                  >  1:885266/1‑36 (MQ=255)
    agAAAAAACGGCGTTTCTCGGTTATCGGTTTAGatg                                 >  1:5871495/1‑34 (MQ=25)
    agAAAAAACGGCGTTTCTCGGCGATCGGTTTAGAAg                                 >  1:3627464/1‑36 (MQ=255)
      aaaaagCGGCGTTTCTCGGCGATCGGTTTAGAAGtt                               <  1:4831264/36‑1 (MQ=255)
      aaaaaaCGGCGTTTCTCGGCGATCGGTTTAGAAGtt                               <  1:4646215/36‑1 (MQ=255)
       aaaaaCGGCGTTTCTCGGCGATCGGTTTAGAAGTta                              <  1:1409200/36‑1 (MQ=255)
         aaaCGGCGTTTTTCGGCGATCGGTTTAGAAGTtata                            <  1:5522772/36‑1 (MQ=255)
         aaaCGGCGTTTCTCGGCGATCGGTTTTGAAGTtata                            >  1:2320526/1‑36 (MQ=255)
         aaaCGGCGTTTCTCGGCGATCGGTTTAGAAGTtata                            <  1:3204224/36‑1 (MQ=255)
         aaaCGGCGTTTCTCGGCGATCGGGTTAGAAGTtata                            >  1:5276503/1‑36 (MQ=255)
           aCGGCGTTTCTCGGCGATCGGTTTTGAAGTtatata                          >  1:2387622/1‑36 (MQ=255)
             ggCGTTTCTCGGCGATCGGGTTTTAAGTTATATAag                        >  1:3728322/1‑36 (MQ=25)
               cGTTTCTCGGCGATCGGGTTTGAAGTTATATAagag                      >  1:3080539/1‑36 (MQ=37)
                    ctcGGCGATCGGTTTTGAAGTTATATAAGAGCGGGt                 <  1:5545693/36‑1 (MQ=255)
                    ctcGGCGATCGGTTTAGAAGTTATATAAGAGCGGGt                 <  1:3939484/36‑1 (MQ=255)
                    ctcGGCGATCGGTTTAGAAGTTATATAAGAGCGGGt                 <  1:2198026/36‑1 (MQ=255)
                    ctcGGCGATCGGTTTAGAAGTTATATAAGAGCGGGt                 <  1:1319688/36‑1 (MQ=255)
                     tcGGCGATCGGTTTAGAAGTTATATAAGAGCGGGta                <  1:3269019/36‑1 (MQ=255)
                       ggCGATCGGTTTAGAAGTTATATAAGAGCGGGtatc              <  1:4280708/36‑1 (MQ=255)
                       ggCGATCGGCTTAGAAGTTATATAAGAGCGGGtatc              >  1:3810546/1‑36 (MQ=255)
                            tCGGTTTAGAAGTTATATAAGAGCGGGTATCTATCa         <  1:1913313/36‑1 (MQ=255)
                             cGGTTTAGAAGTTATATAAGAGCGGGTATCTATCaa        >  1:2676264/1‑36 (MQ=255)
                               gTTTAGAAGTTATATAAGAGCGGGTATCTATCaaaa      <  1:369641/36‑1 (MQ=255)
                                  tAGAAGTTATATAAGAGCGGGTATCTATCaaaaaaa   <  1:4372805/36‑1 (MQ=255)
                                  tAGAAGTTATATAAGAGCGGGTATCTATCaaaaaaa   <  1:3262963/36‑1 (MQ=255)
                                   aGAAGTTATATAAGAGCGGGTATCTATCAAAAAAag  >  1:5701170/1‑36 (MQ=255)
                                   aGAAGTTATATAAGAGCGGGTATCTATCAAAAAAag  >  1:2360601/1‑36 (MQ=255)
                                   |                                   
CAGCAGAAAAAACGGCGTTTCTCGGCGATCGGTTTAGAAGTTATATAAGAGCGGGTATCTATCAAAAAAAG  >  REL606/715719‑715789

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