Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,394,163 G→A L111L (CTG→TTG)  ogt ← methylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,394,1630GA97.1% 127.9 / NA 35L111L (CTG→TTG) ogtmethylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/13);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAG  >  NZ_CP009273/1394079‑1394246
                                                                                    |                                                                                   
aaGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGc                                                                                >  2:373345/1‑90 (MQ=255)
    tACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTa                                                                            <  2:435908/90‑1 (MQ=255)
            tGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTa                                                                    >  1:397776/1‑90 (MQ=255)
             gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                                   >  2:329255/1‑90 (MQ=255)
             gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                                   >  2:675301/1‑90 (MQ=255)
             gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                                   >  1:391904/1‑90 (MQ=255)
                  tGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGcc                                                              <  1:1003665/90‑1 (MQ=255)
                  tGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGcc                                                              >  2:583887/1‑90 (MQ=255)
                    ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCg                                                            <  1:472975/90‑1 (MQ=255)
                    ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCg                                                            <  1:24328/90‑1 (MQ=255)
                       ttCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGACCGAAATGCATTACCTGCCCGCAg                                                         >  1:715659/1‑90 (MQ=255)
                          gATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAgggg                                                      >  1:1178894/1‑90 (MQ=255)
                          gATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAgggg                                                      >  2:1042962/1‑90 (MQ=255)
                             ccGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATa                                                   <  1:1020847/90‑1 (MQ=255)
                               gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                  >  2:568301/1‑89 (MQ=255)
                               gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                  <  1:568301/89‑1 (MQ=255)
                                   gCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGt                                             <  1:710417/90‑1 (MQ=255)
                                       cACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAgtg                                         <  1:291191/90‑1 (MQ=255)
                                                                  ggCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCg              <  2:499854/90‑1 (MQ=255)
                                                                   gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                  <  2:885310/53‑1 (MQ=255)
                                                                   gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                  >  1:885310/1‑53 (MQ=255)
                                                                   gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt             >  1:977224/1‑90 (MQ=255)
                                                                   gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt             >  2:159177/1‑90 (MQ=255)
                                                                   gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt             >  2:442806/1‑90 (MQ=255)
                                                                   gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt             >  1:566246/1‑90 (MQ=255)
                                                                   gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt             >  1:1000441/1‑90 (MQ=255)
                                                                     ccAATTGCCCAGCCAGTTCGCCGTAATGCATTACCTGCCCGCAGGGGGTAGTGCGTAGTGGTTTCCAGACTTCGCGCTGAAATGGGGTcc           <  1:1037588/90‑1 (MQ=255)
                                                                          tGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTccccccc      <  1:95568/90‑1 (MQ=255)
                                                                             tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcccc         <  2:533594/84‑1 (MQ=255)
                                                                             tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcccc         >  1:533594/1‑84 (MQ=255)
                                                                             tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTa   >  2:1029325/1‑90 (MQ=255)
                                                                             tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTa   >  2:189928/1‑90 (MQ=255)
                                                                             tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTa   >  1:760844/1‑90 (MQ=255)
                                                                             tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTa   >  1:464598/1‑90 (MQ=255)
                                                                              cAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAg  <  2:397776/90‑1 (MQ=255)
                                                                                    |                                                                                   
AAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAG  >  NZ_CP009273/1394079‑1394246

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