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% 149.2 / ‑3.0 44T87A (ACG→GCG) ogtO‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (17/26);  minor base G (0/1);  total (17/27)
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.35e-01

TGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTT  >  NC_000913/1399840‑1400095
                                                                                                                                           |                                                                                                                     
tGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt                                                                                                                        <  1:463561/139‑1 (MQ=255)
 gggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c                                                                                                                      <  1:629804/139‑1 (MQ=255)
          tCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccgt                                                                                                             >  1:1811571/1‑136 (MQ=255)
            cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTGCCCCCCCGTAg                                                                                                            <  1:292185/139‑1 (MQ=255)
              tttGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCa                                                                                                          >  1:527800/1‑139 (MQ=255)
                tGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccg                                                                                                              <  2:348653/132‑3 (MQ=255)
                tGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccg                                                                                                              >  1:348653/1‑130 (MQ=255)
                   cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                     <  2:527800/139‑1 (MQ=255)
                      aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                  <  1:185511/139‑1 (MQ=255)
                           cGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTAt                                                                                             <  2:1190341/139‑1 (MQ=255)
                                cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                     <  2:985518/126‑1 (MQ=255)
                                cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                     >  1:985518/1‑126 (MQ=255)
                                    ccgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGt                                                                                               <  1:728754/128‑1 (MQ=255)
                                    ccgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGt                                                                                               >  2:728754/1‑128 (MQ=255)
                                     cgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCt                                                                                   <  2:1029901/139‑1 (MQ=255)
                                      gcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                     <  1:1955386/103‑1 (MQ=255)
                                      gcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                     >  2:1955386/1‑103 (MQ=255)
                                             gACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTAc                                                                           <  1:1573193/139‑1 (MQ=255)
                                                 gCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                     <  2:7130/109‑1 (MQ=255)
                                                 gCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                     >  1:7130/1‑109 (MQ=255)
                                                  cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c                                                                                                                      <  2:1670751/90‑1 (MQ=255)
                                                  cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑c                                                                                                                      >  1:1670751/1‑90 (MQ=255)
                                                    cAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                          >  2:1181315/1‑117 (MQ=255)
                                                    cAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGAATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGGAGTGGGAAGCGTATCaa                                                                                          <  1:1181315/117‑1 (MQ=255)
                                                       ttGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccg                                                                                                              <  2:1553927/93‑3 (MQ=255)
                                                       ttGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccg                                                                                                              >  1:1553927/1‑91 (MQ=255)
                                                            cAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCAc                                                            >  1:1090307/1‑139 (MQ=255)
                                                                  gTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCt                                                      >  1:1345559/1‑139 (MQ=255)
                                                                   ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                     <  2:1811571/139‑1 (MQ=255)
                                                                   ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                     <  1:1987097/139‑1 (MQ=255)
                                                                    tGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                  <  1:808748/93‑1 (MQ=255)
                                                                    tGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                  >  2:808748/1‑93 (MQ=255)
                                                                            aTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                            <  2:1345559/139‑1 (MQ=255)
                                                                             tGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                     >  2:1177576/1‑64 (MQ=255)
                                                                             tGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                     <  1:1177576/64‑1 (MQ=255)
                                                                              gCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAAc                                          <  2:1158763/139‑1 (MQ=255)
                                                                                      tGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                  <  2:1219197/139‑1 (MQ=255)
                                                                                       gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                     <  1:417245/54‑1 (MQ=255)
                                                                                       gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                     >  2:417245/1‑54 (MQ=255)
                                                                                        cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGTTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGAtt                                >  2:1228820/1‑139 (MQ=255)
                                                                                         ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc                                                                                                                    >  2:1957541/1‑53 (MQ=255)
                                                                                         ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccc                                                                                                                    <  1:1957541/53‑1 (MQ=255)
                                                                                            cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtggtg                            >  1:937364/1‑139 (MQ=255)
                                                                                            cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccgta                                                                                                            >  1:1325776/1‑54 (MQ=255)
                                                                                            cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccgta                                                                                                            <  2:1325776/58‑5 (MQ=255)
                                                                                                   tAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtgg                              <  2:1273335/130‑1 (MQ=255)
                                                                                                   tAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtgg                              >  1:1273335/1‑130 (MQ=255)
                                                                                                   tAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGaaa                     >  1:816129/1‑139 (MQ=255)
                                                                                                     gTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                   <  2:1631893/139‑1 (MQ=255)
                                                                                                     gTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                   <  1:1860399/139‑1 (MQ=255)
                                                                                                        cGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                <  2:1090307/139‑1 (MQ=255)
                                                                                                         gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                  >  1:983008/1‑56 (MQ=255)
                                                                                                         gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                  <  2:983008/56‑1 (MQ=255)
                                                                                                         gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                            >  2:118854/1‑110 (MQ=255)
                                                                                                         gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                            <  1:118854/110‑1 (MQ=255)
                                                                                                         gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCt               >  1:888816/1‑139 (MQ=255)
                                                                                                            gtgtTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCAt            <  1:1855160/139‑1 (MQ=255)
                                                                                                               ttttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                     >  2:411471/1‑95 (MQ=255)
                                                                                                               ttttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                     <  1:411471/95‑1 (MQ=255)
                                                                                                               ttttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGc         <  2:888816/139‑1 (MQ=255)
                                                                                                                      gACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCttt  <  2:1095335/139‑1 (MQ=255)
                                                                                                                       aCTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                  <  2:1007859/42‑1 (MQ=255)
                                                                                                                       aCTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                  >  1:1007859/1‑42 (MQ=255)
                                                                                                                                           |                                                                                                                     
TGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTT  >  NC_000913/1399840‑1400095

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: ‑