Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,032,336 G→C 11.5% intergenic (‑19/+48) dcm ← / ← yedJ DNA‑cytosine methyltransferase/putative HD superfamily phosphohydrolase YedJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,032,3360GC11.5% 92.0 / 4.9 35intergenic (‑19/+48)dcm/yedJDNA‑cytosine methyltransferase/putative HD superfamily phosphohydrolase YedJ
Reads supporting (aligned to +/- strand):  ref base G (18/13);  new base C (2/2);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

TACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCG  >  NC_000913/2032198‑2032474
                                                                                                                                          |                                                                                                                                          
tACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCGGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGAt                                                                                                                                        <  1:433523/143‑1 (MQ=255)
     cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCg                                                                                                                                   <  1:27920/143‑1 (MQ=255)
     cATTAAGCTGCGCCACCAACGTTTTAACATCAGAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCg                                                                                                                                   <  2:512507/143‑1 (MQ=255)
       ttAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGACCGAc                                                                                                                                 >  2:622595/1‑143 (MQ=255)
           gCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                             >  1:437873/1‑143 (MQ=255)
                 ccaccaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTCCATCGATTTCACCGGCCATCTGATATGGCAGACAGGTTAACaa                                                                                                                       >  1:461224/1‑142 (MQ=255)
                       aCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATATATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTAc                                                                                                                 <  2:46420/143‑1 (MQ=255)
                        cgctttGACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                <  2:245858/140‑1 (MQ=255)
                        cGTTTTAACATCATACATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATGTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                <  2:506361/143‑1 (MQ=255)
                                catcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTa                                                                                                        >  1:656559/1‑143 (MQ=255)
                                  tcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACGGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAgg                                                                                                      <  2:631237/143‑1 (MQ=255)
                                      aaaTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCCTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATCATTAGAGACTTTTAaggtga                                                                                                  >  1:405646/1‑140 (MQ=255)
                                       aaTTTGCATCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAc                                                                                                 <  1:117603/143‑1 (MQ=255)
                                                     tctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                                                                   <  2:651323/143‑1 (MQ=255)
                                                       tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGa                                                                                 >  1:438536/1‑143 (MQ=255)
                                                                 ccTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCa                                                                       <  1:185423/143‑1 (MQ=255)
                                                                   tgtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAgg                                                                                                                              <  1:267209/86‑1 (MQ=255)
                                                                   tgtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAgg                                                                                                                              >  2:267209/1‑86 (MQ=255)
                                                                      gCGGCACTCCCGGTGCTGGATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTAt                                                                  <  1:75696/143‑1 (MQ=255)
                                                                               ccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTCAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGCAAACGCGTCTATTACCTtgtg                                                         >  2:298690/1‑143 (MQ=255)
                                                                               ccGGTGCTGTATGAATCGGTCACTGTTCTATTTTCCTGCATAGATTTCCCCGGCCACCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgtg                                                         <  2:411747/143‑1 (MQ=255)
                                                                                 ggTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCCCGCGCTCCGCCCGCGGATGAAAGCGCATCTATCAACTTGTGaa                                                       >  1:116695/1‑142 (MQ=255)
                                                                                     cTGTGTGAATCGGGTTGTGATATAGTTCCCTGGATTGATTTCACCGGCCTTCTGAGGTGGCCGACGGGATAGCGTAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac                                                   <  2:167141/143‑1 (MQ=255)
                                                                                                      tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCg                                  >  2:411739/1‑143 (MQ=255)
                                                                                                      tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCg                                  >  2:179886/1‑143 (MQ=255)
                                                                                                                    gCATAGATTTCCCCGTCGATCTGGGATGGGCGACAGTTTAACATAATTACAGCTATTTAGGGCCACGCGCTTAGTCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACATCTTCATTCTCCCCCGCCAGTTCGGCACTg                    <  2:67995/143‑1 (MQ=255)
                                                                                                                      aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAg                  >  2:45242/1‑143 (MQ=255)
                                                                                                                         gATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCtt               >  2:102495/1‑143 (MQ=255)
                                                                                                                           tttCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCACTCTCCCCCGCCAGCTCGGCACTGAGCTCCg             >  2:308099/1‑143 (MQ=255)
                                                                                                                             tCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGGCACCTTCATCCTCCCCCGCCCGTTAGACCCTGAGCTTCGcc           >  2:112858/1‑143 (MQ=255)
                                                                                                                                 cGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCAtga       >  1:145064/1‑141 (MQ=255)
                                                                                                                                  ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCCGCACTGAGCTTCGCCATaaa      >  2:287217/1‑143 (MQ=255)
                                                                                                                                   gCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATacac     >  1:247666/1‑143 (MQ=255)
                                                                                                                                    ccATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGACACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCGCGCCCAGTTCTGCACTTAGCATCGCCACAAAcg    >  2:653163/1‑142 (MQ=255)
                                                                                                                                      aTCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCCCTGAGCTTCACCATAAACTCg  >  1:580911/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCG  >  NC_000913/2032198‑2032474

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