Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,340,143 G→A P200P (CCG→CCA yhdJ → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,340,1430GA98.4% 731.4 / 0.9 185P200P (CCG→CCAyhdJpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (98/84);  total (100/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

Only 100 of 192 total aligned reads displayed.

GGGTAACGTCTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCCGGT  >  REL606/3340071‑3340215
                                                                        |                                                                        
gggTANNGTCTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAANACCAGaa                                                                       >  2:1067833/1‑76 (MQ=255)
 ggTACCGTCTGGGATTTCCCGCGCGTGCGTTATTTAATGGATAAATATGAGAACCACCCGACGCAAAAACCAGAAg                                                                      <  1:1772848/76‑1 (MQ=255)
 ggTAACGTCTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAg                                                                      >  2:1871400/1‑76 (MQ=255)
   tAACGTCTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGcc                                                                    >  2:37350/1‑76 (MQ=255)
    aaCGTCTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCt                                                                   >  1:1007738/1‑76 (MQ=255)
     aCGTCTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCtt                                                                  <  2:4691802/76‑1 (MQ=255)
      cGTCTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTa                                                                 >  2:5236335/1‑76 (MQ=255)
      cGTCTGGGATTTTCCGCGCGTGCGTTATTTAATGAAACAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTa                                                                 <  1:2761842/76‑1 (MQ=255)
        tcggggATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACt                                                               <  2:2571208/73‑1 (MQ=255)
        tCTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACt                                                               >  1:2540784/1‑76 (MQ=255)
         cTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTg                                                              <  1:5598796/76‑1 (MQ=255)
         cTGGGATTTTCCGCGCGTGCGTTATTTAAGGGATAAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTg                                                              <  2:1217401/76‑1 (MQ=255)
           ggggTTTTCCGCGCGTGCGTTATTTAAGGGATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGaa                                                            <  1:1679436/76‑1 (MQ=255)
           gggATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGaa                                                            >  2:1722826/1‑76 (MQ=255)
             gATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAAc                                                          <  2:2972888/76‑1 (MQ=255)
              aTTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACg                                                         >  1:5890930/1‑76 (MQ=255)
              aTTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACg                                                         <  1:4579906/76‑1 (MQ=255)
              aTTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAAAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACg                                                         <  2:4738694/76‑1 (MQ=255)
              aTTTTCCGCACGTGCTTTATTTAATGTAGAAATAGGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACg                                                         <  1:3138799/76‑1 (MQ=255)
               ttttCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGc                                                        >  1:613667/1‑76 (MQ=255)
               ttttCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGc                                                        >  1:297386/1‑76 (MQ=255)
                 ttCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCat                                                      >  1:3542189/1‑76 (MQ=255)
                   ccGCGCGTGCGTTATTTAATGGATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCatta                                                    <  2:3966330/76‑1 (MQ=255)
                   ccGCGCGTGCGTTATTTAAGGGATGAATATGAAACCCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCatta                                                    <  2:2790998/76‑1 (MQ=255)
                    cgcgcgTGCGTTATTTAATGGATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCattat                                                   <  2:3988662/76‑1 (MQ=255)
                    cgcgcgTGCGTTATTTAATGGATAAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCattat                                                   <  1:5757404/76‑1 (MQ=255)
                        cgtgcgtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATtctc                                               <  2:5649853/76‑1 (MQ=255)
                        cgtgcgtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATtctc                                               >  2:2181675/1‑76 (MQ=255)
                         gtgcgtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCg                                              >  2:1309987/1‑76 (MQ=255)
                         gtgcgtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCg                                              >  1:5609718/1‑76 (MQ=255)
                          tgcgtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGc                                             >  1:5237302/1‑76 (MQ=255)
                            cgtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCct                                           >  1:3035628/1‑76 (MQ=255)
                            cgtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCct                                           >  1:2903552/1‑76 (MQ=255)
                            cgtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCct                                           >  1:4485157/1‑76 (MQ=255)
                            cgtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCct                                           >  2:3801232/1‑76 (MQ=255)
                             gtTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCctc                                          >  2:5863464/1‑76 (MQ=255)
                              ttATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCctct                                         >  1:471888/1‑76 (MQ=255)
                               tATTTAATTGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCtt                                        >  1:1820212/1‑76 (MQ=255)
                               tATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCtt                                        >  1:862414/1‑76 (MQ=255)
                               tATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCtt                                        >  2:1642055/1‑76 (MQ=255)
                                aTTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTc                                       >  2:3068879/1‑76 (MQ=255)
                                aTTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTc                                       >  1:3443791/1‑76 (MQ=255)
                                aTTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTc                                       >  2:4655548/1‑76 (MQ=255)
                                aTTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTc                                       >  2:1713024/1‑76 (MQ=255)
                                aTTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTc                                       >  2:4133560/1‑76 (MQ=255)
                                aTTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTc                                       >  1:428039/1‑76 (MQ=255)
                                aTTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTc                                       >  1:522876/1‑76 (MQ=255)
                                 tttAATGGATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTcc                                      <  2:4453426/76‑1 (MQ=255)
                                 tttAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTcc                                      >  1:1375716/1‑76 (MQ=255)
                                 tttAAGGGATGAATATGAGACCCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTcc                                      <  2:320050/76‑1 (MQ=255)
                                 tttAAGGGATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTcc                                      <  2:5989434/76‑1 (MQ=255)
                                 tttAAGGGAGGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTcc                                      <  2:644113/76‑1 (MQ=255)
                                 ttaaaTGGATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTcc                                      <  2:741356/73‑1 (MQ=255)
                                  ttAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCa                                     >  1:5670646/1‑76 (MQ=255)
                                  ttAATGGATGAATAGGAGACCCCCCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCa                                     <  1:1389119/76‑1 (MQ=255)
                                  ttAAGGGATGAATATGAAAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCa                                     <  2:1367877/76‑1 (MQ=255)
                                   tAATGGATGAATATGAGACCCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCaa                                    <  2:2024814/76‑1 (MQ=255)
                                   tAATGGATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCaa                                    <  1:2687573/76‑1 (MQ=255)
                                   tAATGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCaa                                    >  1:4274284/1‑76 (MQ=255)
                                    aaTGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAAc                                   >  2:4809083/1‑76 (MQ=255)
                                    aaTGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAAc                                   >  2:2895712/1‑76 (MQ=255)
                                     aTGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGACACGCATTATTCTCGCCTCTTCCAAcc                                  >  2:1413223/1‑76 (MQ=255)
                                     aTGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAAcc                                  >  2:5730484/1‑76 (MQ=255)
                                      tGGATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAAccc                                 >  1:3774741/1‑76 (MQ=255)
                                      tGGATGAATATGAAACCCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAAccc                                 <  2:3928607/76‑1 (MQ=255)
                                       ggATGAATATGAGAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCa                                <  2:594485/76‑1 (MQ=255)
                                       ggATGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCa                                >  1:1584278/1‑76 (MQ=255)
                                        gaggAATAGGAAAACCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAg                               <  2:1954416/73‑1 (MQ=255)
                                        gATGAATAGGAGACCCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAg                               <  1:733359/76‑1 (MQ=255)
                                         aTGAATATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAgg                              >  2:4060594/1‑76 (MQ=255)
                                            aaTATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGa                           >  2:618459/1‑76 (MQ=255)
                                            aaTATGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGa                           >  1:962800/1‑76 (MQ=255)
                                             atagGAGACCCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGat                          <  1:2133831/76‑1 (MQ=255)
                                              tatGAGACCCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGata                         <  2:2076305/76‑1 (MQ=255)
                                              tatGAGAACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGata                         >  2:1503452/1‑76 (MQ=255)
                                               atGAGACCCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGatat                        <  2:4341459/76‑1 (MQ=255)
                                               aggagaCCCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGatat                        <  1:3732986/74‑1 (MQ=255)
                                                gaagaCCCCCCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATc                       <  1:537273/74‑1 (MQ=255)
                                                  agaACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGt                     >  1:5107727/1‑76 (MQ=255)
                                                  agaACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGt                     >  2:2322684/1‑76 (MQ=255)
                                                   gaACCACCCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGtt                    >  1:5055866/1‑76 (MQ=255)
                                                     acccccCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTtct                  <  1:1151676/76‑1 (MQ=255)
                                                     acccccCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTtct                  <  1:1433375/76‑1 (MQ=255)
                                                     accaccCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTtct                  >  2:5530038/1‑76 (MQ=255)
                                                      ccccccGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTtctc                 <  1:4025952/73‑1 (MQ=255)
                                                      ccaccCGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTtctc                 >  1:4740679/1‑76 (MQ=255)
                                                       cccccGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCg                <  2:5370004/74‑1 (MQ=255)
                                                         cccGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGAc              >  1:394074/1‑76 (MQ=255)
                                                           cGACGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGAccc            <  1:4923358/76‑1 (MQ=255)
                                                              cGCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGtt         <  2:4938271/76‑1 (MQ=255)
                                                               gCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGttt        <  2:3396841/76‑1 (MQ=255)
                                                               gCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGttt        >  1:905648/1‑76 (MQ=255)
                                                               gCAAAAACCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGttt        <  2:4246023/76‑1 (MQ=255)
                                                                 aaaaaCCAGAAGCCTTACTAAAACGAATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGc      <  1:4216483/76‑1 (MQ=255)
                                                                  aaaaCCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGcc     >  2:255732/1‑76 (MQ=255)
                                                                   aaaCCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCCg    <  2:1268125/76‑1 (MQ=255)
                                                                   aaaCCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCCg    >  2:4921533/1‑76 (MQ=255)
                                                                   aaaCCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCCg    >  2:4423265/1‑76 (MQ=255)
                                                                     ctgggAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCCGGt  <  1:805121/73‑1 (MQ=255)
                                                                     aCCAGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCCGGt  <  1:2430464/76‑1 (MQ=255)
                                                                        |                                                                        
GGGTAACGTCTGGGATTTTCCGCGCGTGCGTTATTTAATGGATGAATATGAGAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCCGGT  >  REL606/3340071‑3340215

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