Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,084,514 T→C 11.8% D137G (GAT→GGT)  yegK ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,084,5140T→C11.8% 4.1 ‑4.1 0.062410yegKhypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (25/29):  total (25/37)
Fisher's exact test strand distribution p-value = 1.70e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGGCGTGAGCGGCAGTTCCAGACCACGACCAAAATCCACCACTACACCGCCGTCGCCAATTTGCATAATC  >  REL606/2084480‑2084549
                                  |                                   
cGGCGTAAGCGGTAGTTCCAGACCACGACCAAAATc                                    <  1:4945915/36‑1 (MQ=38)
  gCGTTAGCGGCAGTTCCAGACCACGACCAAAATcca                                  <  1:7042018/36‑1 (MQ=255)
    gTGAGCGGCAGTTCCAGACCACGACCAAAATccacc                                <  1:1283788/36‑1 (MQ=255)
      gAGCGGCAGTTCCCGACCACGACCAAAATCcaccac                              >  1:1574514/1‑36 (MQ=255)
      gAGCGGCAGTTCCAGACCACGACCAAAATCcgccac                              <  1:2448774/36‑1 (MQ=255)
      gAGCGGCAGTTCCAGACCACGACCAAAATCcaccac                              <  1:856234/36‑1 (MQ=255)
      gAGCGGCAGTTCCAGACCACGACCAAAATCcaccac                              <  1:6479730/36‑1 (MQ=255)
       aGCGGCAGTTCCAGACCACGACCAAAATCCACCACt                             <  1:5333590/36‑1 (MQ=255)
        gCGGCAGTTCCAGACCACGACCAAAATCCACCACTa                            <  1:4649407/36‑1 (MQ=255)
         cGGCAGTTCCAGACCACGACCAAAATCCACCACTac                           >  1:4592648/1‑36 (MQ=255)
         cGGCAGTTCCAGACCACGACCAAAATCCACCACTac                           <  1:3545654/36‑1 (MQ=255)
         cGGCAGTTCCAGACCACGACCAAAATCCACCACTac                           <  1:2077296/36‑1 (MQ=255)
         cGGCAGTTCCAGACCACGACCAAAATCCACCACTac                           <  1:826151/36‑1 (MQ=255)
          ggCAGTTCCAGACCACGACCAAAATCCACCACTaca                          >  1:2079925/1‑36 (MQ=255)
          ggCAGTTCCAGACCACGACCAAAATCCACCACTaca                          <  1:2529145/36‑1 (MQ=255)
           gCAGTTCCAGACCACGACCAAAATCCACCACTacac                         >  1:2605898/1‑36 (MQ=255)
           gCAGTTCCAGACCACGACCAAAATCCACCACTacac                         <  1:5028628/36‑1 (MQ=255)
           gCAGTTCCAGACCACGACCAAAATCCACCACTacac                         >  1:1348284/1‑36 (MQ=255)
           gCAGTTCCAGACCACGACCAAAATCCACCACTacac                         <  1:6226301/36‑1 (MQ=255)
           gCAGTTCCAAACCACGACCAAAATCCACCACTacac                         <  1:5201575/36‑1 (MQ=255)
            cAGTTCCAGACCACGACCAAAATCCACCACTACAcc                        <  1:1631671/36‑1 (MQ=255)
               ttCCAGACCTCGACCAAAATCCACCACTACAccgcc                     <  1:4560076/36‑1 (MQ=255)
               ttCCAGACCACGACCAAAATCCACCACTACAccgcc                     <  1:7093422/36‑1 (MQ=255)
                ccaaGACAACGACAAAAACCCACCACTACAccgccg                    <  1:2261263/33‑1 (MQ=19)
                 ccAGACCACGACCAAAATCCACCACTACAACGCCGt                   >  1:2548545/1‑36 (MQ=255)
                  cAGACCACGACCAAAATCCACCACTACACCGCCGTc                  >  1:6034061/1‑36 (MQ=255)
                  cAGACCACGACCAAAATCCACCACTACACCGCCGTc                  >  1:1880810/1‑36 (MQ=255)
                  cAGACCAAGACCAAAAT.CACCACTACCCCCCCGTcc                 >  1:6588420/1‑35 (MQ=16)
                   aGACCACGACCAAAATCCACCACTACACCGCCGTCg                 <  1:1383396/36‑1 (MQ=255)
                   aGACCACGACCAAAATCCACCACTACACCGCCGTCg                 >  1:3959628/1‑36 (MQ=255)
                   aGACCACGACCAAAATCCACCACTACACCGCCGTCg                 >  1:421443/1‑36 (MQ=255)
                    gACCACGACCAAAATCCACCACTACACCGCCGTCGc                >  1:2093334/1‑36 (MQ=255)
                     aCCACGACCAAAATCCACCACTACACCGCCGTCGcc               >  1:6685437/1‑36 (MQ=255)
                     aCCACAACAAAAACCCACCACTACACCGCCGTCGcc               <  1:4886661/36‑1 (MQ=25)
                      ccACGACCAAAATCCACCACTACACCGCCGTCGCCa              >  1:7117940/1‑36 (MQ=255)
                      ccACGACAAAAATCCACCACTACACCGCCGTCGCCa              <  1:1295618/36‑1 (MQ=255)
                      ccACAACAAAAATCCACCACTACCCCGCCGTCGCCa              <  1:3925723/36‑1 (MQ=25)
                       cACGACCAAAATCCACCACTACACCGCCGTCGCCac             >  1:3468951/1‑35 (MQ=255)
                       cACGACCAAAATCCACCACTACACCGCCGTCGCCaa             >  1:3465406/1‑36 (MQ=255)
                       cACCACAAAAATCCACCACTACTCCGCCGTCGCCaa             <  1:4025837/36‑1 (MQ=25)
                       aaCGACAAAAATCCACCACTACACCGCCGTCGCCaa             <  1:5551129/35‑1 (MQ=255)
                        aCGACCAAAATCCACCACTACACCGCCGTCGCCAAt            >  1:839488/1‑36 (MQ=255)
                        aCGACCAAAATCCACCACTACACCGCCGTCGCCAAt            >  1:6063731/1‑36 (MQ=255)
                        aCGACAAAAATCCACCACTACACCGCCGTCGCCAAt            <  1:1914043/36‑1 (MQ=255)
                        aCGACAAAAACCCACCACTGCACCGCCGTCGCCAAt            <  1:6609754/36‑1 (MQ=25)
                        aCGACAAAAACCCACCACTACCCCGCCGTCGCCAAt            <  1:3801346/36‑1 (MQ=25)
                         agACAAAAACCCCCCACTACACCGCCGTCGCCAAtt           <  1:3734783/35‑1 (MQ=25)
                          gACCAAAATCCACCACTACACCGCCGTCGCCAAttt          >  1:997034/1‑36 (MQ=255)
                          gACCAAAATCCACCACTACACCGCCGTCGCCAAttt          >  1:5246472/1‑36 (MQ=255)
                          gACCAAAATCCACCACTACACCGCCGTCGCCAAttt          >  1:4632980/1‑36 (MQ=255)
                          gACCAAAATCCACCACTACACCGCCGTCGCCAAttt          >  1:1860563/1‑36 (MQ=255)
                          gACAAAAACCCACCACTACACCGCCGTCGCCAAttt          <  1:4529913/36‑1 (MQ=37)
                           acaaaaaTCCACCACTACACCGCCGTCGCCAATTTg         <  1:3604864/33‑1 (MQ=255)
                           acaaaaaTCCACCACTACACCGCCGTCGCCAATTTg         <  1:1529634/33‑1 (MQ=255)
                           aCCAAAATCCACCACTACACCGCCGTCGCCAATTTg         >  1:1993584/1‑36 (MQ=255)
                            ccAAAATCCACCACTACACCGCCGTCGCCAATTTGc        >  1:4581649/1‑36 (MQ=255)
                            ccAAAATCCACCACTACACCGCCGTCGCCAATTTGc        <  1:2228808/36‑1 (MQ=255)
                            ccAAAACCCACCACTACACCGCCGTCGCCAATTTGc        <  1:4940892/36‑1 (MQ=255)
                             aaaaaCCCACCACTACACCGCCGTCGCCAATTTGCa       <  1:1725649/35‑1 (MQ=255)
                              aaaaTCCACCACTACACCGCCGTCGCCAATTTGCAt      <  1:434876/36‑1 (MQ=255)
                                 aTCCACCACTACACCGCCGTCGCCAATTTGCATAAt   >  1:3599971/1‑36 (MQ=255)
                                 aTCAACCACTACACCGCCGTCGCCAATTTGCATAAt   <  1:5741061/36‑1 (MQ=255)
                                  tCCACCACTACACCGCCGTCGCCAATTTGCATAATc  <  1:1133018/36‑1 (MQ=255)
                                  tCCACCACTACACCGCCGTCGCCAATTTGCATAATc  >  1:52472/1‑36 (MQ=255)
                                  tCCACCACAGCACCGCCGTCGCCAATTTGCATAATc  <  1:7318830/36‑1 (MQ=37)
                                  |                                   
CGGCGTGAGCGGCAGTTCCAGACCACGACCAAAATCCACCACTACACCGCCGTCGCCAATTTGCATAATC  >  REL606/2084480‑2084549

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