Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,308,106 G→A 100% A6A (GCG→GCA yhcB → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,308,1060GA100.0% 396.5 / NA 116A6A (GCG→GCAyhcBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (57/59);  total (57/59)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 116 total aligned reads displayed.

CTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCG  >  REL606/3308058‑3308153
                                                |                                               
cTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACt                                               <  1:479419/51‑1 (MQ=255)
 tACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTa                                              >  1:4144844/1‑51 (MQ=255)
 tACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTa                                              <  1:6038378/51‑1 (MQ=255)
 tACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTa                                              >  1:10498147/1‑51 (MQ=255)
  aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             >  1:3397451/1‑51 (MQ=255)
  aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             <  1:4558825/51‑1 (MQ=255)
  aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             >  1:6911426/1‑51 (MQ=255)
  aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             <  1:1949863/51‑1 (MQ=255)
  aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             <  1:4269558/51‑1 (MQ=255)
   cGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAAt                                            >  1:7377723/1‑51 (MQ=255)
    ggAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAAtt                                           <  1:5059443/51‑1 (MQ=255)
     gAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTg                                          >  1:5374061/1‑51 (MQ=255)
     gAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTg                                          <  1:10101889/51‑1 (MQ=255)
     gAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTg                                          <  1:7192504/51‑1 (MQ=255)
       aaCAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTggg                                        >  1:6062233/1‑51 (MQ=255)
       aaCAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTggg                                        <  1:10334365/51‑1 (MQ=255)
       aaCAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTggg                                        <  1:4423004/51‑1 (MQ=255)
        aCAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGt                                       >  1:5051921/1‑51 (MQ=255)
         cAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGtt                                      >  1:4948040/1‑51 (MQ=255)
         cAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGtt                                      >  1:6768410/1‑51 (MQ=255)
         cAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGtt                                      >  1:1958060/1‑51 (MQ=255)
          aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     <  1:5448649/51‑1 (MQ=255)
          aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     >  1:2237172/1‑51 (MQ=255)
          aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     >  1:9207752/1‑51 (MQ=255)
          aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     >  1:4825146/1‑51 (MQ=255)
          aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     >  1:889268/1‑51 (MQ=255)
          aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     >  1:5439388/1‑51 (MQ=255)
           aaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAg                                    <  1:3268049/51‑1 (MQ=255)
           aaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAg                                    >  1:5662585/1‑51 (MQ=255)
           aaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAg                                    >  1:1412667/1‑51 (MQ=255)
            aaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgt                                   >  1:7236411/1‑51 (MQ=255)
            aaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgt                                   <  1:1394189/51‑1 (MQ=255)
            aaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgt                                   <  1:3916041/51‑1 (MQ=255)
             aGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                  >  1:30401/1‑51 (MQ=255)
             aGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                  >  1:4475438/1‑51 (MQ=255)
             aGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                  <  1:5826930/51‑1 (MQ=255)
             aGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                  >  1:327147/1‑51 (MQ=255)
              gACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcg                                 >  1:7981265/1‑51 (MQ=255)
               aCAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcgt                                <  1:6802174/51‑1 (MQ=255)
                cAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcgtc                               <  1:3870763/51‑1 (MQ=255)
                cAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcgtc                               <  1:6574489/51‑1 (MQ=255)
                 aaCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGtcgtcg                              >  1:9888199/1‑51 (MQ=255)
                 aaCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGtcgtcg                              <  1:5184232/51‑1 (MQ=255)
                 aaCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGtcgtcg                              <  1:9105478/51‑1 (MQ=255)
                 aaCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGtcgtcg                              <  1:3258823/51‑1 (MQ=255)
                  aCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCgg                             <  1:10691680/51‑1 (MQ=255)
                  aCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCgg                             >  1:3102172/1‑51 (MQ=255)
                  aCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCgg                             >  1:6138568/1‑51 (MQ=255)
                   cGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGc                            <  1:8951543/51‑1 (MQ=255)
                    gggAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGca                           >  1:3950778/1‑51 (MQ=255)
                       agaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatca                        <  1:6010/51‑1 (MQ=255)
                       agaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatca                        <  1:10108917/51‑1 (MQ=255)
                       agaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatca                        >  1:6051567/1‑51 (MQ=255)
                       agaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatca                        <  1:10828078/51‑1 (MQ=255)
                        gaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatcat                       <  1:10090365/51‑1 (MQ=255)
                        gaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatcat                       <  1:590000/51‑1 (MQ=255)
                         aTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCatt                      <  1:4675863/51‑1 (MQ=255)
                         aTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCatt                      >  1:3896898/1‑51 (MQ=255)
                         aTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCatt                      >  1:9183342/1‑51 (MQ=255)
                           gTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatc                         <  1:9379349/46‑1 (MQ=255)
                           gTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatc                         <  1:5880414/46‑1 (MQ=255)
                           gTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCattat                    >  1:4931266/1‑51 (MQ=255)
                            ttCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCattatt                   >  1:8714388/1‑51 (MQ=255)
                             tCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTg                  <  1:10251305/51‑1 (MQ=255)
                             tCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTg                  <  1:10076468/51‑1 (MQ=255)
                              cATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTgg                 <  1:5517785/51‑1 (MQ=255)
                               aTGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGt                >  1:9941823/1‑51 (MQ=255)
                                tGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTg               >  1:1515610/1‑51 (MQ=255)
                                tGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTg               <  1:1693295/51‑1 (MQ=255)
                                 gACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGc              <  1:8495917/51‑1 (MQ=255)
                                  aCCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCt             >  1:2059598/1‑51 (MQ=255)
                                  aCCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCt             <  1:5697433/51‑1 (MQ=255)
                                   ccTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtg            <  1:485044/51‑1 (MQ=255)
                                   ccTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCt             <  1:8467786/50‑1 (MQ=255)
                                    cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgt           >  1:3718431/1‑51 (MQ=255)
                                    cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgt           >  1:1320321/1‑51 (MQ=255)
                                    cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgt           >  1:9885820/1‑51 (MQ=255)
                                    cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgt           <  1:7743009/51‑1 (MQ=255)
                                    cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgt           >  1:7441404/1‑51 (MQ=255)
                                     tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          >  1:5976856/1‑51 (MQ=255)
                                     tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          >  1:9047226/1‑51 (MQ=255)
                                     tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          >  1:722898/1‑51 (MQ=255)
                                     tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          >  1:8584457/1‑51 (MQ=255)
                                     tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          <  1:1101398/51‑1 (MQ=255)
                                     tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          >  1:1529778/1‑51 (MQ=255)
                                      gggAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTgg         <  1:6939116/51‑1 (MQ=255)
                                      gggAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTgg         >  1:10607672/1‑51 (MQ=255)
                                       ggAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGc        <  1:5018023/51‑1 (MQ=255)
                                         aaTATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCa      >  1:2871948/1‑51 (MQ=255)
                                         aaTATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCa      <  1:2461494/51‑1 (MQ=255)
                                          atatGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCAt     >  1:9423750/1‑51 (MQ=255)
                                          atatGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCAt     >  1:7188263/1‑51 (MQ=255)
                                          atatGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCAt     >  1:3952858/1‑51 (MQ=255)
                                          atatGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCAt     >  1:10388008/1‑51 (MQ=255)
                                           tatGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATg    <  1:10103634/51‑1 (MQ=255)
                                            atGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGc   <  1:7064885/51‑1 (MQ=255)
                                            atGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGc   <  1:9852818/51‑1 (MQ=255)
                                             tGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCg  <  1:2998446/51‑1 (MQ=255)
                                             tGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCg  >  1:494198/1‑51 (MQ=255)
                                             tGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCg  <  1:5462152/51‑1 (MQ=255)
                                                |                                               
CTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCG  >  REL606/3308058‑3308153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: