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.3% 136.1 / ‑4.4 37L111L (CTG→TTG) ogtmethylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/19);  total (17/20)
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.21e-01

AGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGC  >  NZ_CP009273/1394080‑1394247
                                                                                   |                                                                                    
aGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGcc                                                                                >  1:1375702/1‑90 (MQ=255)
   tACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTa                                                                             <  1:1072251/90‑1 (MQ=255)
   tACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTa                                                                             <  2:1061655/90‑1 (MQ=255)
          aTGCTGATGGGATTTGATCCGTTTGCCGCACCAAGGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCGGTAATGCAtt                                                                      >  2:767365/1‑90 (MQ=255)
          aTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTa                                                                             >  2:358339/1‑83 (MQ=255)
          aTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTa                                                                             <  1:358339/83‑1 (MQ=255)
          aTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCAtt                                                                      >  1:1022633/1‑90 (MQ=255)
           tGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTa                                                                     >  2:1566056/1‑90 (MQ=255)
            gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                                    >  2:30206/1‑90 (MQ=255)
                   ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCg                                                             <  1:271522/90‑1 (MQ=255)
                     aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCa                                                           <  2:1272807/90‑1 (MQ=255)
                      ttCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAg                                                          <  1:87864/90‑1 (MQ=255)
                      ttCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAg                                                          <  2:1375702/90‑1 (MQ=255)
                         gATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAgggg                                                       >  1:412149/1‑90 (MQ=255)
                          aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGa                                                      <  1:30206/90‑1 (MQ=255)
                             cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAg                                                   <  1:220592/90‑1 (MQ=255)
                              gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGt                                                  >  2:523480/1‑90 (MQ=255)
                              gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGt                                                  >  2:195444/1‑90 (MQ=255)
                                  gCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGt                                              <  2:811826/90‑1 (MQ=255)
                                    cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAg                                            <  1:523480/90‑1 (MQ=255)
                                           aCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGGGTTTTc                                     >  2:431388/1‑90 (MQ=255)
                                                aCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCTTTACCTg                                                                 >  1:65874/1‑57 (MQ=255)
                                                aCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCTTTACCTg                                                                 <  2:65874/57‑1 (MQ=255)
                                                aCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGAc                                >  1:1432507/1‑90 (MQ=255)
                                                aCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGAc                                >  1:621097/1‑90 (MQ=255)
                                                  gcgccgcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACtt                              >  1:242056/1‑90 (MQ=255)
                                                     ccgcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTcgc                           <  2:1530913/90‑1 (MQ=255)
                                                     ccgcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTcgc                           <  2:328889/90‑1 (MQ=255)
                                                     ccgcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTcgc                           <  2:103832/90‑1 (MQ=255)
                                                       gcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCgcgc                         <  2:380286/90‑1 (MQ=255)
                                                       gcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCgcgc                         <  1:813345/90‑1 (MQ=255)
                                                                  gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt              >  1:687283/1‑90 (MQ=255)
                                                                  gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt              >  2:861714/1‑90 (MQ=255)
                                                                           cTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGt     <  2:242056/90‑1 (MQ=255)
                                                                           cTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGt     <  2:1338847/90‑1 (MQ=255)
                                                                            tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTa    >  2:1424112/1‑90 (MQ=255)
                                                                              aGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGc  <  2:71797/90‑1 (MQ=255)
                                                                                   |                                                                                    
AGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGC  >  NZ_CP009273/1394080‑1394247

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: