Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,558,775 Δ1 bp 100% coding (663/1425 nt) hsdS ← specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,558,7750T.100.0% 594.8 / NA 120coding (663/1425 nt)hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (49/71);  total (49/71)
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 152 total aligned reads displayed.

GTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTA  >  REL606/4558725‑4558821
                                                  |                                              
gTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACt                                                <  1:8477298/51‑1 (MQ=255)
gTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACt                                                <  1:3477024/51‑1 (MQ=255)
gTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCCAATTCAACt                                                <  1:5047122/51‑1 (MQ=255)
 tttCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACtt                                               <  1:6780805/51‑1 (MQ=255)
 tttCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACtt                                               >  1:5146644/1‑51 (MQ=255)
 tttCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACtt                                               >  1:7124747/1‑51 (MQ=255)
 tttCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACtt                                               >  1:2516427/1‑51 (MQ=255)
 tttCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACtt                                               <  1:4337371/51‑1 (MQ=255)
  ttCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttt                                              <  1:4249845/51‑1 (MQ=255)
  ttCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttt                                              <  1:8188681/51‑1 (MQ=255)
  ttCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttt                                              <  1:2397079/51‑1 (MQ=255)
    cACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttt                                            >  1:9883731/1‑51 (MQ=255)
    cACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttt                                            >  1:5590176/1‑51 (MQ=255)
     aCTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACtttttt                                           <  1:8641226/51‑1 (MQ=255)
     aCTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACtttttt                                           >  1:8435933/1‑51 (MQ=255)
     aCTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACtttttt                                           >  1:9036157/1‑51 (MQ=255)
        ccAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttatt                                        <  1:6873637/51‑4 (MQ=255)
        ccAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttatt                                        >  1:193919/1‑48 (MQ=255)
        ccAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttatt                                        <  1:4042109/51‑4 (MQ=255)
         cAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttatta                                       <  1:865899/51‑5 (MQ=255)
         cAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttatta                                       <  1:4461888/51‑5 (MQ=255)
         cAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttatta                                       <  1:6485989/51‑5 (MQ=255)
         cAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttatta                                       <  1:6388855/51‑5 (MQ=255)
          aGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTtattat                                     <  1:3653897/51‑1 (MQ=255)
            tcttctGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATaa                                   <  1:1028565/51‑1 (MQ=255)
             cttctGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATtata                                    <  1:6072412/49‑1 (MQ=39)
             cttctGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATaaa                                  >  1:8071729/1‑51 (MQ=255)
              ttctGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATaaaa                                 >  1:7400811/1‑51 (MQ=255)
              ttctGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATaaaa                                 >  1:7252793/1‑51 (MQ=255)
               tctGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAAt                                <  1:1207396/51‑1 (MQ=255)
               tctGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAAt                                <  1:2859164/51‑1 (MQ=255)
               tctGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAAt                                <  1:5507478/51‑1 (MQ=255)
                ctGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATc                               <  1:3077559/51‑1 (MQ=255)
                 tGGAATAGGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCt                              <  1:23629/51‑1 (MQ=255)
                 tGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCt                              <  1:3642940/51‑1 (MQ=255)
                 tGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCt                              <  1:5609483/51‑1 (MQ=255)
                 tGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCt                              <  1:3731827/51‑1 (MQ=255)
                  ggAATAAGGATTTTGTATTTATCTAATTCAACTTTTTTATTATAAAATCtt                             <  1:1803319/51‑1 (MQ=255)
                  ggAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCtt                             >  1:3796447/1‑51 (MQ=255)
                  ggAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCtt                             >  1:5734516/1‑51 (MQ=255)
                  ggAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCtt                             <  1:1690206/51‑1 (MQ=255)
                   gAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCttt                            <  1:1803317/51‑1 (MQ=255)
                   gAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCttt                            <  1:10555536/51‑1 (MQ=255)
                   gAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCttt                            >  1:4074558/1‑51 (MQ=255)
                   gAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCttt                            <  1:192469/51‑1 (MQ=255)
                   gAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCttt                            >  1:983550/1‑51 (MQ=255)
                   gAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCttt                            <  1:8827633/51‑1 (MQ=255)
                   gAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCttt                            <  1:3473825/51‑1 (MQ=255)
                    aaTAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTa                           <  1:6604376/51‑1 (MQ=255)
                     aTAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAt                          <  1:10260295/51‑1 (MQ=255)
                     aTAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAt                          <  1:9705575/51‑1 (MQ=255)
                     aTAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAt                          <  1:4768210/51‑1 (MQ=255)
                      tAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtc                         <  1:3571870/51‑1 (MQ=255)
                      tAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtc                         >  1:4886661/1‑51 (MQ=255)
                       aaGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtct                        >  1:2787288/1‑51 (MQ=255)
                       aaGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtct                        <  1:4528967/51‑1 (MQ=255)
                       aaGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtct                        <  1:5213544/51‑1 (MQ=255)
                        aGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtctt                       <  1:8372378/51‑1 (MQ=255)
                        aGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtctt                       <  1:6606648/51‑1 (MQ=255)
                        aGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtctt                       <  1:3710040/51‑1 (MQ=255)
                         ggATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtcttc                      <  1:184261/51‑1 (MQ=255)
                         ggATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtcttc                      <  1:4598415/51‑1 (MQ=255)
                          gATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtcttct                     >  1:797863/1‑51 (MQ=255)
                           aTTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATcttctt                    <  1:1096862/51‑1 (MQ=255)
                           aTTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATcttctt                    >  1:7142989/1‑51 (MQ=255)
                            ttttGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCttt                   <  1:261845/51‑1 (MQ=255)
                            ttttGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCttt                   <  1:10203320/51‑1 (MQ=255)
                            ttttGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCttt                   <  1:758446/51‑1 (MQ=255)
                            ttttGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCttt                   >  1:2488595/1‑51 (MQ=255)
                            ttttGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCttt                   <  1:4920617/51‑1 (MQ=255)
                              ttGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTgg                 >  1:1860262/1‑51 (MQ=255)
                               tGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGt                >  1:1292742/1‑51 (MQ=255)
                               tGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGt                <  1:8961458/51‑1 (MQ=255)
                               tGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGt                >  1:6021919/1‑51 (MQ=255)
                                gTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTc               <  1:3172240/51‑1 (MQ=255)
                                gTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTc               >  1:4036114/1‑51 (MQ=255)
                                 tAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCa              >  1:6546519/1‑51 (MQ=255)
                                 tAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCa              >  1:2854707/1‑51 (MQ=255)
                                 tAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCa              >  1:7711940/1‑51 (MQ=255)
                                  aGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCaa             >  1:3603826/1‑51 (MQ=255)
                                  aGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCaa             <  1:8733396/51‑1 (MQ=255)
                                    ttATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAAcc           >  1:4492068/1‑51 (MQ=255)
                                    ttATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAAcc           >  1:6675520/1‑51 (MQ=255)
                                    ttATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAAcc           <  1:8228430/51‑1 (MQ=255)
                                     tATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCt          >  1:8156118/1‑51 (MQ=255)
                                      aTCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTa         <  1:4632213/51‑1 (MQ=255)
                                      aTCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTa         <  1:9371445/51‑1 (MQ=255)
                                      aTCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTa         >  1:3937107/1‑51 (MQ=255)
                                      aTCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTa         <  1:3870254/51‑1 (MQ=255)
                                       tCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTAc        <  1:948757/51‑1 (MQ=255)
                                        cTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTAcc       >  1:5873275/1‑51 (MQ=255)
                                        cTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTAcc       >  1:6369259/1‑51 (MQ=255)
                                         tAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCa      <  1:3820842/51‑1 (MQ=255)
                                         tAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCa      >  1:10191751/1‑51 (MQ=255)
                                         tAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCa      >  1:10608516/1‑51 (MQ=255)
                                         tAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCa      >  1:4643599/1‑51 (MQ=255)
                                          aaTTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAg     <  1:10918905/51‑1 (MQ=255)
                                            ttCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGtt   >  1:5674559/1‑51 (MQ=255)
                                            ttCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGtt   <  1:2445628/51‑1 (MQ=255)
                                             tCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTa  >  1:1515733/1‑51 (MQ=255)
                                                  |                                              
GTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTA  >  REL606/4558725‑4558821

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