Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,851,887 G→A R479H (CGT→CAT)  mutS → DNA mismatch repair protein MutS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,851,8870GA100.0% 66.2 / NA 19R479H (CGT→CAT) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/6);  total (13/6)

CCGCGAGCGTGAACGTACCGGCCTGGACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCGTGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTACATCATTCCAGAGCT  >  NZ_CP009273/2851807‑2851968
                                                                                |                                                                                 
ccGCGAGCGTGAACGTACCGGCCTGGACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCATGGGCAAAg                                                                          >  1:353060/1‑90 (MQ=255)
ccGCGAGCGTGAACGTACCGGCCTGGACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCATGGGCAAAg                                                                          >  1:114827/1‑90 (MQ=255)
           aaCGTACCGGCCTGGACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCATGGGCAAAGCCATCTGGCAc                                                               <  1:15118/90‑1 (MQ=255)
           aaCGTACCGGCCTGGACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCATGGGCAAAGCCATCTGGCAc                                                               <  1:437794/90‑1 (MQ=255)
                 ccGGCCTGGACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCATGGGCAAAGCCATCTGGCACCCATCa                                                         >  2:138568/1‑90 (MQ=255)
                 ccGGCCTGGACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCATGGGCAAAGCCATCTGGCACCCATCa                                                         >  2:105650/1‑90 (MQ=255)
                        ggACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACAt                                                  >  2:500853/1‑90 (MQ=255)
                        ggACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACAt                                                  >  2:500305/1‑90 (MQ=255)
                                                  cGGTGCACGGCTACTACATTCAAATCAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCg                        <  1:79736/90‑1 (MQ=255)
                                                              actacATTCAAATCAGCCATGGGCAAAGCCATCTGGCAccc                                                             <  1:443517/41‑1 (MQ=255)
                                                              actacATTCAAATCAGCCATGGGCAAAGCCATCTGGCAccc                                                             >  2:443517/1‑41 (MQ=255)
                                                              actacATTCAAATCAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTAcatca            >  2:17095/1‑90 (MQ=255)
                                                              actacATTCAAATCAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTAcatca            >  2:463459/1‑90 (MQ=255)
                                                              actacATTCAAATCAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTAcatca            >  2:617100/1‑90 (MQ=255)
                                                                   aTTCAAATCAGCCATGGGCAAAGCCATCTGGCCCCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTACATCATTCCa       <  1:463459/90‑1 (MQ=255)
                                                                        aaTCAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTACATCATTCCAGAGCt  >  1:281488/1‑90 (MQ=255)
                                                                        aaTCAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTACATCATTCCAGAGCt  >  1:274005/1‑90 (MQ=255)
                                                                           cAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTACATCATTCCaga     <  2:230912/84‑1 (MQ=255)
                                                                           cAGCCATGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTACATCATTCCaga     >  1:230912/1‑84 (MQ=255)
                                                                                |                                                                                 
CCGCGAGCGTGAACGTACCGGCCTGGACACGCTGAAAGTTGGCTTTAATGCGGTGCACGGCTACTACATTCAAATCAGCCGTGGGCAAAGCCATCTGGCACCCATCAACTACATGCGTCGCCAGACGCTGAAAAACGCCGAGCGCTACATCATTCCAGAGCT  >  NZ_CP009273/2851807‑2851968

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: