Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,561,146 G→A intergenic (‑123/+78) hsdM ← / ← hsdR DNA methylase M/endonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,561,1460GA100.0% 883.5 / NA 215intergenic (‑123/+78)hsdM/hsdRDNA methylase M/endonuclease R
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (113/102);  total (113/102)

Only 100 of 224 total aligned reads displayed.

GAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGGAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTAT  >  REL606/4561073‑4561218
                                                                         |                                                                        
gaagaaATTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAg                                                                        >  1:6217358/1‑76 (MQ=255)
 aagaaATTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGc                                                                       >  2:4309618/1‑76 (MQ=255)
  agaaATTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCg                                                                      <  1:7289393/76‑1 (MQ=255)
  agaaATTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCg                                                                      >  1:1753340/1‑76 (MQ=255)
   gaaCTTGCCGAGAGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGt                                                                     <  1:6547345/76‑1 (MQ=21)
   gaaAGTGCCGAATGTAATATACACAAAATCAGTCAGATTTCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGt                                                                     <  2:4647251/76‑1 (MQ=18)
    aaaTTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTa                                                                    <  2:3017882/76‑1 (MQ=255)
    aaaTTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTa                                                                    <  2:2820609/76‑1 (MQ=255)
     aaTTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTAc                                                                   <  2:2126890/76‑1 (MQ=255)
      aTTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACa                                                                  <  2:3384228/76‑1 (MQ=255)
        tGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGa                                                                >  1:5562053/1‑76 (MQ=255)
          ccGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAg                                                              >  1:5213813/1‑76 (MQ=255)
            gAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTa                                                            >  1:484263/1‑76 (MQ=255)
            gAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTa                                                            >  2:332934/1‑76 (MQ=255)
             aGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTAc                                                           >  2:4570804/1‑76 (MQ=255)
             aGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTAc                                                           >  2:7622935/1‑76 (MQ=255)
              gggTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACt                                                          <  2:722533/76‑1 (MQ=255)
                gTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTg                                                        <  2:464808/76‑1 (MQ=255)
                gTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTg                                                        >  1:5887638/1‑76 (MQ=255)
                  aaTATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTCCTTGAc                                                      >  1:1662030/1‑76 (MQ=255)
                  aaTATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGAc                                                      >  1:1880624/1‑76 (MQ=255)
                   atatACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGAcc                                                     <  2:3949004/76‑1 (MQ=255)
                   atatACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGAcc                                                     >  1:4487678/1‑76 (MQ=255)
                   atatACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGAcc                                                     <  1:3340453/76‑1 (MQ=255)
                    tataCACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCg                                                    >  2:5096034/1‑76 (MQ=255)
                    tataCACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCg                                                    >  1:6698327/1‑76 (MQ=255)
                     ataCACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCgg                                                   >  1:6335400/1‑76 (MQ=255)
                      taCACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCggg                                                  <  1:29436/76‑1 (MQ=255)
                      taCACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCggg                                                  >  1:6112804/1‑76 (MQ=255)
                       acacAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCgggg                                                 <  2:5228670/76‑1 (MQ=255)
                       acacAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCgggg                                                 >  1:6966141/1‑76 (MQ=255)
                       acacAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCgggg                                                 <  2:6923139/76‑1 (MQ=255)
                       acacAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCgggg                                                 <  2:4388894/76‑1 (MQ=255)
                       acacAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCgggg                                                 >  2:3643424/1‑76 (MQ=255)
                        cacaAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGt                                                <  2:6024883/76‑1 (MQ=255)
                        cacaAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGt                                                <  2:4294041/76‑1 (MQ=255)
                        cacaAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGt                                                >  1:2422769/1‑76 (MQ=255)
                        cacaAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGt                                                <  2:816512/76‑1 (MQ=255)
                        cacaAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGt                                                <  1:1862532/76‑1 (MQ=255)
                        cacaAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGt                                                >  1:4390016/1‑76 (MQ=255)
                         acaAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTg                                               <  2:869383/76‑1 (MQ=255)
                         acaAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTg                                               <  1:7109249/76‑1 (MQ=255)
                           aaaaTCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTCCTTGACCGGGTTGaa                                             <  2:1880605/76‑1 (MQ=25)
                           aaaaTCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGaa                                             >  2:4428220/1‑76 (MQ=255)
                           aaaaTCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGAGaa                                             >  2:3036415/1‑76 (MQ=255)
                            aaaTCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAAc                                            <  2:3526783/76‑1 (MQ=255)
                              aTCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGa                                          <  2:6710213/76‑1 (MQ=255)
                               tCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAg                                         >  1:7533285/1‑76 (MQ=255)
                               tCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAg                                         <  2:5797357/76‑1 (MQ=255)
                               tCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAg                                         <  2:5006026/76‑1 (MQ=255)
                                 aTTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAgcg                                       <  2:2327760/76‑1 (MQ=255)
                                 aTTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAgcg                                       <  1:7158352/76‑1 (MQ=255)
                                  ttCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAgcgc                                      <  2:2424428/76‑1 (MQ=255)
                                   tCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCa                                     <  1:7200830/76‑1 (MQ=255)
                                     aGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGc                                   >  1:3460426/1‑76 (MQ=255)
                                     aGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGc                                   >  1:596367/1‑76 (MQ=255)
                                     aGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGc                                   >  1:6371450/1‑76 (MQ=255)
                                      gATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGcc                                  >  1:2203651/1‑76 (MQ=255)
                                      gATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGcc                                  <  1:4452453/76‑1 (MQ=255)
                                       aTTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCa                                 >  2:2248591/1‑76 (MQ=255)
                                       aTTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCa                                 <  2:2374845/76‑1 (MQ=255)
                                       aTTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCa                                 >  1:1986709/1‑76 (MQ=255)
                                        ttGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCaa                                <  1:316582/76‑1 (MQ=255)
                                        ttGCAACAAGGCGTCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCaa                                <  2:6887615/76‑1 (MQ=38)
                                         tGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAAc                               >  1:1486185/1‑76 (MQ=255)
                                         tGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAAc                               >  1:6315621/1‑76 (MQ=255)
                                         tGCATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGCGAACGAGCGCAGCCAAc                               >  1:261304/1‑76 (MQ=255)
                                           cATCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGc                             >  1:6795389/1‑76 (MQ=255)
                                            aTCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACgaa                            >  2:2937821/1‑74 (MQ=255)
                                             tCAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCag                           >  2:6658992/1‑76 (MQ=255)
                                              cAAGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCaga                          >  2:4829211/1‑76 (MQ=255)
                                               aaGGCGGCAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCagag                         <  1:1578019/76‑1 (MQ=255)
                                                  gcggcAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCa                      <  2:6405979/76‑1 (MQ=255)
                                                  gcggcAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCa                      <  2:1559741/76‑1 (MQ=255)
                                                  gcggcAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCa                      <  2:3324675/76‑1 (MQ=255)
                                                   cggcAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAg                     <  1:405046/76‑1 (MQ=255)
                                                      cAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCt                  >  2:3430689/1‑76 (MQ=255)
                                                      cAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCt                  >  1:2516053/1‑76 (MQ=255)
                                                      cAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCt                  <  2:7639458/76‑1 (MQ=255)
                                                      cAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCt                  >  2:1693368/1‑76 (MQ=255)
                                                      cAAGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCt                  >  1:4780919/1‑76 (MQ=255)
                                                       aaGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTg                 <  2:2738448/76‑1 (MQ=255)
                                                       aaGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTg                 >  1:1065301/1‑76 (MQ=255)
                                                        aGTGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGa                <  2:4947868/76‑1 (MQ=255)
                                                          tGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAg              <  2:1120666/76‑1 (MQ=255)
                                                          tGAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAg              >  2:6968453/1‑76 (MQ=255)
                                                           gAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAgg             >  1:6762005/1‑76 (MQ=255)
                                                           gAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAgg             >  1:2814896/1‑76 (MQ=255)
                                                           gAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAgg             >  1:3857296/1‑76 (MQ=255)
                                                           gAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAgg             >  1:512295/1‑76 (MQ=255)
                                                           gAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAgg             >  2:5635609/1‑76 (MQ=255)
                                                           gAGCGAATCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAgg             >  1:872307/1‑76 (MQ=255)
                                                            aGCGAAACCCCGGAAGCGTACAGAAGAACTTGACCGGGATGACCGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGa            <  1:782250/76‑1 (MQ=255)
                                                                aaTCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGaa        >  1:4709786/1‑76 (MQ=255)
                                                                aaTCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGaa        <  2:34555/76‑1 (MQ=255)
                                                                aaTCCCCGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGaa        <  1:1304075/76‑1 (MQ=255)
                                                                   ccccGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAgtg     >  2:644682/1‑76 (MQ=255)
                                                                   ccccGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAgtg     >  1:3667216/1‑76 (MQ=255)
                                                                   ccccGGAAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAgtg     <  2:2563826/76‑1 (MQ=255)
                                                                      cGGAAGCGTACAGAAGTACTTGACCGGGGTGAGCGAGCCCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGtat  <  2:2891524/76‑1 (MQ=37)
                                                                         |                                                                        
GAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGATTGCATCAAGGCGGCAAGTGAGCGAATCCCCGGGAGCGTACAGAAGTACTTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTAT  >  REL606/4561073‑4561218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑