Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,399,985 (C)7→9 coding (252/516 nt) ogt ← methylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,399,9781.C100.0% 60.5 / NA 19T87A (ACG→GCG) ogtmethylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (8/11);  total (8/11)
*NC_0009131,399,9782.C94.7% 57.8 / ‑2.9 19coding (259/516 nt)ogtmethylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (8/10);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCG  >  NC_000913/1399846‑1400080
                                                                                                                                     ||                                                                                                       
tCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt                                                                                                           >  1:389770/1‑133 (MQ=255)
   aTCCGTTTGTGGCACGAACGGCAGTCCCCGCGTCAGGACGGGCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTGGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑‑ccc                                                                                                     <  2:250420/133‑1 (MQ=255)
   aTCCGTTTCCCGCACCAACGGCCCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑‑ccc                                                                                                     <  1:207458/133‑1 (MQ=255)
       gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑‑ccccccc                                                                                                 >  2:462646/1‑133 (MQ=255)
          tGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑‑cccccccccg                                                                                              <  2:328264/133‑4 (MQ=255)
           gCCGCACCAACGGCACGGCCCGCGCCAGGACGGCCCAATTGCTCAGCCAGGTGGGCGTAATGGTTTATCTTTCCGCGGGGGTTACCCCATAGTGGTGTCCAGTCTTCGCCCTGAACAGGTGc‑‑‑cccccccctgt                                                                                             >  2:127462/1‑129 (MQ=255)
             cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑‑CCCCCCCGTAGc                                                                                            <  1:497702/133‑1 (MQ=255)
                    aCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTggg                                                                                      <  1:88634/133‑1 (MQ=255)
                      ggCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCCGTAGCAGTGGGa                                                                                     <  1:358649/133‑1 (MQ=255)
                          cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATCCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGt                                                                                <  2:381240/133‑1 (MQ=255)
                           acgccgcaccACGACCGCGCAATTGATCTACCAGTTGCCGGTAACCCATTACCTGCCTGCTGGAGATAGTGCGTAGTGTTTTACAGTCTTCGCGCTGAACTGGCATC‑‑CCCCCCCGTAGCAGTGGAAAGCGTAt                                                                              <  2:475902/132‑1 (MQ=255)
                               cgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                           <  1:403005/133‑1 (MQ=255)
                                      ggACGGCCCAATTGCGCAGCCAGTGGGCTGTATTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCt                                                                    <  2:192513/133‑1 (MQ=255)
                                           gCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGAt                                                               <  1:127462/133‑1 (MQ=255)
                                            cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGAtt                                                              >  1:267627/1‑133 (MQ=255)
                                                           aGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTc                                               <  1:462646/133‑1 (MQ=255)
                                                                                 gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCCGTAGCAGGGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAAcc                          >  1:23957/1‑133 (MQ=255)
                                                                                 gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCg                         >  2:495216/1‑133 (MQ=255)
                                                                                   ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGTAAAATATTCACGAAGCTTGTCGCTTAAACCGcc                       <  2:408324/133‑1 (MQ=255)
                                                                                      cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTgg                    >  2:276830/1‑133 (MQ=255)
                                                                                                    tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGaaa      >  2:31230/1‑133 (MQ=255)
                                                                                                     aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAAt     >  1:439641/1‑133 (MQ=255)
                                                                                                      gtgtTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg    <  2:267627/133‑1 (MQ=255)
                                                                                                       tgtTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgc   >  2:18340/1‑133 (MQ=255)
                                                                                                                 aCTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcg  >  2:357084/1‑124 (MQ=255)
                                                                                                                 aCTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcg  <  1:357084/124‑1 (MQ=255)
                                                                                                                                     ||                                                                                                       
TCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCG  >  NC_000913/1399846‑1400080

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