Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,582,333 T→C D239G (GAC→GGC)  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,3330TC100.0% 129.0 / NA 38D239G (GAC→GGC) hsdMtype I restriction enzyme EcoKI methyltransferase component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/24);  total (14/24)

GGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTC  >  NC_000913/4582199‑4582460
                                                                                                                                      |                                                                                                                               
ggTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAgg                                                                                                                            <  2:910058/140‑1 (MQ=255)
  tGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGtt                                                                                                                          <  2:148832/140‑1 (MQ=255)
  tGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGtt                                                                                                                          <  2:194319/140‑1 (MQ=255)
  tGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGtt                                                                                                                          <  1:242199/140‑1 (MQ=255)
  tGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGtt                                                                                                                          <  1:595295/140‑1 (MQ=255)
   gAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTg                                                                                                                         >  2:861794/1‑140 (MQ=255)
    aaCAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGc                                                                                                                        <  1:735532/140‑1 (MQ=255)
           gTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAATGGCGACAATATGCGCCCTCGGCAGGTTTTCACCCGCGGCACCCAGAGAGTTGTCCCGAAGGGTTGCGGCGCCGTGGCCGAGGTTGACTTCaaa                                                                                                                 >  2:735532/1‑139 (MQ=255)
              cGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGGTAGTGGCGACAACTGGTGCCTTCCGGACGTTTTCTCCCGCGCCGACCACGGGGTTGCCCAGACGGGTTGCCCCCCCGTGGCCGAGGGTGCCTTCAATATc                                                                                                              >  2:1082545/1‑140 (MQ=255)
                   aaTGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTgca                                                                                                         <  2:87979/140‑1 (MQ=255)
                        tGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAgcagg                                                                                                    <  1:861794/140‑1 (MQ=255)
                           tGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAgcaggcag                                                                                                 <  1:335880/140‑1 (MQ=255)
                               tGCGGCGCTGCCAAACGGCGGGTTAGCGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATc                                                                                                              >  2:510818/1‑123 (MQ=255)
                               tGCGGCGCTGCCAAACGGCGGGTTAGCGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATc                                                                                                              <  1:510818/123‑1 (MQ=255)
                                              cggcggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAgacgac                                                                              <  2:97470/140‑1 (MQ=255)
                                                 cggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGGGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGAAGGCAGTTCATCAGTGCCAGACGACggg                                                                           >  2:262674/1‑140 (MQ=255)
                                                   ggTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTg                                                                         >  1:182759/1‑140 (MQ=255)
                                                   ggTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTg                                                                         <  1:1084502/140‑1 (MQ=255)
                                                    gTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTCCCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGc                                                                        <  2:98513/140‑1 (MQ=255)
                                                           gCGACAATAAGCGCCTTCGGCAGGGTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGGTTGCGCCGCCGCGGCCGAGGTTGCCCTCAATATCGTGCAGCAGGCAGTTCCTCAGTGGCCGACGACGGGGGCCGGGCAc                                                                 >  2:1084502/1‑140 (MQ=255)
                                                              aCAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAg                                                              >  1:325312/1‑140 (MQ=255)
                                                                            cGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATaa                                                >  2:258567/1‑140 (MQ=255)
                                                                            cGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATaa                                                <  1:258567/140‑1 (MQ=255)
                                                                            cGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTAGAGGCCGATaa                                                >  2:1093691/1‑140 (MQ=255)
                                                                            cGGCAGGTTTTCACCGTCGCTACCCAGAGCGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATaa                                                >  2:1043035/1‑140 (MQ=255)
                                                                                    tttCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCgg                                        <  1:1043035/140‑1 (MQ=255)
                                                                                            tCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTg                                >  1:616716/1‑140 (MQ=255)
                                                                                                      gagTGTTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGc                      <  2:616716/140‑1 (MQ=255)
                                                                                                          gtTGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgtgt                  >  1:972757/1‑140 (MQ=255)
                                                                                                           ttGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTc                 >  2:577689/1‑140 (MQ=255)
                                                                                                           ttGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTc                 <  2:607065/140‑1 (MQ=255)
                                                                                                            tGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTg                       >  1:720549/1‑133 (MQ=255)
                                                                                                            tGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTg                       <  2:720549/133‑1 (MQ=255)
                                                                                                            tGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgtgt                  <  1:1044563/138‑1 (MQ=255)
                                                                                                            tGCCCAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCg                >  1:956365/1‑140 (MQ=255)
                                                                                                            tGCCCAGACGGATTGCGCCGCCGTAGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgtgt                  >  2:1044563/1‑138 (MQ=255)
                                                                                                              cccAGACGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGcc              >  1:570147/1‑140 (MQ=255)
                                                                                                                   aCGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCaa         <  2:485177/140‑1 (MQ=255)
                                                                                                                   aCGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCaa         >  1:862432/1‑140 (MQ=255)
                                                                                                                    cGGATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAg        <  2:723132/140‑1 (MQ=255)
                                                                                                                     ggATTGCGCCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTAGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGTGTGTCGCCATCAAgg       <  1:1093691/140‑1 (MQ=255)
                                                                                                                          gcgcCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtcgtc  <  2:956365/140‑1 (MQ=255)
                                                                                                                          gcgcCGCCGTGGCCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtcgtc  <  1:547769/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                               
GGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTC  >  NC_000913/4582199‑4582460

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