Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,789,554 A→T 17.1% T155S (ACC→TCC)  yqcE → predicted transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,789,5540A→T17.1% 2.0 ‑4.3 ‑0.031463yqcEpredicted transporter
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (8/15):  total (14/17)
Fisher's exact test strand distribution p-value = 1.31e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.40e-01

CGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATG  >  REL606/2789521‑2789587
                                 |                                 
cgcgGCGTCGGTGTAATGTCGCTGGCGGTGTTTAcc                                 <  1:5679039/36‑1 (MQ=255)
  cggcgTCGGTGTCATGTCGCTGGCGGTGTTTACCAt                               <  1:3798211/36‑1 (MQ=255)
  cggcgTCGGTGTAATGTCGCTGGCGGTGTTTACCAt                               <  1:3317741/36‑1 (MQ=255)
  cggcgTCGGTGTAATGTCGCTGGCGGTGTTTACCAt                               <  1:2388760/36‑1 (MQ=255)
   ggcgTCGGTGTAATGTCGCTGGCGGGGTTTATCAtt                              >  1:5302660/1‑35 (MQ=37)
     cgTCGGTGTAATGTCGCTGGTGGTGTTTACCAtgtg                            <  1:2165820/36‑1 (MQ=255)
      gTCGGTGTAATGTCGCTGGCGGTGTTTTCCCTTTgg                           >  1:5378261/1‑36 (MQ=25)
      gTCGGTGTAATGTCGCTGGCGGTGTTTTCCATGTgg                           >  1:2299823/1‑36 (MQ=255)
      gTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTgg                           <  1:5420371/36‑1 (MQ=255)
        cGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGt                         <  1:845824/36‑1 (MQ=255)
         ggTGTAATGTCGCTGGCGGTGTTTACCATGTTGGTc                        >  1:5094038/1‑36 (MQ=255)
          gtgtGACGTCGCTGGCGGTGTTTCCCATGTGGGTCt                       <  1:3152914/36‑1 (MQ=21)
          gtgtAATGTCGCTGGCGGTGTTTTGCATTTGGGGCt                       >  1:3025907/1‑36 (MQ=16)
          gtgtAATGTCGCTGGCGGTGTTTACCATGTGGGTCt                       <  1:808562/36‑1 (MQ=255)
          gtgtAATGTCGCTGGCGGTGGTTTTCATGTGGGttt                       >  1:3985174/1‑34 (MQ=21)
             tAATGTCGCTGGCGGTGTTTACCATGTGGGTCttta                    >  1:204361/1‑35 (MQ=255)
              aaTGTCGCTGGCGGTGTTTTCCATGTGGGTCttttt                   >  1:5156302/1‑35 (MQ=255)
               aTGTCGCTGGCGGTGTTTGCCATTTGGGTCtttttt                  >  1:1540468/1‑34 (MQ=25)
                 gTCGCTGGCGGTGTTTGCCATGTGGGTCTTTTCTCg                <  1:887940/36‑1 (MQ=255)
                   cGCTGGCGGTGTTTACCATGTTGGTCTTTTCTCGCt              <  1:994654/36‑1 (MQ=255)
                   cGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCt              <  1:978039/36‑1 (MQ=255)
                    gCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCtt             <  1:1615524/36‑1 (MQ=255)
                         cGGTGTTTACCATGTGGGTCTTTTCTCGC.TTGCAcc       <  1:5604375/36‑1 (MQ=38)
                          ggTGTTTACCATGTGGGTCTTTTCTCGCTTTGCAcc       >  1:4787571/1‑36 (MQ=255)
                           gtgtTTACCATGTGGGTCTTTTCTCGCTTTGCACCg      <  1:1854584/36‑1 (MQ=255)
                            tgtTTACCATGTGGGTCTTTTCTCGCTTTGCACCgg     >  1:5138402/1‑36 (MQ=255)
                             gtTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGa    >  1:1088672/1‑36 (MQ=255)
                              tttACCATGTGGGTCTTTTCTCGCTTTGCACCGGAt   <  1:4654302/36‑1 (MQ=255)
                              tttACCATGTGGGTCTGTTCTCGCTTTGCACCGGAt   >  1:1290752/1‑36 (MQ=255)
                               ttACCATGTGGGTCTTTTCTCGCTTTGCACCGGATg  <  1:1427065/36‑1 (MQ=255)
                               ttACCATGCGGGTCTTTTCTCGATTTGCACCGGATg  >  1:1784951/1‑36 (MQ=37)
                                 |                                 
CGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATG  >  REL606/2789521‑2789587

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