Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,582,392 G→C 10.6% P219P (CCC→CCG)  hsdM ← type I restriction enzyme EcoKI methyltransferase component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,582,3920GC10.6% 143.2 / 8.1 47P219P (CCC→CCG) hsdMtype I restriction enzyme EcoKI methyltransferase component
Reads supporting (aligned to +/- strand):  ref base G (18/24);  new base C (3/2);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTG  >  NC_000913/4582247‑4582535
                                                                                                                                                 |                                                                                                                                               
gcggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCa                                                                                                                                             >  1:566011/1‑150 (MQ=255)
 cggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCAc                                                                                                                                            <  1:472778/150‑1 (MQ=255)
  gggTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCAcc                                                                                                                                           <  2:484795/150‑1 (MQ=255)
   ggTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCa                                                                                                                                          >  1:711850/1‑150 (MQ=255)
    gTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAg                                                                                                                                         >  1:50519/1‑150 (MQ=255)
         tGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCCGGCACCCGtcgtc                                                                                                                                    <  1:42312/150‑5 (MQ=255)
          ggCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGTGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAg                                                                                                                                   >  1:227153/1‑150 (MQ=255)
           gCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGGCGAGGTGGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAgg                                                                                                                                  <  2:749453/150‑1 (MQ=255)
                    tGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAAc                                                                                                                         >  2:985364/1‑150 (MQ=255)
                             ggCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGATGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTgg                                                                                                                <  2:1421915/150‑1 (MQ=255)
                                  gTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTg                                                                                                           <  1:106606/150‑1 (MQ=255)
                                    tttCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGa                                                                                                         <  2:907370/150‑1 (MQ=255)
                                          cGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCt                                                                                                   <  2:924889/150‑1 (MQ=255)
                                             cGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCg                                                                                                <  1:1208648/150‑1 (MQ=255)
                                             cGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCg                                                                                                <  2:566011/150‑1 (MQ=255)
                                                   ccAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCCGGCa                                                                                                                                             >  2:48103/1‑99 (MQ=255)
                                                   ccAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCCGGCa                                                                                                                                             <  1:48103/99‑1 (MQ=255)
                                                     agagTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCa                                                                                        >  2:1201521/1‑150 (MQ=255)
                                                            tGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgt                                                                                 >  1:1091570/1‑150 (MQ=255)
                                                             gCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtc                                                                                >  2:390938/1‑150 (MQ=255)
                                                              cccAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtcg                                                                               <  1:401119/150‑1 (MQ=255)
                                                                cAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtcgtc                                                                             >  1:464511/1‑150 (MQ=255)
                                                                 aGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTcc                                                                            >  1:241496/1‑150 (MQ=255)
                                                                  gACGGATTGCGCCGCCGTGGTCGAGGTGGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCa                                                                           <  2:565036/150‑1 (MQ=255)
                                                                        ttGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCAt                                                                     <  2:337147/150‑1 (MQ=255)
                                                                         tGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCAtt                                                                    >  1:175713/1‑150 (MQ=255)
                                                                          gcgcCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGGCGTCCAGATCATTg                                                                   <  1:980670/150‑1 (MQ=255)
                                                                                ccgTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGc                                                             <  2:510165/150‑1 (MQ=255)
                                                                                        gAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAAc                                                     >  1:561105/1‑150 (MQ=255)
                                                                                            ttGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAg                                                 <  2:227153/150‑1 (MQ=255)
                                                                                             tGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGc                                                <  1:985364/150‑1 (MQ=255)
                                                                                                    aaTATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGc                                         <  2:175713/150‑1 (MQ=255)
                                                                                                         cGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtcaa                                    >  2:895926/1‑150 (MQ=255)
                                                                                                           tGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtcaatc                                  <  2:787716/150‑1 (MQ=255)
                                                                                                            gcaccaGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACAt                                                   >  1:1180373/1‑132 (MQ=255)
                                                                                                            gcaccaGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACAt                                                   <  2:1180373/132‑1 (MQ=255)
                                                                                                             cagcagGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtcaatcaa                                >  2:879792/1‑150 (MQ=255)
                                                                                                               gcaggcagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCaaaa                              <  1:1154778/150‑1 (MQ=255)
                                                                                                                 aggcagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAg                            <  1:390938/150‑1 (MQ=255)
                                                                                                                      gTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCccgc                       >  1:745349/1‑150 (MQ=255)
                                                                                                                         catcaGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGt                    >  2:689442/1‑150 (MQ=255)
                                                                                                                            caGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTAcc                 <  2:1091570/150‑1 (MQ=255)
                                                                                                                              gTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTg               >  1:1230023/1‑150 (MQ=255)
                                                                                                                                     cgaactggTGCCCGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCgg        >  2:358625/7‑150 (MQ=255)
                                                                                                                                      cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCggg       <  2:464511/150‑1 (MQ=255)
                                                                                                                                       aactggTGCCCGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGt      >  1:207707/5‑150 (MQ=255)
                                                                                                                                           ggTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTg  <  2:241496/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTG  >  NC_000913/4582247‑4582535

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