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,1630GA100.0% 102.7 / NA 27L111L (CTG→TTG) ogtmethylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/18);  total (9/18)

GGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTA  >  NZ_CP009273/1394076‑1394245
                                                                                       |                                                                                  
ggCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAAtt                                                                                  >  2:1173/1‑90 (MQ=255)
ggCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAAtt                                                                                  <  1:735084/90‑1 (MQ=255)
ggCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAAtt                                                                                  <  1:41823/90‑1 (MQ=255)
          gacgaTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATg                                                                        <  1:698896/90‑1 (MQ=255)
             gaTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCAt                                                                     <  2:389157/90‑1 (MQ=255)
               tGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTa                                                                   >  1:308316/1‑90 (MQ=255)
               tGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTa                                                                   >  1:141888/1‑90 (MQ=255)
                gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                                  >  2:671606/1‑90 (MQ=255)
                 cTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAcc                                                                 <  2:166847/90‑1 (MQ=255)
                   gATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTg                                                               <  1:325783/90‑1 (MQ=255)
                        gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCg                                                           <  1:416576/89‑1 (MQ=255)
                        gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCg                                                           >  2:416576/1‑89 (MQ=255)
                        gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGc                                                          >  2:726673/1‑90 (MQ=255)
                         aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACTTGCCCGCa                                                         <  1:890824/90‑1 (MQ=255)
                         aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCa                                                         <  2:4407/90‑1 (MQ=255)
                          ttCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAg                                                        <  1:767120/90‑1 (MQ=255)
                          ttCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAg                                                        >  2:359326/1‑90 (MQ=255)
                            cGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAggg                                                      <  1:619447/90‑1 (MQ=255)
                              aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGa                                                    <  1:624320/90‑1 (MQ=255)
                                 cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                 <  2:18339/90‑1 (MQ=255)
                                      gCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGt                                            <  2:814140/90‑1 (MQ=255)
                                          cACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAgtg                                        <  2:197659/90‑1 (MQ=255)
                                          cACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAgtg                                        <  1:500504/90‑1 (MQ=255)
                                          cACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAgtg                                        <  1:373121/90‑1 (MQ=255)
                                                                      gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt            >  2:858600/1‑90 (MQ=255)
                                                                               cTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGt   <  2:597804/90‑1 (MQ=255)
                                                                                tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTa  >  2:602057/1‑90 (MQ=255)
                                                                                       |                                                                                  
GGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTA  >  NZ_CP009273/1394076‑1394245

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