Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,399,978 T→C T87A (ACG→GCG)  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,9780TC98.1% 173.9 / ‑3.1 52coding (259/516 nt)ogtmethylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/30);  minor base . (1/0);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

GGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAG  >  NC_000913/1399841‑1400113
                                                                                                                                          |                                                                                                                                       
gggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGccc                                                                                                                                       >  2:820517/1‑140 (MQ=255)
 ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGcccc                                                                                                                                      <  2:917705/140‑1 (MQ=255)
     tCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGcccccccc                                                                                                                                  >  1:1301516/1‑140 (MQ=255)
      cGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCg                                                                                                                                 <  1:820517/140‑1 (MQ=255)
        aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTa                                                                                                                               >  2:1180765/1‑140 (MQ=255)
        aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTa                                                                                                                               >  2:1153087/1‑140 (MQ=255)
               tGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTgg                                                                                                                        >  1:806678/1‑140 (MQ=255)
                 ccGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                      <  2:3613/140‑1 (MQ=255)
                 ccGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                      <  2:729611/140‑1 (MQ=255)
                 ccGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                      <  1:430526/140‑1 (MQ=255)
                 ccGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                      <  2:1111769/140‑1 (MQ=255)
                   gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAg                                                                                                                    <  2:770829/140‑1 (MQ=255)
                   gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAg                                                                                                                    <  1:325471/140‑1 (MQ=255)
                   gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAg                                                                                                                    <  2:196338/140‑1 (MQ=255)
                     aCCAACGGCACGCGCCGCGCCAGGACGCCCCAATTGCTCAGCCAGCTGGCCGGAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCTTAGCACAGGGAAGCg                                                                                                                  >  2:435115/1‑140 (MQ=255)
                      ccAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGt                                                                                                                 <  2:1338669/140‑1 (MQ=255)
                      ccAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGt                                                                                                                 <  2:254987/140‑1 (MQ=255)
                      ccAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGt                                                                                                                 <  2:560188/140‑1 (MQ=255)
                      ccAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGt                                                                                                                 <  2:923710/140‑1 (MQ=255)
                      ccAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGt                                                                                                                 <  1:990484/140‑1 (MQ=255)
                        aaCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTAt                                                                                                               <  1:774313/140‑1 (MQ=255)
                         aCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAg                                                                                                                    >  1:837360/1‑134 (MQ=255)
                         aCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAg                                                                                                                    <  2:837360/134‑1 (MQ=255)
                           ggCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                            >  1:641399/1‑140 (MQ=255)
                              aCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                      <  2:36709/127‑1 (MQ=255)
                              aCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                      >  1:36709/1‑127 (MQ=255)
                               cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                      <  2:412724/126‑1 (MQ=255)
                               cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                      >  1:412724/1‑126 (MQ=255)
                               cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                            >  1:596775/1‑136 (MQ=255)
                               cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                            <  2:596775/136‑1 (MQ=255)
                                 cgccgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGc                                                                                                      >  1:1274253/1‑140 (MQ=255)
                                                gCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGAtt                                                                                               >  2:1015060/1‑132 (MQ=255)
                                                gCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGAtt                                                                                               <  1:1015060/132‑1 (MQ=255)
                                                  ccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCaaaa                                                                                     >  2:988482/1‑140 (MQ=255)
                                                                               aTTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGCGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGcc                                                        <  2:1301516/140‑1 (MQ=255)
                                                                                   ccTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                                    <  1:988482/140‑1 (MQ=255)
                                                                                    cTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGAt                                                   >  2:768844/1‑140 (MQ=255)
                                                                                     tgcgggCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGAtt                                                  >  1:165180/6‑140 (MQ=255)
                                                                                      gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtg                                                 >  1:332974/1‑140 (MQ=255)
                                                                                      gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtg                                                 >  1:1232508/1‑140 (MQ=255)
                                                                                       cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtgg                                                >  2:601803/1‑140 (MQ=255)
                                                                                       cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtgg                                                >  2:1134827/1‑140 (MQ=255)
                                                                                                 aTAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGaaa                                       <  1:558992/140‑1 (MQ=255)
                                                                                                     tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                                  <  2:45566/140‑1 (MQ=255)
                                                                                                        gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCAt                              >  2:141547/1‑140 (MQ=255)
                                                                                                            tgtTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACTGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGc                           >  2:684718/1‑140 (MQ=255)
                                                                                                             gtTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGcc                          <  2:1201250/140‑1 (MQ=255)
                                                                                                               tttCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCctt                        >  1:859829/1‑140 (MQ=255)
                                                                                                                   cagcttccgcgCTAAAAGGCCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg                   <  1:747563/133‑1 (MQ=255)
                                                                                                                          cgcgCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATaa             >  2:874514/1‑140 (MQ=255)
                                                                                                                          cgcgCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATaa             >  1:1054598/1‑140 (MQ=255)
                                                                                                                              cTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGa         <  2:1054598/140‑1 (MQ=255)
                                                                                                                                   aTGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc    <  2:641399/140‑1 (MQ=255)
                                                                                                                                   aTGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc    <  1:601803/140‑1 (MQ=255)
                                                                                                                                   aTGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGACAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc    <  1:1180765/140‑1 (MQ=255)
                                                                                                                                    tGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCca   <  1:44144/140‑1 (MQ=255)
                                                                                                                                     ggCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCcag  <  2:332974/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                       
GGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAG  >  NC_000913/1399841‑1400113

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: