Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,125,004 A→T 7.7% D167V (GAT→GTT)  dam → DNA adenine methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,125,0040AT 7.7% 145.2 / 4.4 52D167V (GAT→GTT) damDNA adenine methylase
Reads supporting (aligned to +/- strand):  ref base A (21/27);  new base T (2/2);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

AGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGAC  >  W3110S.gb/4124939‑4125073
                                                                 |                                                                     
aGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCa                                                                  <  1:988715/71‑1 (MQ=255)
aGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCa                                                                  <  1:214166/71‑1 (MQ=255)
  cAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATg                                                                >  1:1422269/1‑71 (MQ=255)
    gagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgc                                                              >  1:970576/1‑71 (MQ=255)
    gagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgc                                                              >  1:122689/1‑71 (MQ=255)
    gagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgc                                                              <  1:247652/71‑1 (MQ=255)
    gagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgc                                                              >  1:1623611/1‑71 (MQ=255)
         gTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATg                                                                >  1:1442625/1‑64 (MQ=255)
            tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCa                                                       <  1:80863/70‑1 (MQ=255)
              aCTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCa                                                                  >  1:242685/1‑57 (MQ=255)
               cTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                    >  1:944759/1‑70 (MQ=255)
                 tCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCa                                                       <  1:614080/65‑1 (MQ=255)
                  cGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatga                                                  <  1:866728/69‑1 (MQ=255)
                  cGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGatgatg                                                >  1:811378/1‑71 (MQ=255)
                  cGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGatgatg                                                >  1:1551104/1‑71 (MQ=255)
                   gCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGc                                               <  1:525201/71‑1 (MQ=255)
                   gCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGc                                               <  1:364511/71‑1 (MQ=255)
                   gCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGc                                               <  1:1026484/71‑1 (MQ=255)
                    cTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCa                                              <  1:492047/71‑1 (MQ=255)
                       aaaaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAGATGATGCATcc                                           >  1:1647/1‑71 (MQ=255)
                         aaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATgg                                                               >  1:1768586/1‑49 (MQ=255)
                         aaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgt                                         >  1:911472/1‑71 (MQ=255)
                             cgcAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGc                                                                   <  1:1075103/41‑1 (MQ=255)
                             cgcAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGc                                                                   <  1:1456160/41‑1 (MQ=255)
                                   aTGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgtc                                        <  1:1291729/62‑1 (MQ=255)
                                   aTGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgtc                                        <  1:131841/62‑1 (MQ=255)
                                     gCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                    <  1:533654/48‑1 (MQ=255)
                                      ccTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGa                             <  1:939547/70‑1 (MQ=255)
                                      ccTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGa                             <  1:1518072/70‑1 (MQ=255)
                                        tttttCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAgat                                                    <  1:97811/45‑1 (MQ=255)
                                        tttttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                    <  1:863045/45‑1 (MQ=255)
                                        tttttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                    <  1:1744619/45‑1 (MQ=255)
                                        tttttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcg                                          <  1:1232236/55‑1 (MQ=255)
                                        tttttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATcc                          >  1:1648617/1‑71 (MQ=255)
                                          tttCTATTGTGAGTCTTACGCCGTTAGCATGgcgcgcg                                                         >  1:1348668/1‑38 (MQ=255)
                                          tttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTa                                   >  1:396103/1‑60 (MQ=255)
                                           ttCTATTGTGAGTCTTACGCCGATAGCATGgcgcgcg                                                         <  1:681088/37‑1 (MQ=255)
                                            tCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTa                                   >  1:254677/1‑58 (MQ=255)
                                            tCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTa                                   >  1:548531/1‑58 (MQ=255)
                                            tCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCt                      >  1:243929/1‑71 (MQ=255)
                                               aTTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAg                                                      <  1:904902/36‑1 (MQ=255)
                                                ttGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCtt                     <  1:981549/68‑1 (MQ=255)
                                                  gtgAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                    <  1:1510567/35‑1 (MQ=255)
                                                  gtgAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATcc                          <  1:611385/61‑1 (MQ=255)
                                                  gtgAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATcc                          <  1:826251/61‑1 (MQ=255)
                                                  gtgAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCa                >  1:1712052/1‑71 (MQ=255)
                                                     aGTCTTACGCCGATAGCATGGCGCGCGCAGATGATTCATCCGTCGTCTATTGCGATCCGCCTTATGCACCg             <  1:876344/71‑1 (MQ=255)
                                                       tCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGcc                       <  1:403474/59‑1 (MQ=255)
                                                         ttACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGt         >  1:508485/1‑71 (MQ=255)
                                                           aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCt       >  1:1007686/1‑71 (MQ=255)
                                                                gATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGAt                            >  1:1157685/1‑45 (MQ=255)
                                                                gATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGAc  >  1:846156/1‑71 (MQ=255)
                                                                 |                                                                     
AGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGAC  >  W3110S.gb/4124939‑4125073

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