Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,978,074 C→G 13.8% L200V (CTG→GTG)  yedY → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,978,0740C→G13.8% 3.2 ‑4.1 ‑0.054598yedYhypothetical protein
Reads supporting (aligned to +/- strand):  new base (7/1):  ref base (17/28):  total (25/29)
Fisher's exact test strand distribution p-value = 1.70e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.74e-01

TGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGAC  >  REL606/1978039‑1978109
                                   |                                   
tGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGc                                     >  1:4010256/1‑36 (MQ=255)
 gccaggccagCAGGACCGCTTTATCGGCGGCGGGCt                                    <  1:2296076/36‑1 (MQ=255)
  ccaggccagCGGGACCGCTTTATCGGCGGCGGGCTg                                   <  1:6364992/36‑1 (MQ=255)
  ccaggccagCAGGACCGCTTTATCGGCGGCGGGGTg                                   >  1:5078708/1‑36 (MQ=255)
  ccaggacagCAGGACCGCTTTATCGGCGGCCGGGTg                                   >  1:5072916/1‑36 (MQ=25)
   caggccagCAGGACCGCTTTATCGGCGGCGGGCTGa                                  <  1:6442286/36‑1 (MQ=255)
    aggccagCAGGACCGCTTTATCGGCGGCGGGCTGaa                                 >  1:7748581/1‑36 (MQ=255)
     ggccagCAGGACCGCTTTATCGGCGGCGGGCTGaaa                                >  1:1612604/1‑36 (MQ=255)
     ggccagCAGGACCGCTTTAACGGCGGGGGGGTGaaa                                >  1:4479022/1‑36 (MQ=25)
       ccagCAGGACCGCTTTATCGGCGGCGGGGTGgaaaa                              >  1:2759157/1‑34 (MQ=25)
       ccagCAGGACCGCTTTATCGGCGGCGGGCTGAAata                              <  1:7227632/36‑1 (MQ=255)
        cagcagGACCGCTTTATCGGCGGCGGGGTGAAatat                             >  1:889447/1‑36 (MQ=255)
         agcagGACCGCTTTATCGGCGGCGGGCTTAAATATc                            >  1:7991245/1‑36 (MQ=255)
          gcagGACCGCTTTATCGGCGGCGGGCTGAAATATcc                           <  1:2194909/36‑1 (MQ=255)
          gcagGACCGCTTTATCGGCGGCGGGCTGAAATATcc                           <  1:5972856/36‑1 (MQ=255)
          gcagGACCGCTTTATCGGCGGCGGGCTGAAATATcc                           <  1:1184184/36‑1 (MQ=255)
           cagGACCGCTTTATCGGCGGCGGGCTGGAATATCCt                          >  1:3220127/1‑36 (MQ=255)
           cagGACCGCTTTATCGGCGGCGGGCTGAAATATCCt                          <  1:2603789/36‑1 (MQ=255)
           cagGACCGCTTTATCGGCGGCGGGATGGAATATTCt                          >  1:4817617/1‑36 (MQ=25)
             ggACCGCTTTATCGGCGGCGGGCTGAAATATCCTTa                        <  1:225293/36‑1 (MQ=255)
              gACCGCTTTATCGGCGGCGGGGTGAAATATCCTTAt                       >  1:3150624/1‑36 (MQ=255)
               aCCGCTTTATCGGCGGCGGGCTGAAATATCCTTATg                      <  1:7964067/36‑1 (MQ=255)
               aCCGCTTTATCGGCGGCGGGCTGAAATATCCTTATg                      >  1:2286053/1‑36 (MQ=255)
                 cGCTTTATCGGCGGCGGGCTGAAATATCCTTATGTc                    <  1:5931806/36‑1 (MQ=255)
                  gCTTTATCGGCGGCGGGCTGAAATATCCTTATGTCg                   <  1:4445535/36‑1 (MQ=255)
                  gCTTTATCGGCGGCGGGCTGAAATATCCTTATGTCg                   <  1:8470804/36‑1 (MQ=255)
                   cTTTATCGGCGGCGGGCTGAAATATCCTTATGTCGa                  <  1:1962992/36‑1 (MQ=255)
                    tttATCGGCGGCGGGCTGAAATATCCTTATGTCGaa                 <  1:1808076/36‑1 (MQ=255)
                     ttATCGGCGGCGGGGTGAAATATTCTTTTGTTGAAg                >  1:5732420/1‑36 (MQ=16)
                     ttATCGGCGGCGGGCTGGAAAATTCTTATGTCGCAg                >  1:5332081/1‑36 (MQ=16)
                     ttATCGGCGGCGGGCTGAAATATCCTTATGTTGAAg                >  1:6946001/1‑36 (MQ=255)
                      tATCGGCGGCGGGCTGAAATATCCTTATGTCGAAgg               <  1:1719513/36‑1 (MQ=255)
                       aTCGGCGGCGGGCTGAAATATCCTTATGTCGAAGGa              <  1:8305741/36‑1 (MQ=255)
                        tCGGCGGCGGGCTGAAATATCCTTATGTCGAAGGAt             <  1:1155482/36‑1 (MQ=255)
                        tCGGCGGCGGGCTGAAATATCCTTATGTCGAAGGAt             <  1:2304163/36‑1 (MQ=255)
                        tCGGCGGCGGGCTGAAATATCCTTATGTCGAAGGAt             >  1:4177634/1‑36 (MQ=255)
                         cggcggcggGGTGAAATATCCTTATGTCGAAGGAtt            <  1:7358653/36‑1 (MQ=255)
                         cggcggcggGCTGAAATATCCTTACGGCGAAGGAtt            >  1:5922295/1‑36 (MQ=37)
                          ggcggcggGCTGAAATATCCTTATGTCGAAGGATTg           >  1:1698850/1‑36 (MQ=255)
                          ggcggcggGCTGAAATATCCTTATGTCGAAGGATTg           <  1:5750484/36‑1 (MQ=255)
                          ggcggcggGCTGAAATATCCTTATGTCGAAGGATTg           <  1:7810744/36‑1 (MQ=255)
                          ggcggcggGCTGAAATATCCTTATGTCGAAGGATTg           >  1:6680825/1‑36 (MQ=255)
                           gcggcggGCTGAAATATCCTTATGTCGAAGGATTGc          <  1:1357105/36‑1 (MQ=255)
                           gcggcggGCTGAAATATCCTTATGTCGAAGGATTGc          <  1:4313535/36‑1 (MQ=255)
                           gcggcggGCTGAAATATCCTTATGTCGAAGGATTGc          <  1:4215999/36‑1 (MQ=255)
                            cggcggGCTGAAATATCCTTATGTCGAAGGATTGCg         <  1:3945359/36‑1 (MQ=255)
                             ggcggGCTGAAATATCCTTATGTCGAAGGATTGCGt        >  1:7442448/1‑36 (MQ=255)
                             ggcggGCTGAAATATCCTTATGGCGAAGGATTGCGt        >  1:48844/1‑36 (MQ=255)
                               cggGCTGAAATATCCTTATGTCGAAGGATTGCGtct      >  1:6794682/1‑36 (MQ=255)
                               cggGCTGAAATATCCTTATGTCGAAGGATTGCGtct      <  1:5387610/36‑1 (MQ=255)
                                gggCTGAAATATCCTTATGTCGAAGGATTGCGtctc     >  1:620168/1‑36 (MQ=255)
                                 ggCTGAAATATCCTTATGTCGAAGGATTGCGTCTcg    <  1:3661282/36‑1 (MQ=255)
                                  gCTGAAATATCCTTATGTCGAAGGATTGCGTCTcga   >  1:4456/1‑36 (MQ=255)
                                  gCTGAAATATCCTTATGTCGAAGGATTGCGTCTcga   <  1:5939272/36‑1 (MQ=255)
                                  gCTGAAATATCCTTATGACGAAGGATTGCGTCTcga   >  1:117462/1‑36 (MQ=255)
                                   cTGAAATATCCTTATGTCGAAGGATTGCGTCTcgac  <  1:698242/36‑1 (MQ=255)
                                   |                                   
TGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGAC  >  REL606/1978039‑1978109

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