Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,532,534 G→C 15.6% 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,5340GC15.6% 52.7 / 4.3 24A217G (GCC→GGC) yjhPKpLE2 phage‑like element; putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base G (10/10);  new base C (2/2);  total (12/12)
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.66e-01

CGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAA  >  NC_000913/4532402‑4532675
                                                                                                                                    |                                                                                                                                             
cGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTCGGTCACGTAGCGTTTCGGCGCGCTGTTTAGTTCGGCCCTGACTTCGGCGACGACGTc                                                                                                                                     >  1:58791/1‑143 (MQ=255)
   tatTTCCAGCCAGAAATACCGGGAAGCTGGTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagg                                                                                                                                  <  2:446554/143‑3 (MQ=255)
   tagttCCACCCAGAAATACCGGGAACCTGTTGTTGCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATc                                                                                                                                  <  2:450803/140‑1 (MQ=255)
   gatTTCGCGCCAGAGATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTGCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATc                                                                                                                                  <  1:637699/142‑1 (MQ=255)
     tttCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCACAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAATTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                >  1:63599/1‑143 (MQ=255)
     tgtCCGGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                <  2:283526/141‑1 (MQ=255)
              agcgaTACCGGGAACCCGTTGTTACCGGGCGATGAGTGCGAACGCGCGCCAGCCAAAGCGTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCg                                                                                                                       <  2:204691/139‑1 (MQ=255)
                   gaCCGGGAACTTGTTGTTACTGGGCGATGATTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                                                                  <  2:27574/142‑1 (MQ=255)
                    aCCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcca                                                                                                                 <  1:67375/143‑1 (MQ=255)
                     ccGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccag                                                                                                                >  2:54710/1‑143 (MQ=255)
                                                               cAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGGGGCTTCATACCGGTc                                                                      >  2:21791/1‑143 (MQ=255)
                                                                    aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCCTTTTGCGGCTTCATACCGGGCCCCGc                                                                 >  2:152398/1‑143 (MQ=255)
                                                                    aaaGCATTCCCGCGCGTAGGTCACGTAGAGTTTCGGCGCTATTTTGAGTTCGGCCCTGACTGCGGCGGCGGAGTCATCGACTGGAGGCGCGTCCAGCCAACGGCGCAGGGTCAGCCATCTTGCGGCTTCATACCGGTCCCAGc                                                                 >  2:499595/1‑143 (MQ=255)
                                                                          gtCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                                           <  2:631229/142‑1 (MQ=255)
                                                                           gcccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                                                          <  1:133916/142‑1 (MQ=255)
                                                                             ccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTggc                                                        >  1:579381/1‑142 (MQ=255)
                                                                             ccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGt                                                        >  1:460611/1‑143 (MQ=255)
                                                                                                cGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTccac                                     >  1:242248/1‑143 (MQ=255)
                                                                                                         gcgATGTTTAGTTCGGCCCTGAATACGGCGGCGAAGTAAGCGTCAGGAATCGCTTCCAGCCCGCGGCGCATGGTCAGCAATTTTGAGGCGGCATACGGGGCCCAGCCCTCCGGGTCTTCGGCTACCTTCCCCACCACGTCAtc                            >  2:76862/1‑142 (MQ=255)
                                                                                                             tGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCGTTGTCGGTAGCCTCATCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCg                        <  1:649621/143‑1 (MQ=255)
                                                                                                                        gCCCTGACTTCGGATGCGATGTCTACGTCAGCATTCGCTTCCAGCCAGCGGCGCATGGTCATCCATTTTGCGGCTTCATACCGTTCCCAGCCTTCCTGGTCTGCCATTACCATTTCCACCACGTCATAGCCGAGATCGTCGaa             <  2:103007/143‑1 (MQ=255)
                                                                                                                          gatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACg           >  2:518852/3‑143 (MQ=255)
                                                                                                                          gatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGACCGCCGAACg           >  1:534358/3‑143 (MQ=255)
                                                                                                                                 tCGCCGGCGAAGTCATCGTCTGGATTCGCCGCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGAc    <  1:404705/143‑1 (MQ=255)
                                                                                                                                   ggcggcGAAGGCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATCCCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACaa  <  2:239826/143‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                             
CGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAA  >  NC_000913/4532402‑4532675

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: