Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,279,523 A→G 19.7% L527P (CTC→CCC)  yfaS ← putative membrane protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,279,5230A→G19.7% 7.3 ‑2.5 ‑0.0341580yfaSputative membrane protein
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (17/10):  total (24/10)
Fisher's exact test strand distribution p-value = 7.82e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.80e-01

AACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGTCACATGACCGCTGGCGTCGGACACC  >  REL606/2279490‑2279558
                                 |                                   
aaCAAATAGCGGCTCGGTTTATCGGCGGCGGGgaga                                   <  1:3933482/36‑1 (MQ=255)
aaCAAATAGCGGCTCGGTTTATCGGCGGCGGGgaga                                   <  1:3360901/36‑1 (MQ=255)
  cAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGAtt                                 <  1:7605641/36‑1 (MQ=255)
   aaaTAGCGGCTCGGTTTATCGGCGGGGGGGGGGttt                                >  1:5950260/1‑36 (MQ=25)
   aaaTAGCGGCTCGGTTTATCGGCGGGGGGGGGGttt                                >  1:1935547/1‑36 (MQ=25)
   aaaTAGCGGCTCGGTTTATCGGCGGGGGGGGGGttt                                >  1:2792341/1‑36 (MQ=25)
        gCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGt                           >  1:6256665/1‑36 (MQ=37)
        gCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGt                           <  1:262214/36‑1 (MQ=255)
         cGGCTCGGTTTATCGGCGGCGGGGGGGTTTAAAGtt                          >  1:7559982/1‑35 (MQ=21)
          ggCTCGGTTTATCGGCGGCGGGGGGGTTTAAAGTca                         >  1:4691344/1‑36 (MQ=25)
           gCTCGGTTTATCGGCGGCGGGGAGATTTAACGTcac                        <  1:7193246/36‑1 (MQ=255)
            cTCGGTTTATCGGCGGCGGGGAGATTTAACGTcaca                       >  1:818239/1‑36 (MQ=255)
             tCGGTTTATCGGCGGCGGGGAGATTTAACGTCACAt                      <  1:5278445/36‑1 (MQ=255)
              cGGTTTATCGGCGGCGTGGAGATTTAACGTCACATg                     >  1:2510261/1‑36 (MQ=255)
              cGGTTTATCGGCGGCGGGGAGGTTTTAAGTCACATg                     >  1:3273064/1‑36 (MQ=25)
              cGGTTTATCGGCGGCGGGGAGATTTAACGTCACATg                     >  1:1480255/1‑36 (MQ=255)
                gTTTATCGGCGGAGGGGAGATTTAACGTGACAGGat                   >  1:3664369/1‑35 (MQ=25)
                 tttATCGGCGGCGGGGAGATTTAACGTCACATGAcc                  <  1:2979066/36‑1 (MQ=255)
                 tttATCGGCGGCGGGGAGATTTAACGTCACATGAcc                  <  1:3936106/36‑1 (MQ=255)
                  ttATCGGCGGCGGGGAGATTTAACGTCACATGACCg                 <  1:4128408/36‑1 (MQ=255)
                  ttATCGGCGGCGGGGAGATTTAACGTCACATGACCg                 >  1:7640990/1‑36 (MQ=255)
                  ttATCGGCGGCGGGGAGATTTAACGTCAAATGACCg                 >  1:2901012/1‑36 (MQ=255)
                   tATCGGCGGCGGGGGGATTTAAAGTGACATGACCGc                >  1:7124323/1‑36 (MQ=25)
                    aTCGGCGGCGGGGAGATTTAACGTCACATGACCGCt               <  1:1406993/36‑1 (MQ=255)
                     tCGGCGGCGGGGAGATTTAACGTCACATGACCGCTg              >  1:4356761/1‑36 (MQ=255)
                      cggcggcggGGAGATTTAACGTCACATGACCGCTgg             >  1:8043623/1‑36 (MQ=255)
                       ggcggcggGGAGATTTAACGTCACATGACCGCTGGc            >  1:5749566/1‑36 (MQ=255)
                       ggcggcggGGAGATTTAACGTCACATGACCGCTGGc            <  1:4543977/36‑1 (MQ=255)
                        gcggcggGGAGATTTAACGTCACATGACCGCTGGCg           >  1:6767448/1‑36 (MQ=255)
                         cggcggGGAGATTTAACGTCACATGACCGCTGGCGt          <  1:508627/36‑1 (MQ=255)
                           gcggGGAGATTTAACGTCACATGACCGCTGGCGTCg        >  1:1529092/1‑36 (MQ=255)
                            cggGGAGATTTAACGTCACATGACCGCTGGCGTCgg       >  1:4822926/1‑36 (MQ=255)
                             ggggAGATTTAACGTCACATGACCGCTGGCGTCGGa      >  1:1173595/1‑36 (MQ=255)
                             ggggAGATTTAACGTCACATGACCGCTGGCGTCGGa      >  1:7587298/1‑36 (MQ=255)
                             ggggAGATTTAACGTCACATGAACGCTGGCGTCGGa      >  1:3960479/1‑36 (MQ=255)
                               ggAGATTTAACGTCACATGACCGCTGGAGTCGGaca    >  1:6842200/1‑36 (MQ=255)
                                gagaTTTAACGTCACATGACCGCTGGCGTCGGacac   >  1:2908452/1‑36 (MQ=255)
                                gagaTTTAACGTCACATGACCGCTGGCGTCGGacac   >  1:2533892/1‑36 (MQ=255)
                                gagaTTTAACGTCACATGACCGCTGGCGTCGGacac   <  1:1293121/36‑1 (MQ=255)
                                gagaTTTAACGTCACATGACCGCTGGCGTCGGacac   >  1:6396751/1‑36 (MQ=255)
                                 agaTTTAACGTCACATGACCGCTGGCGTCGGACAcc  >  1:4680841/1‑36 (MQ=255)
                                 |                                   
AACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGTCACATGACCGCTGGCGTCGGACACC  >  REL606/2279490‑2279558

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG