Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,573,574 G→C 13.8% G83G (GGC→GGG yhiQ ← predicted SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,573,5740G→C13.8% 2.5 ‑2.5 0.042249yhiQpredicted SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (14/21):  total (15/27)
Fisher's exact test strand distribution p-value = 7.85e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCA  >  REL606/3573539‑3573607
                                   |                                 
gaaTCCACCACATCCGGCAAATAATCGCCTTTAATg                                   <  1:7006984/34‑1 (MQ=255)
gCATCCACCACATCCGGCAAATAATCGCCTTTAATg                                   >  1:1149301/1‑36 (MQ=255)
 cATCCACCACATCCGGCAAATAATCGCCTTTAATGc                                  <  1:1621788/36‑1 (MQ=255)
  aTCCACCACATCCGGCAAATTATTGCCTTTAAGGcc                                 >  1:3154195/1‑36 (MQ=25)
  aTCCACCACATCCGGCAAATAATCGCCTTTAATTcc                                 >  1:3162691/1‑36 (MQ=255)
  aTCCACCACATCCGGCAAATAATCGCCTTTAATGcc                                 >  1:1442205/1‑36 (MQ=255)
  aTCCACCACATCCGGCAAATAATCGCCTTTAATGcc                                 <  1:6553640/36‑1 (MQ=255)
  aTCCACCACATCCGGCAAATAATCGCCTTTAATGcc                                 <  1:5579254/36‑1 (MQ=255)
    ccaccaCATCCGGCAAATAATCGCCTTTAATGCCCa                               >  1:1183582/1‑36 (MQ=255)
     caccacATCCGGCAAATAATCGCCTTTAATGCCCAc                              >  1:5690508/1‑36 (MQ=255)
     caccacATCCGGCAAATAATCGCCTTTAATGCCCAc                              <  1:4387211/36‑1 (MQ=255)
      accacATCCGGCAAATAATCGCCTTTAATGCCCAcc                             <  1:1853385/36‑1 (MQ=255)
       ccacATCCGGCAAATAATCGCCTTTAATGCCCACCg                            <  1:7403465/36‑1 (MQ=255)
        cacaTCCGGCAAATAATCGCCTTTAATGCCCACCGc                           >  1:1734443/1‑36 (MQ=255)
         acaTCCGGCAAATAATCGCCTTTAATGCCCACCGCt                          <  1:3879592/36‑1 (MQ=255)
         acaTCCGGCAAATAATCGCCTTTAATGCCCACCGCt                          <  1:5392223/36‑1 (MQ=255)
           aTCCGGCAAATAATCGCCTTTAATGCCCACCGCttt                        >  1:2615804/1‑36 (MQ=255)
           aTCCGGCAAATAATCGCCTTTAATGCCCACCGCttt                        <  1:6859318/36‑1 (MQ=255)
           aTCCGGCAAATAATCGCCTTTAATGCCCACCGCttt                        <  1:363830/36‑1 (MQ=255)
              cgtcAAATAATCGCCTTTAATGCCCACCGCTTTCGc                     <  1:5731521/33‑1 (MQ=255)
              cGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGc                     >  1:1266312/1‑36 (MQ=255)
              cGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGc                     >  1:5933802/1‑36 (MQ=255)
               ggCAAATAATCGCCTTTAATGCCCACCGCTTTCGcc                    <  1:3264907/36‑1 (MQ=255)
               ggCAAATAATCGCCTTTAATGCCCACCGCTTTCGcc                    <  1:6067734/36‑1 (MQ=255)
                gCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCa                   <  1:2514385/36‑1 (MQ=255)
                gCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCa                   <  1:1478360/36‑1 (MQ=255)
                  aaaTAATCGCCTTTAATGCCCACCGCTTTCGCCAcc                 <  1:5654888/36‑1 (MQ=255)
                  aaaTAATCGCCTTTAATGCCCACCGCTTTCGCCAcc                 <  1:1625030/36‑1 (MQ=255)
                  aaaTAATCACCTTTAATGCCCTCGGCTTTCGCCAcc                 <  1:5652434/36‑1 (MQ=21)
                   aataatCGCCTTTAATGCCCACCGCTTTCGCCACCg                <  1:3457621/36‑1 (MQ=255)
                      aatCGCCTTTAATGCCCACCGCTTTCGCCACCGCCt             >  1:6626326/1‑36 (MQ=255)
                      aatCGCCTTTAATGCCCACCGCTTTCGCCACCGCCt             >  1:6532402/1‑36 (MQ=255)
                         cGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCAc          <  1:7513991/36‑1 (MQ=255)
                         cGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCAc          <  1:4809955/36‑1 (MQ=255)
                          gCCTTAAATCCCCACCGCTTTCCCCACCGCCTCAcc         <  1:5404239/36‑1 (MQ=25)
                          cccTTTAACCCCCCCCGTTTCCGCCACCGCCTCAcc         <  1:2782578/35‑1 (MQ=16)
                           ccTTAAATCCCCCCCGCTTTCGCCACCGCCTCACcg        <  1:4851784/36‑1 (MQ=25)
                            cTTTAACCCCCCCCGCTTTCCCCACCGCCTCACcgc       <  1:3824477/36‑1 (MQ=19)
                            cTTAAATCCCCCCCGCTTACACCACCGCCTCACcgc       <  1:2022114/36‑1 (MQ=16)
                              ttAATGCCCACCGCTTTCGCCACCGCCTCACCGCGa     >  1:5771435/1‑36 (MQ=255)
                                aaTCCCCACCGCTTCCGCCACCGCCTCACCGCGacc   <  1:1546145/36‑1 (MQ=37)
                                 aTGCCCACCGCTTTCGCCACCGCCTCACCGCGaccc  >  1:3249715/1‑35 (MQ=255)
                                 aTGCCCACCGCTTTCGCCACCGCCTCACCGCGacca  >  1:6375982/1‑36 (MQ=255)
                                 aTGCCCACCGCTTTCGCCACCGCCTCACCGCGacca  >  1:675952/1‑36 (MQ=255)
                                   |                                 
GCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCA  >  REL606/3573539‑3573607

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