Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,030,947 C→T 6.8% A457A (GCG→GCA)  dcm ← DNA cytosine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,030,9470CT 6.8% 182.3 / 2.7 59A457A (GCG→GCA) dcmDNA cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (23/32);  new base T (2/2);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGC  >  NC_000913/2030807‑2031092
                                                                                                                                            |                                                                                                                                                 
ccAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGtt                                                                                                                                            <  8:356767/148‑1 (MQ=255)
 cAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGttt                                                                                                                                           >  7:435492/1‑148 (MQ=255)
 cAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGttt                                                                                                                                           >  8:121307/1‑148 (MQ=255)
  aGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGttt                                                                                                                                           >  1:163089/1‑147 (MQ=255)
     cccTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGAtt                                                                                                                                       <  5:274610/148‑1 (MQ=255)
     cccTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGAtt                                                                                                                                       <  5:268404/148‑1 (MQ=255)
        tGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGAttttt                                                                                                                                    >  1:468312/1‑148 (MQ=255)
         gCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTg                                                                                                                                   >  1:310449/1‑148 (MQ=255)
         gCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGAATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTg                                                                                                                                   >  6:408575/1‑148 (MQ=255)
           ccGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGc                                                                                                                                 <  2:232958/148‑1 (MQ=255)
            cGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCTCGCGTGGCAATAGCGCGCATATTTTTGCTGCGTGTGGCCTTATCGTTAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGGCACCGCCTGTTTGATTTTTGGCt                                                                                                                                <  2:418367/148‑1 (MQ=255)
              gTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCa                                                                                                                              <  6:315840/148‑1 (MQ=255)
                    aGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGt                                                                                                                        >  7:353490/1‑148 (MQ=255)
                     gACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGGGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGtt                                                                                                                       <  8:317281/148‑1 (MQ=255)
                     gACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTACTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTg                                                                                                                                          <  2:236074/129‑1 (MQ=255)
                     gACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTACTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTg                                                                                                                                          >  1:236074/1‑129 (MQ=255)
                      aCTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGttt                                                                                                                      <  8:353490/148‑1 (MQ=255)
                      aCTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGGGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGttt                                                                                                                      <  1:334433/148‑1 (MQ=255)
                          gCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGcc                                                                                                                  <  2:392628/148‑1 (MQ=255)
                                cgcTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACcgcg                                                                                                            <  8:435492/148‑1 (MQ=255)
                                        tATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTgg                                                                                                                                  >  3:242200/1‑118 (MQ=255)
                                        tATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTgg                                                                                                                                  <  4:242200/118‑1 (MQ=255)
                                        tATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAg                                                                                                      <  2:468312/146‑1 (MQ=255)
                                        tATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGAc                                                                                                    >  3:142343/1‑148 (MQ=255)
                                            gCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCgg                                                                                                 <  6:59390/147‑1 (MQ=255)
                                            gCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGc                                                                                                <  2:256524/148‑1 (MQ=255)
                                                     cgcgTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGaccacc                                                                                       <  2:211132/148‑1 (MQ=255)
                                                      gcgTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGACCGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCGCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCg                                                                                      <  1:293417/148‑1 (MQ=255)
                                                        gTGGCAATCGCGCGAATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTTCTGACGCAACTCCACCGCCTGTTTGATATTTGGCTCAAGCAGTTTTGCCAACGCGTCAAAGACCGGCACGACCACCGAg                                                                                    <  5:439969/148‑1 (MQ=255)
                                                         tGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGt                                                                                   <  2:310449/148‑1 (MQ=255)
                                                           gCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTa                                                                                 <  3:370403/148‑1 (MQ=255)
                                                           gCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTa                                                                                 <  7:121307/148‑1 (MQ=255)
                                                                cgcgcgCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGAc                                                                                                    <  8:379943/124‑1 (MQ=255)
                                                                cgcgcgCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGAc                                                                                                    >  7:379943/1‑124 (MQ=255)
                                                                 gcgcgcATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGaa                                                                            <  2:163089/147‑1 (MQ=255)
                                                                      cATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGc                                                                       >  3:189311/1‑147 (MQ=255)
                                                                        tatTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGAt                                                                    <  5:408575/148‑1 (MQ=255)
                                                                         atTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGa                                                                     >  7:23053/1‑146 (MQ=255)
                                                                            tttGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAgg                                                                 <  3:378829/147‑1 (MQ=255)
                                                                             ttGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGc                                                                <  8:23053/147‑1 (MQ=255)
                                                                                     aggggCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTc                                                       <  6:17584/145‑1 (MQ=255)
                                                                                           cTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACTGCCTGTTTGAttttt                                                                                                                                    <  5:48895/65‑1 (MQ=255)
                                                                                           cTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACTGCCTGTTTGAttttt                                                                                                                                    >  6:48895/1‑65 (MQ=255)
                                                                                            ttATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTg                                                             >  1:177923/1‑135 (MQ=255)
                                                                                            ttATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTg                                                             <  2:177923/135‑1 (MQ=255)
                                                                                             tATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCg                                               >  1:421071/1‑148 (MQ=255)
                                                                                             tATCGTGAACGTCGGCCATGTTGTACCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGa                                                                             >  4:172425/1‑118 (MQ=255)
                                                                                             tATCGTGAACGTCGGCCATGTTGTACCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGa                                                                             <  3:172425/118‑1 (MQ=255)
                                                                                                  agAACGTCGGCCAAGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTATTCGCTCAAGCAGTTTTGACACCGCGGCAAAGACCGGCACGACCACCGAGTTACAGAACTGGCGATAGGCCTGACTGTCCGAAACCGGAATa                                          >  6:297530/2‑148 (MQ=255)
                                                                                                                gttgtGCCTCTTGCTGACGCAACGCCACTGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAAc                                                                           >  7:431110/1‑101 (MQ=255)
                                                                                                                gttgtGCCTCTTGCTGACGCAACGCCACTGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAAc                                                                           <  8:431110/101‑1 (MQ=255)
                                                                                                                gttgtGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTtc                            >  7:331599/1‑148 (MQ=255)
                                                                                                                         ctTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGgcgc                   >  3:237698/1‑148 (MQ=255)
                                                                                                                         ctTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCACGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGgcgc                   >  5:422595/1‑148 (MQ=255)
                                                                                                                                 cGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGTGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCaaaa            <  4:400600/147‑1 (MQ=255)
                                                                                                                                 cGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGTGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCaaaa            >  3:400600/1‑147 (MQ=255)
                                                                                                                                     aCGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATa                                                                   >  7:339507/1‑88 (MQ=255)
                                                                                                                                     aCGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATa                                                                   <  8:339507/88‑1 (MQ=255)
                                                                                                                                          aCCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgc  >  1:220549/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCGCGCGCATATTTTTGCTGCGAGTGGCCTTATCGTGAACGTCGGCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGC  >  NC_000913/2030807‑2031092

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