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.C97.7% 143.8 / ‑3.7 43coding (259/516 nt)ogtO‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (18/24);  total (18/25)
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.99e-01

GGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑CCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAT  >  NC_000913/1399842‑1400100
                                                                                                                                         |                                                                                                                           
ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           <  1:533569/139‑1 (MQ=255)
ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           <  1:79971/139‑1 (MQ=255)
ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           <  1:581052/139‑1 (MQ=255)
     cGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑cccc                                                                                                                     <  1:1088323/139‑1 (MQ=255)
       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑cccccg                                                                                                                   >  1:721209/1‑137 (MQ=255)
       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑cccccg                                                                                                                   >  2:1147828/1‑137 (MQ=255)
         ccGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑cccccgta                                                                                                                 <  2:653852/139‑5 (MQ=255)
          cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑cccccgtag                                                                                                                <  1:411357/139‑6 (MQ=255)
           gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGc                                                                                                                <  2:721209/139‑1 (MQ=255)
           gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTAAAATGGCGTCCCC‑CCCCGTAGc                                                                                                                <  1:271687/139‑1 (MQ=255)
                  gCACCAACGGCACACGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGa                                                                                                         <  2:872533/139‑1 (MQ=255)
                   cACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGaa                                                                                                        >  1:727811/1‑139 (MQ=255)
                    aCCAACGGCGCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑cccc                                                                                                                     <  2:729674/124‑1 (MQ=255)
                    aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑cccc                                                                                                                     >  1:729674/1‑124 (MQ=255)
                    aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           >  2:1630446/1‑119 (MQ=255)
                    aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           <  1:1630446/119‑1 (MQ=255)
                       aaCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGt                                                                                                    <  1:1147828/139‑1 (MQ=255)
                       aaCAGCACGCGCCGCGCCAGGACGGTCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGt                                                                                                    <  1:1069006/139‑1 (MQ=255)
                            cACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCaa                                                                                               <  1:96600/139‑1 (MQ=255)
                             aCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCaat                                                                                              >  2:1480366/1‑139 (MQ=255)
                                                cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTAt                                                                                                  <  1:781157/116‑1 (MQ=255)
                                                cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTAt                                                                                                  >  2:781157/1‑116 (MQ=255)
                                                                gTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCt                                                           >  2:995128/1‑139 (MQ=255)
                                                                 ttGGCCGTAATGCATTTCCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                          <  1:995128/139‑1 (MQ=255)
                                                                 ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                          <  2:727811/139‑1 (MQ=255)
                                                                 ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                          >  1:303899/1‑139 (MQ=255)
                                                                           tGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaaa                                                <  2:766570/139‑1 (MQ=255)
                                                                               ttACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGc                                            >  2:1318938/1‑139 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAg                                                                                                       >  1:420684/1‑74 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAg                                                                                                       <  2:420684/74‑1 (MQ=255)
                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑cccccg                                                                                                                   >  2:945940/1‑58 (MQ=255)
                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑cccccg                                                                                                                   <  1:945940/60‑3 (MQ=255)
                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑c                                                                                                                        <  1:896141/55‑1 (MQ=255)
                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc‑c                                                                                                                        >  2:896141/1‑55 (MQ=255)
                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           <  1:941407/53‑1 (MQ=255)
                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           <  1:1157262/53‑1 (MQ=255)
                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           >  2:1157262/1‑53 (MQ=255)
                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           >  2:941407/1‑53 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           <  1:590695/52‑1 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                           >  2:590695/1‑52 (MQ=255)
                                                                                              ggATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAg                             <  1:78958/139‑1 (MQ=255)
                                                                                                   gTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAATATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                       >  1:189688/1‑124 (MQ=255)
                                                                                                   gTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAATATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                       <  2:189688/124‑1 (MQ=255)
                                                                                                   gTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                        <  2:488538/139‑1 (MQ=255)
                                                                                                       gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAg                                                                                                       <  2:291254/56‑1 (MQ=255)
                                                                                                       gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAg                                                                                                       >  1:291254/1‑56 (MQ=255)
                                                                                                       gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCaata                                                                                             <  1:535646/66‑1 (MQ=255)
                                                                                                       gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCaata                                                                                             >  2:535646/1‑66 (MQ=255)
                                                                                                        tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGa                                                                                                         <  2:1415080/53‑1 (MQ=255)
                                                                                                        tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGa                                                                                                         <  1:1276398/53‑1 (MQ=255)
                                                                                                        tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGa                                                                                                         >  1:1415080/1‑53 (MQ=255)
                                                                                                        tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGa                                                                                                         >  2:1276398/1‑53 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCaat                                                                                              >  2:1100967/1‑63 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCaat                                                                                              <  1:1100967/63‑1 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                                 <  1:480380/108‑1 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGATAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                                 >  2:480380/1‑108 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑CCCACCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCa                  <  1:1001749/139‑1 (MQ=255)
                                                                                                          gtgtTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGa                                                                                                         >  2:829742/1‑51 (MQ=255)
                                                                                                          gtgtTTTCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGa                                                                                                         <  1:829742/51‑1 (MQ=255)
                                                                                                             ttttCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                     <  2:612429/132‑1 (MQ=255)
                                                                                                             ttttCCAGACTTCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                     >  1:612429/1‑132 (MQ=255)
                                                                                                                       ttCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTAc                                                                                <  1:1541570/63‑1 (MQ=255)
                                                                                                                       ttCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTAc                                                                                >  2:1541570/1‑63 (MQ=255)
                                                                                                                       ttCGCGCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCg    >  2:1218463/1‑139 (MQ=255)
                                                                                                                         cgcgCTGAAATGGCGTCCCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAt  <  1:1218463/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                           
GGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑CCC‑CCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAT  >  NC_000913/1399842‑1400100

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