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,1630GA95.2% 68.5 / ‑3.9 21L111L (CTG→TTG) ogtmethylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (9/11);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC  >  NZ_CP009273/1394078‑1394239
                                                                                     |                                                                            
cAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTgg                                                                          <  1:805951/90‑1 (MQ=255)
  aGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGcc                                                                        >  2:553115/1‑90 (MQ=255)
            aTGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCAtt                                                              <  2:887283/90‑1 (MQ=255)
             tGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTa                                                             >  2:87315/1‑90 (MQ=255)
              gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                            >  1:445931/1‑90 (MQ=255)
              gCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTAc                                                            >  2:789902/1‑90 (MQ=255)
                 gATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTg                                                         <  1:7414/90‑1 (MQ=255)
                   tGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGcc                                                       <  2:14907/90‑1 (MQ=255)
                      gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGc                                                    >  1:362621/1‑90 (MQ=255)
                      gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGc                                                    <  1:789902/90‑1 (MQ=255)
                       aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCa                                                   >  2:569312/1‑90 (MQ=255)
                                        cACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAgtg                                  <  2:849526/90‑1 (MQ=255)
                                        cACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAgtg                                  <  2:962021/90‑1 (MQ=255)
                                                         gcgccAGGACGGCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCgcgc                 <  2:506169/90‑1 (MQ=255)
                                                                    gCCCGATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGGTAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt      >  1:997834/1‑90 (MQ=255)
                                                                    gCCCAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt      >  2:140100/1‑90 (MQ=255)
                                                                      ccAATTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc    <  1:974638/90‑1 (MQ=255)
                                                                        aaTTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcccc  >  2:140880/1‑90 (MQ=255)
                                                                        aaTTGCTCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcccc  <  2:22081/90‑1 (MQ=255)
                                                                              tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCgcgc                 >  2:464084/1‑69 (MQ=255)
                                                                              tCAGCCAATTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCgcgc                 <  1:464084/69‑1 (MQ=255)
                                                                                     |                                                                            
CAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC  >  NZ_CP009273/1394078‑1394239

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