Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 539,303 A→G C97R (TGC→CGC)  ECB_00490 ← putative DNA N‑6‑adenine‑methyltransferase of bacteriophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606539,3030AG98.2% 659.4 / 0.2 166C97R (TGC→CGC) ECB_00490putative DNA N‑6‑adenine‑methyltransferase of bacteriophage
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (85/78);  total (85/81)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 173 total aligned reads displayed.

ATGCAGCCTCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCAGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACGTGTTTCCCGGT  >  REL606/539228‑539375
                                                                           |                                                                         
    aGCCTCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGCGGCGgtgt                                                                       >  1:1015858/1‑76 (MQ=255)
      ccTCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGgtgtgt                                                                     >  1:5378456/1‑76 (MQ=255)
      ccTCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGgtgtgt                                                                     >  2:304695/1‑76 (MQ=255)
      ccTCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGgggtgt                                                                     >  1:3722271/1‑76 (MQ=255)
       cTCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTa                                                                   <  1:1679991/76‑1 (MQ=255)
        tCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTAc                                                                  <  1:1209708/76‑1 (MQ=255)
        tCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTAc                                                                  <  1:3615366/76‑1 (MQ=255)
        tCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTAc                                                                  <  2:510037/76‑1 (MQ=255)
         cccGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACg                                                                 >  2:985209/1‑76 (MQ=255)
         cccGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACg                                                                 >  2:23538/1‑76 (MQ=255)
         cccGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACg                                                                 >  1:1418917/1‑76 (MQ=255)
         cccGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACg                                                                 <  2:5998726/76‑1 (MQ=255)
          ccGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACgg                                                                >  1:748695/1‑76 (MQ=255)
          ccGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGGGTGTACgg                                                                >  2:2843233/1‑76 (MQ=255)
             aCGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACggggg                                                             >  1:3444114/1‑76 (MQ=255)
             aCGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACggggg                                                             >  1:1667189/1‑76 (MQ=255)
              cGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGt                                                            <  2:1929662/76‑1 (MQ=255)
               gggCAAGAATCCTTGAGCCGAACGCCATCACCTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTa                                                           <  1:813624/76‑1 (MQ=255)
                ggCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTAt                                                          <  1:2143609/76‑1 (MQ=255)
                 gCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATg                                                         >  2:3799409/1‑76 (MQ=255)
                 gCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATg                                                         >  2:2078530/1‑76 (MQ=255)
                 gCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATAAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATg                                                         >  1:2486307/1‑76 (MQ=255)
                   aaGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATggg                                                       >  2:3866418/1‑76 (MQ=255)
                    aGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGt                                                      <  1:3210172/76‑1 (MQ=255)
                     gAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTg                                                     <  2:1495869/76‑1 (MQ=255)
                      aaTCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGGGc                                                    >  2:2892007/1‑76 (MQ=255)
                       aTCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTggc                                                   >  1:5893717/1‑74 (MQ=255)
                       aTCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGcc                                                   <  2:5199949/76‑1 (MQ=255)
                       aTCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGGGTGTACGGGGGGATGGGGGcc                                                   >  1:3220805/1‑76 (MQ=255)
                        tCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCa                                                  <  2:2616044/76‑1 (MQ=255)
                        tCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGGGTGGACGGGGGGATGGGTGCCa                                                  >  1:1840771/1‑76 (MQ=255)
                         ccTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAg                                                 >  1:449509/1‑76 (MQ=255)
                         ccTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAg                                                 <  1:4962800/76‑1 (MQ=255)
                         ccTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGGGCCAg                                                 >  2:5827217/1‑76 (MQ=255)
                         ccTTGAGCCAAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAg                                                 <  1:4020319/76‑1 (MQ=255)
                         ccTGGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAg                                                 <  2:63380/76‑1 (MQ=255)
                          cTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGa                                                <  1:691247/76‑1 (MQ=255)
                          cTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGGGTGGACGGGGGAATGGGTGCCAGa                                                >  1:923396/1‑76 (MQ=255)
                           ttGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGAt                                               <  2:1076543/76‑1 (MQ=255)
                           ttGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGAt                                               <  1:1300904/76‑1 (MQ=255)
                             gAGCCGAACAGCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATaa                                             <  2:1503565/76‑1 (MQ=255)
                              aGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAAt                                            <  2:3444826/76‑1 (MQ=255)
                              aGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAAt                                            <  1:884387/76‑1 (MQ=255)
                              aGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAAt                                            >  1:5047693/1‑76 (MQ=255)
                                 cGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATTAg                                         >  2:5460118/1‑76 (MQ=255)
                                 cGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAg                                         <  1:6295895/76‑1 (MQ=255)
                                 cGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGGGGGTACGGGGGTATGGGGGCCAACTAATGAg                                         >  2:4690006/1‑76 (MQ=255)
                                   aaCGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCa                                       <  2:4361067/76‑1 (MQ=255)
                                   aaCGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAATGAGcc                                       >  2:2354224/1‑75 (MQ=255)
                                   aaCGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGAGGGGATGGGTTCCAGATAAAGAGCa                                       >  2:162259/1‑76 (MQ=255)
                                    aCGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCAt                                      >  1:2196226/1‑76 (MQ=255)
                                    aCGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCAt                                      <  1:4113602/76‑1 (MQ=255)
                                      gCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCg                                    <  1:3804261/76‑1 (MQ=255)
                                      gCCATCAACTCTCCACGATCAACGGTCGTAAAGCGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCg                                    >  1:5127554/1‑76 (MQ=255)
                                       ccATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGa                                   <  2:5716197/76‑1 (MQ=255)
                                       ccATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGa                                   >  2:1251662/1‑76 (MQ=255)
                                       ccATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAATGAGCATCGa                                   >  1:502113/1‑76 (MQ=255)
                                       ccATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGGGTACGGGGGGATGGGTGCCAGGTAATGAGCACCgg                                   >  2:1397644/1‑75 (MQ=255)
                                        cATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAg                                  <  2:3674678/76‑1 (MQ=255)
                                        cATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAg                                  <  2:6159046/76‑1 (MQ=255)
                                        cATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAg                                  <  1:5647194/76‑1 (MQ=255)
                                        cATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAATGAGCATCGAg                                  >  2:857337/1‑76 (MQ=255)
                                        cATCAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAATGAGCATCGAg                                  >  2:4756947/1‑76 (MQ=255)
                                           cAACTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCt                               >  1:6270235/1‑76 (MQ=255)
                                           cAACTATCCACGATCAACGGTCGTAAAGTGGCGGTGGGTACGGGGGTATGGGGGCCAGATGATAGGCATCGAGCCt                               >  2:6286133/1‑76 (MQ=255)
                                            aaCTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTAACGGGGGCAGGGGGCCAGATAATCAACATCCAGcgct                              >  2:2441984/1‑73 (MQ=255)
                                            aaCTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGGCAGATAATAAGCACCGAGCCtt                              >  2:3730347/1‑76 (MQ=255)
                                            aaCTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCtt                              >  1:2314745/1‑76 (MQ=255)
                                            aaCTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCtt                              >  2:226338/1‑76 (MQ=255)
                                             aCTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCttt                             <  1:3352235/76‑1 (MQ=255)
                                             aCTCTCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCttt                             <  2:3313816/76‑1 (MQ=255)
                                              ctctCCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTtta                            <  2:2516637/76‑1 (MQ=255)
                                               tctcCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTttat                           >  1:532419/1‑76 (MQ=255)
                                                ctcCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTttatt                          <  1:1480874/76‑1 (MQ=255)
                                                ctcCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGGGCCAGATAATGAGCATCGAGCCTttatt                          >  1:113933/1‑76 (MQ=255)
                                                ctcCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAATGAGCATCGAGCCTttatt                          >  1:3379710/1‑76 (MQ=255)
                                                 tcCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTttatta                         <  2:417443/76‑1 (MQ=255)
                                                 tcCACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTttatta                         >  1:4235356/1‑76 (MQ=255)
                                                  ccACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTtattat                        >  1:4952996/1‑76 (MQ=255)
                                                  ccACGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTtattat                        <  2:1895425/76‑1 (MQ=255)
                                                     cGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTc                     <  1:1615699/76‑1 (MQ=255)
                                                     cGATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGGGCCAGATAATGAGCATCGAGCCTTTATTATTTc                     >  1:5053062/1‑76 (MQ=255)
                                                      gATCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTcc                    >  1:3749123/1‑76 (MQ=255)
                                                       aTCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCa                   <  1:2614277/76‑1 (MQ=255)
                                                       aTCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCa                   <  1:981726/76‑1 (MQ=255)
                                                       aTCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCa                   <  2:4620678/76‑1 (MQ=255)
                                                       aTCAACGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGGGCCAGATAATGAGCATCGAGCCTTTATTATTTTCa                   >  1:2894550/1‑76 (MQ=255)
                                                          aaCGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTg                <  1:595464/76‑1 (MQ=255)
                                                           aCGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTTCACTGa               >  2:4076460/1‑76 (MQ=255)
                                                           aCGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGa               <  2:4291151/76‑1 (MQ=255)
                                                            cGGTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGAc              >  1:3193022/1‑76 (MQ=255)
                                                              gTCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCCTTATTATTTCCACTGACgt            >  2:5252899/1‑76 (MQ=255)
                                                              gTCGTAAAGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACgt            >  2:905515/1‑76 (MQ=255)
                                                               tCGTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACgtg           <  2:3417665/76‑1 (MQ=255)
                                                                 gTAAAGTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACGTGtt         >  2:5676460/1‑76 (MQ=255)
                                                                 gTAAAGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACGTGtt         >  2:964516/1‑76 (MQ=255)
                                                                   aaaGTGGCGGTGTGTACGGGGGGATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACGTGTTTc       >  2:6159913/1‑76 (MQ=255)
                                                                      gTGGCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACGTGTTTCCCg    >  1:3207581/1‑76 (MQ=255)
                                                                        ggCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACGTGTTTCCCGGt  <  1:3719058/76‑1 (MQ=255)
                                                                        ggCGGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACGTGTTTCCCGGt  >  1:2677930/1‑76 (MQ=255)
                                                                           |                                                                         
ATGCAGCCTCCCGACGGGCAAGAATCCTTGAGCCGAACGCCATCAACTCTCCACGATCAACGGTCGTAAAGTGGCAGTGTGTACGGGGGTATGGGTGCCAGATAATGAGCATCGAGCCTTTATTATTTCCACTGACGTGTTTCCCGGT  >  REL606/539228‑539375

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