Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,532,534 G→C 10.0% A217G (GCC→GGC)  yjhP ← KpLE2 phage‑like element; putative methyltransferase YjhP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,5340GC10.0% 133.3 / 3.1 50A217G (GCC→GGC) yjhPKpLE2 phage‑like element; putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base G (26/18);  new base C (3/3);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCC  >  NC_000913/4532396‑4532670
                                                                                                                                          |                                                                                                                                        
caGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCACAGCATTCCCGCGCGCAGGTCACGCAGCGTCCCGGCGCGACGCCTCGTTCCGCCCTCCCCCCggcggc                                                                                                                                      >  2:212580/2‑143 (MQ=255)
   ccTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTACGGCCCTGACTTCGGCGGCGaa                                                                                                                                   >  2:792674/1‑143 (MQ=255)
     tCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGCTGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGt                                                                                                                                 >  1:924402/1‑143 (MQ=255)
     tCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCAAAAGCATTCCCGCGCGTAGGTCACGTAGCGCTTCGGCGCGATGCTTAGTCCGGCCCTGACTTCGGAGGCGCAGt                                                                                                                                 >  2:284017/1‑143 (MQ=255)
        atatTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCAt                                                                                                                              >  1:154301/1‑143 (MQ=255)
        atatTTACAGCCAGTATGACCGGGAACCTGTTGTTACCGGGCGATGAGTGCTAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGGTTCGGCGCGATGTATAGTTCGGCCCTGACTTCGGCGGCGAAGTCAt                                                                                                                              <  2:146455/143‑1 (MQ=255)
          atggccagccagAAATACCGGGAACCTGTTGTTGCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCaggg                                                                                                                            <  1:94366/139‑4 (MQ=255)
             tCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGCTGTTTAGTTCGGCCCTGACTTCGGCGGCGAGGTccgcgacc                                                                                                                         >  2:318600/1‑136 (MQ=255)
                 gccagAAATAGCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTCGGTCACGTAGCGTTTCGGCGCGAGGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGAt                                                                                                                     >  2:162946/1‑143 (MQ=255)
                  ccagAAATACCGGGCACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTCGTTCGGCCCTGACTCCCGCGGCGACGTCGTCGTCAGGAtt                                                                                                                    >  2:910350/1‑143 (MQ=255)
                  ccagAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGAtt                                                                                                                    <  1:847657/143‑1 (MQ=255)
                            cGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTGCGGCCGCGAAGTCATCGTCAGGATTCGCTGccagc                                                                                                          >  2:912151/1‑143 (MQ=255)
                               gAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagccag                                                                                                       <  1:229212/143‑1 (MQ=255)
                                            cggggCGATGGGTGCGGGCACGCCCCAGCCAGAGCGTTCCGGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGGAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCa                                                                                          <  2:630125/141‑1 (MQ=255)
                                              gggCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGc                                                                                        <  1:147237/143‑1 (MQ=255)
                                                gCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCa                                                                                      >  1:509364/1‑143 (MQ=255)
                                                      cgcgCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                                                                <  1:166568/140‑1 (MQ=255)
                                                        tGCGAGCACGCCCCAGCCAAAGCAGTCCGGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCgg                                                                              <  2:688358/143‑1 (MQ=255)
                                                        tGCGACCACGCCCCAGCCAACGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGGCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCgg                                                                              <  2:439440/143‑1 (MQ=255)
                                                        tGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGCTGTTTAGTTCGGCCCTGCCTTCGCCGGCGAAGTAATCGTCAGGATTCGCTTCCAGACCGCGGCGCATGGTCAGCCATCTTGCgg                                                                              >  1:426097/1‑143 (MQ=255)
                                                           gAACACGCCCCAGCCAAAGCATTCCCGGGCGGAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGAGAAGCCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                                                           <  1:794209/143‑1 (MQ=255)
                                                           gAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                                                           <  1:556437/143‑1 (MQ=255)
                                                                 gCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGATTCGGCGCGATGTTTGGTTCGGCCCTGACTTCGGCGGCGAGGTCAGCGTCAGGATTCGCATGCGGCCAGAGGCGGGGAGTCAGCCGTTGTGCGGCTGCATAcc                                                                     >  2:59042/1‑143 (MQ=255)
                                                                     cAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCATCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTc                                                                 >  1:684921/1‑143 (MQ=255)
                                                                        ccAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAGGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCCGCCATTTTGCGGCTTCATACCGGTCCCa                                                              >  1:326267/1‑143 (MQ=255)
                                                                          aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCCTCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCCGc                                                            >  1:793383/1‑143 (MQ=255)
                                                                          aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGTCTGCATACCGGTCCCGGc                                                            >  1:547657/1‑143 (MQ=255)
                                                                          aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGc                                                            >  1:509978/1‑143 (MQ=255)
                                                                          aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCAGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGCAGTCATCGTCAGGATGCGCTTCCAGCCAGCGGCGCCTGATCAGCCATTGTGCGGCTTAATACCGGTCCCaac                                                            >  2:648409/1‑141 (MQ=255)
                                                                           aaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCGCAGcc                                                           >  1:905581/1‑143 (MQ=255)
                                                                            aGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCt                                                          >  1:913991/1‑143 (MQ=255)
                                                                                  cccGCGCGTAGTCCACGTAGCGTTTCGGCGCGATGGTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTgg                                                    <  2:542920/143‑1 (MQ=255)
                                                                                     gcgcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt                                                 <  1:533598/143‑1 (MQ=255)
                                                                                       gcgTAGGTTAGGTGTCGTTTCGGCGCGAGGCTTAGTTCGCCCCTGTCTGCGGCGGCGAAGTCATCGTCAGGATTCCCTTCCAGCCCGCGGCGCATGGTCAGCCATTTCGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGc                                               <  2:735433/143‑1 (MQ=255)
                                                                                                aCGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGAGTCGCTTCCAGCCAGCGGCGCATGGTCAGCCGTTTTGCGGCTTCATACCGGTCCGCGCCTTCCTGGTCCTCCAGTACCaa                                      >  2:334557/1‑142 (MQ=255)
                                                                                                    aGCGTTGCGGCGCCATGTGTCGTTCGGCCCTGACCTCGCCGGCGCAGCCATCGTCAGGATTCGCTTCCCGCCAGCGGCGCCTGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGCACCATTTcc                                  <  2:744784/143‑1 (MQ=255)
                                                                                                     gCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAGTTTGCGGCTTCATACCGGTCCCCGCATTCCTGGTCTGCCAGTACCATTTcca                                 >  2:751800/1‑143 (MQ=255)
                                                                                                       gTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTccacc                               >  1:715354/1‑143 (MQ=255)
                                                                                                                cgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCCCGGCATAg                      >  1:662080/1‑143 (MQ=255)
                                                                                                                    gTTTAGCTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCga                  <  1:654038/143‑1 (MQ=255)
                                                                                                                     tttAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCgag                 <  2:333308/143‑1 (MQ=255)
                                                                                                                       tAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAt               <  2:154301/143‑1 (MQ=255)
                                                                                                                           tCGGCCCTGCCTCCGCGGGGCAATTGATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAtcgtc           <  2:700692/143‑1 (MQ=255)
                                                                                                                           tCGGCCCTGACTTCGGCGGCGAAGTCAACGTCGGCGTGCTCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAtcgtc           <  1:637092/143‑1 (MQ=255)
                                                                                                                           tCGGCCCGGACTTCCGCGGCGAAGTCAGCGTCCGGATTCGCTTCCAGCCAGCGGCGCATGGTCCGTCCTCTTGCTTCTTCCTACCCGGCCCAGCCGTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAtcgtc           <  1:793699/143‑1 (MQ=255)
                                                                                                                            tgacgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCCGCCAGTACCATTTCCGCCACGGCATAGCCGCGCtcggcg          >  1:370516/7‑143 (MQ=255)
                                                                                                                             ggCCCTGACTTCGGCGGCGAAGTCATCGTCCGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCCTTTTGCGGCTTCATACCGGTCCCAGCCGTACTGGTCTGCCGGTCCCATTTCCCGCACGTCCTCGCCGCGAtctgcct         >  1:665643/1‑138 (MQ=255)
                                                                                                                               cccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAtcgtcggag       >  1:677655/1‑140 (MQ=255)
                                                                                                                               cccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCGTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCGACCACGTCATAGCGGAGATCGTCGAAc       >  1:83238/1‑143 (MQ=255)
                                                                                                                                ccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAGCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCCCGTCATAGCCGAGATCGTCGAACg      >  2:759699/1‑143 (MQ=255)
                                                                                                                                    aCTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCGTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATACACGAGGTAGTCGAAccatcc  >  2:322153/1‑138 (MQ=255)
                                                                                                                                          |                                                                                                                                        
AAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCC  >  NC_000913/4532396‑4532670

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