Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,535,089 (A)6→5 coding (9/597 nt) yhhF → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,535,0840A.94.9% 172.1 / 3.8 39coding (4/597 nt)yhhFpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (13/24);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

AGTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGG  >  REL606/3535010‑3535156
                                                                          |                                                                        
aGTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaa                                                                         >  2:600973/1‑76 (MQ=255)
aGTCTATCACTTAACCTAATTCACATCACCGCCAGCAAGTATGTGTTATCTGGCGGATTGANCAATTTATCATGaa                                                                         >  1:1144387/1‑76 (MQ=255)
 gTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaa                                                                        <  2:1224731/76‑1 (MQ=255)
 gTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaa                                                                        >  1:937743/1‑76 (MQ=255)
 gTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaa                                                                        >  1:352215/1‑76 (MQ=255)
    tATCACTTAACTTAATTCACATCACTGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaac                                                                     >  2:382944/1‑75 (MQ=255)
    tATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaac                                                                     <  1:441032/76‑2 (MQ=255)
     aTCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaacc                                                                    <  2:391716/76‑3 (MQ=255)
     aTCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaacc                                                                    >  2:169113/1‑74 (MQ=255)
      tCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaaccg                                                                   >  2:532247/1‑73 (MQ=255)
      tCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaaccg                                                                   >  2:795113/1‑73 (MQ=255)
         cttaacttaaTTCACATCACCGCCTGGAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAAt                                                               <  2:77038/76‑1 (MQ=255)
           taacttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCa                                                             <  1:611612/76‑1 (MQ=255)
            aacttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCAt                                                            <  2:1205989/76‑1 (MQ=255)
            aacttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCAt                                                            >  1:563858/1‑76 (MQ=255)
              cttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTc                                                          <  1:1061965/76‑1 (MQ=255)
              cttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTc                                                          >  2:631831/1‑76 (MQ=255)
               ttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTcc                                                         <  1:539490/76‑1 (MQ=255)
                 aaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCgg                                                       >  1:252262/1‑76 (MQ=255)
                   ttCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCa                                                     <  1:650592/76‑1 (MQ=255)
                   ttCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCa                                                     <  1:650604/76‑1 (MQ=255)
                    tCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAg                                                    >  1:430344/1‑76 (MQ=255)
                       caTCACCGCCTGCANGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCgg                                                 >  1:3186/1‑76 (MQ=255)
                        aTCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGc                                                <  2:657742/76‑1 (MQ=255)
                                   cAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCat                                     >  2:735498/1‑76 (MQ=255)
                                   cAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCat                                     >  2:156883/1‑76 (MQ=255)
                                     aGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCatta                                   <  1:489840/76‑1 (MQ=255)
                                     aGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCatta                                   <  1:489828/76‑1 (MQ=255)
                                      gTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCattat                                  >  2:181750/1‑76 (MQ=255)
                                         tgtgTTATCTGGCGGATTGAGCAATTGATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTANTgg                               <  1:1204118/76‑1 (MQ=255)
                                          ntgtTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGc                              >  1:957496/2‑76 (MQ=255)
                                          gtgtTATCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGc                              <  1:944812/76‑1 (MQ=255)
                                             ttATCTGGCGGATTGAGCAATTTATCATNTAAAACCGAAACATTCCGGGANCGTCCAAATCCGCATTATTGGCggg                           >  1:195049/1‑76 (MQ=255)
                                               aTCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCa                         <  1:344517/76‑1 (MQ=255)
                                                tCTGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAg                        <  1:651686/76‑1 (MQ=255)
                                                  tGGCGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTg                      <  2:1144387/76‑1 (MQ=255)
                                                  tGGCGGANGGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAggg                      >  2:1075336/1‑74 (MQ=255)
                                                   ggCGGATTGAGCACTTTATCNTGAAAAACCGAATCANTCCGGCAGCGGCCAACTCCGCATTATTGGCGGGCAGTgg                     <  2:1336303/76‑1 (MQ=255)
                                                     cGGATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCg                   >  1:764349/1‑76 (MQ=255)
                                                     cGGATTGAGCAATTTATCATGAAAAACCGAAACATACCGGCAGCGGCCAAATCCGCATTATTGGCGGGAAGTGGCg                   <  2:783914/76‑1 (MQ=255)
                                                      ggATTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGa                  <  2:443742/76‑1 (MQ=255)
                                                        aTTGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAgg                <  2:181987/76‑1 (MQ=255)
                                                          tGAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGcc              <  2:550342/76‑1 (MQ=255)
                                                           gAGCAATTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCg             <  1:709653/76‑1 (MQ=255)
                                                               aaTTTATCATGAAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTaa          <  2:1269921/76‑1 (MQ=255)
                                                               aaTTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGGGGCGAGGCCGCaaa         >  2:876147/1‑76 (MQ=255)
                                                                aTTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAAc        <  1:600896/76‑1 (MQ=255)
                                                                aTTTATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAAc        <  1:511523/76‑1 (MQ=255)
                                                                  ttATCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGGGGCCGCAAAcac      >  2:240405/1‑74 (MQ=255)
                                                                     tCATGAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCg   <  2:352215/76‑1 (MQ=255)
                                                                      cATGAAAAAACGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGGCCTAAACTccccg  >  2:451021/1‑74 (MQ=255)
                                                                          |                                                                        
AGTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGG  >  REL606/3535010‑3535156

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