Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 731,455 T→C C625R (TGT→CGT)  rhsC → rhs element protein RhsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913731,4550TC100.0% 59.3 / NA 20C625R (TGT→CGT) rhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/6);  total (14/6)

TAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTG  >  NC_000913/731324‑731520
                                                                                                                                   |                                                                 
tAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACg                                                           <  2:109753/140‑1 (MQ=255)
       cGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCg                                                    >  2:427192/1‑140 (MQ=255)
                 acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACAGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          <  2:361939/140‑1 (MQ=255)
                   gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          <  1:247255/138‑1 (MQ=255)
                   gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          >  2:5982/1‑138 (MQ=255)
                   gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          <  1:5982/138‑1 (MQ=255)
                   gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          >  2:247255/1‑138 (MQ=255)
                        gCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATg                                   >  2:402388/1‑140 (MQ=255)
                                 cgcgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          >  2:772443/1‑124 (MQ=255)
                                 cgcgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          <  1:772443/124‑1 (MQ=255)
                                 cgcgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAt                          >  2:666633/1‑140 (MQ=255)
                                   cgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAtgc                        >  1:178937/1‑140 (MQ=255)
                                                aCGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc           >  2:113862/1‑140 (MQ=255)
                                                aCGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc           >  1:355730/1‑140 (MQ=255)
                                                    cgccgGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATcc       >  2:213203/1‑140 (MQ=18)
                                                    cgccgGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATcc       >  1:561794/1‑140 (MQ=18)
                                                    cgccgGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATcc       >  2:622732/1‑140 (MQ=17)
                                                         gTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg  >  1:751322/1‑140 (MQ=18)
                                                                 accaccGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc    <  2:421935/130‑1 (MQ=17)
                                                                 accaccGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc    >  1:421935/1‑130 (MQ=18)
                                                                                                                                   |                                                                 
TAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTG  >  NC_000913/731324‑731520

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