Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,423,484 T→A pseudogene (284/444 nt)
C29S (TGC→AGC) 
rzpR →
rzoR →
pseudogene, Rac prophage, Bacteriophage Rz lysis protein family,Phage or Prophage Related, putative Rac prophage endopeptidase, completely contained in another CDS
Rac prophage, putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,423,4840TA100.0% 71.8 / NA 26pseudogene (284/444 nt)
C29S (TGC→AGC) 
rzpR
rzoR
pseudogene, Rac prophage, Bacteriophage Rz lysis protein family,Phage or Prophage Related, putative Rac prophage endopeptidase, completely contained in another CDS
Rac prophage, putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (12/14);  total (12/14)

CTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATA  >  NC_000913/1423369‑1423619
                                                                                                                   |                                                                                                                                       
cTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCg                                                                                                                  >  2:237167/1‑139 (MQ=17)
 tCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCgg                                                                                                                 <  1:672731/139‑1 (MQ=18)
      gCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTgg                                                                                                            <  2:631855/139‑1 (MQ=18)
             acacGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGc                                                                                                     >  1:1082441/1‑139 (MQ=17)
                        ttAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTcccccc                                                                                          >  2:618041/1‑139 (MQ=17)
                           gCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGAc                                                                                       >  2:506560/1‑139 (MQ=18)
                                                 atgCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGGCTCCGGCGTGGATAATGCAGCCTCCCGCCGACTGGCAGACACCGCTGAACGGGa                                                                 <  1:677838/139‑1 (MQ=12)
                                                   gCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCt                                                                                                                     >  2:557261/1‑85 (MQ=18)
                                                   gCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCt                                                                                                                     <  1:557261/85‑1 (MQ=25)
                                                              tgatgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                                    >  2:224812/1‑139 (MQ=18)
                                                                 tgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCgg                                                 >  2:178165/1‑139 (MQ=25)
                                                                  gTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCggg                                                <  2:122205/139‑1 (MQ=37)
                                                                      ccGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACg                                            <  2:1163491/139‑1 (MQ=255)
                                                                            gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACg                                            <  1:416113/133‑1 (MQ=255)
                                                                            gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACg                                            >  2:416113/1‑133 (MQ=255)
                                                                             tcgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGa                                                                                        >  2:225105/1‑88 (MQ=18)
                                                                             tcgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGa                                                                                        <  1:225105/88‑1 (MQ=18)
                                                                                                      gTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTg            <  2:1082441/139‑1 (MQ=255)
                                                                                                      gTCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTg            <  1:237167/139‑1 (MQ=255)
                                                                                                       tCTGTCAGTCAGAGCTTGAAGCCACCACCTCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGaa                                              <  1:905916/104‑1 (MQ=255)
                                                                                                       tCTGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGaa                                              >  2:905916/1‑104 (MQ=255)
                                                                                                         tGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACg                                            >  2:388230/1‑104 (MQ=255)
                                                                                                         tGTCAGTCAGAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACg                                            <  1:388230/104‑1 (MQ=255)
                                                                                                           tcagtcagAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTcccccc                                                                                          >  2:473325/1‑56 (MQ=255)
                                                                                                           tcagtcagAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTcccccc                                                                                          <  1:473325/56‑1 (MQ=255)
                                                                                                           tcagtcagAGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAg       <  1:224812/139‑1 (MQ=255)
                                                                                                                cagaGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATa  <  1:721210/139‑1 (MQ=255)
                                                                                                                   |                                                                                                                                       
CTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATA  >  NC_000913/1423369‑1423619

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