Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 573,440 C→T R24C (CGC→TGC)  rusA → DLP12 prophage, endonuclease RUS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913573,4400CT100.0% 166.1 / NA 48R24C (CGC→TGC) rusADLP12 prophage, endonuclease RUS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/20);  total (28/20)

AGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCC  >  NC_000913/573305‑573575
                                                                                                                                       |                                                                                                                                       
aGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGtgca                                                                                                                                      >  2:79603/1‑139 (MQ=255)
  aaTGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGtgcacg                                                                                                                                    <  2:636749/139‑1 (MQ=255)
               aGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCag                                                                                                                       <  2:667356/139‑1 (MQ=255)
                          gaACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGggcaggca                                                                                                            <  1:309462/139‑1 (MQ=255)
                                      tCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAAc                                                                                                >  2:951065/1‑139 (MQ=255)
                                                    ggggTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAACCAt                                                                                  >  2:867263/1‑139 (MQ=255)
                                                               ggCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCaa                                                                       <  1:1045101/139‑1 (MQ=255)
                                                                 cGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATg                                                                     >  1:524214/1‑139 (MQ=255)
                                                                   tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGa                                                                                       <  1:184347/119‑1 (MQ=255)
                                                                   tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGa                                                                                       >  2:184347/1‑119 (MQ=255)
                                                                   tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                                   >  1:999914/1‑139 (MQ=255)
                                                                   tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                                   >  1:600379/1‑139 (MQ=255)
                                                                     aaTACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAAt                                                                                     >  1:17172/1‑119 (MQ=255)
                                                                     aaTACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAAt                                                                                     <  2:17172/119‑1 (MQ=255)
                                                                       tcccTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGat                                                               <  2:1082954/137‑1 (MQ=255)
                                                                             cAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                                         <  1:124075/139‑1 (MQ=255)
                                                                              aGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCCCGTCAGCGCAGCGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTCGATATCGGcc                                                        >  2:539027/1‑139 (MQ=255)
                                                                              aGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGcc                                                        >  1:91639/1‑139 (MQ=255)
                                                                               gCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCt                                                       >  1:45374/1‑139 (MQ=255)
                                                                                cATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTg                                                      >  2:490743/1‑139 (MQ=255)
                                                                                 aTCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTgg                                                     >  1:101542/1‑139 (MQ=255)
                                                                                  tCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGc                                                    >  2:503623/1‑139 (MQ=255)
                                                                                   cacaTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                                   <  2:342841/123‑1 (MQ=255)
                                                                                   cacaTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                                   >  1:342841/1‑123 (MQ=255)
                                                                                   cacaTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCt                                                   >  1:501176/1‑139 (MQ=255)
                                                                                   cacaTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCt                                                   >  1:28092/1‑139 (MQ=255)
                                                                                   cacaTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCt                                                   >  2:1070718/1‑139 (MQ=255)
                                                                                   cacaTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCt                                                   >  1:775545/1‑139 (MQ=255)
                                                                                       ttACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGc                                               <  2:924106/139‑1 (MQ=255)
                                                                                              ggCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGaa                                        <  1:704211/139‑1 (MQ=255)
                                                                                               gCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                                        <  2:769082/106‑1 (MQ=255)
                                                                                               gCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                                        >  1:769082/1‑106 (MQ=255)
                                                                                               gCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGaaa                                       >  2:1111336/1‑139 (MQ=255)
                                                                                                 cTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAAt                                     <  1:977915/139‑1 (MQ=255)
                                                                                                    cGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCg                                  >  2:972804/1‑139 (MQ=255)
                                                                                                               cGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGcc                       <  1:79603/139‑1 (MQ=255)
                                                                                                                 cTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATc                                                                                    >  2:627775/1‑76 (MQ=255)
                                                                                                                 cTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATc                                                                                    <  1:627775/76‑1 (MQ=255)
                                                                                                                   aTTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                                   <  2:1066981/91‑1 (MQ=255)
                                                                                                                   aTTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                                   >  1:1066981/1‑91 (MQ=255)
                                                                                                                     tACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                                        <  2:96705/84‑1 (MQ=255)
                                                                                                                     tACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                                        >  1:96705/1‑84 (MQ=255)
                                                                                                                        cGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATa                                                                                                          >  2:82440/1‑47 (MQ=255)
                                                                                                                        cGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATa                                                                                                          <  1:82440/47‑1 (MQ=255)
                                                                                                                          ccATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAt            >  1:343495/1‑139 (MQ=255)
                                                                                                                            aTAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                                         <  1:29165/92‑1 (MQ=255)
                                                                                                                            aTAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                                         >  2:29165/1‑92 (MQ=255)
                                                                                                                              aaTCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAg                           >  2:358618/1‑120 (MQ=255)
                                                                                                                              aaTCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAg                           <  1:358618/120‑1 (MQ=255)
                                                                                                                                    gggTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAtcgccgtcgcc  <  2:999914/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
AGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCC  >  NC_000913/573305‑573575

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