Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,399,985 (C)7→8 coding (252/516 nt) ogt ← O‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,399,9781.C100.0% 118.9 / NA 35T87A (ACG→GCG) ogtO‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (17/18);  total (17/18)

GGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCC  >  NC_000913/1399842‑1400111
                                                                                                                                         |                                                                                                                                     
ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc                                                                                                                                     <  2:1322197/139‑1 (MQ=255)
           gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                           <  2:422878/139‑1 (MQ=255)
                 cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                                     <  1:874025/139‑1 (MQ=255)
                    aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                                  <  1:520424/139‑1 (MQ=255)
                            aatgcgccgcgccACGACGGCGCAATTGCTCAGCTAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAAAGGAGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                          <  1:1083981/136‑1 (MQ=255)
                              cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                                     <  1:1113625/126‑1 (MQ=255)
                              cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                                     >  2:1113625/1‑126 (MQ=255)
                                  ccgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGc                                                                                                    >  2:936276/1‑139 (MQ=255)
                                                     ttGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc                                                                                                                                     >  2:523563/1‑86 (MQ=255)
                                                     ttGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc                                                                                                                                     <  1:523563/86‑1 (MQ=255)
                                                                     ccGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCg                                                                 >  1:623715/1‑139 (MQ=255)
                                                                       gTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTTTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                                     <  2:869903/85‑1 (MQ=255)
                                                                       gTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                                     >  1:869903/1‑85 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGAt                                                 >  2:1149732/1‑139 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGTCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                          >  2:918382/1‑80 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                                    >  2:1066953/1‑70 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                                    <  1:1066953/70‑1 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                          <  1:918382/80‑1 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACg                                                                           >  1:214229/1‑111 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACg                                                                           <  2:214229/111‑1 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtg                                               >  2:549882/1‑139 (MQ=255)
                                                                                                    tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaataa                                                                                                       <  2:348371/70‑1 (MQ=255)
                                                                                                    tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc                                                                                                                                     <  1:1202923/39‑1 (MQ=255)
                                                                                                    tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc                                                                                                                                     >  2:1202923/1‑39 (MQ=255)
                                                                                                    tGCGTAGTGTTTTCCAGACCTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaataa                                                                                                       >  1:348371/1‑70 (MQ=255)
                                                                                                    ggCGTTGTGTTTTCCAGACTTCGGGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgc                                  <  1:553737/138‑1 (MQ=255)
                                                                                                       gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAAcc                                                         >  1:39581/1‑113 (MQ=255)
                                                                                                       gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAAcc                                                         <  2:39581/113‑1 (MQ=255)
                                                                                                          gtgtTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCAt                            <  1:549882/139‑1 (MQ=255)
                                                                                                              tttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                           >  1:464165/1‑40 (MQ=255)
                                                                                                              tttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                           <  2:464165/40‑1 (MQ=255)
                                                                                                              tttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGcc                        >  2:759478/1‑139 (MQ=255)
                                                                                                                     aCTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                                     <  2:1045807/87‑1 (MQ=255)
                                                                                                                     aCTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                                     >  1:1045807/1‑87 (MQ=255)
                                                                                                                         cgcgCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                                      >  2:810546/1‑66 (MQ=255)
                                                                                                                         cgcgCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                                      <  1:810546/66‑1 (MQ=255)
                                                                                                                         cgcgCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                                  >  1:495770/1‑102 (MQ=255)
                                                                                                                         cgcgCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                                  <  2:495770/102‑1 (MQ=255)
                                                                                                                          gcgcTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATa            >  2:189819/1‑139 (MQ=255)
                                                                                                                                    ggCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc  <  1:1149732/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
GGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCC  >  NC_000913/1399842‑1400111

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