Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,682,286 G→A R43Q (CGA→CAA)  rstA → DNA‑binding transcriptional regulator RstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,682,2860GA100.0% 61.1 / NA 18R43Q (CGA→CAA) rstADNA‑binding transcriptional regulator RstA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/6);  total (12/6)

GCCGCGTACCTGGCAAAACATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGC  >  NC_000913/1682210‑1682362
                                                                            |                                                                            
gCCGCGTACCTGGCAAAACATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCaaga                                                                           <  2:307424/80‑1 (MQ=255)
               aaaCATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGgtg                                                           >  1:307443/1‑81 (MQ=255)
                           cAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGAcat                                                >  2:307452/1‑80 (MQ=255)
                            aGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGAcatca                                              <  2:307425/81‑1 (MQ=255)
                                tACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGc                                           >  2:307447/1‑80 (MQ=255)
                                 aCCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTa                                         >  1:307451/1‑81 (MQ=255)
                                  ccGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTAc                                        <  2:307426/81‑1 (MQ=255)
                                   cGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTAc                                        >  1:307450/1‑80 (MQ=255)
                                            gcgcgGCGACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGa                              <  2:307427/81‑1 (MQ=255)
                                             cgcgGCGACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGAc                             >  1:307453/1‑81 (MQ=255)
                                                   gACCAGGCCGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATg                       <  2:307428/81‑1 (MQ=255)
                                                           cGAAGAAACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTg               >  2:307459/1‑81 (MQ=255)
                                                              agaaACCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGTACGGCATGACCATTTGTCg            >  2:307448/1‑81 (MQ=255)
                                                                aaaCCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGt           >  2:307456/1‑80 (MQ=255)
                                                                 aaCCATTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGa         >  1:307449/1‑81 (MQ=255)
                                                                     aTTTTGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTa     >  1:307441/1‑81 (MQ=255)
                                                                        ttGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTAcgc  >  2:307455/1‑81 (MQ=255)
                                                                         tGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTAcgc  <  2:307430/80‑1 (MQ=255)
                                                                         tGCAAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTAcgc  <  1:307429/80‑1 (MQ=255)
                                                                            |                                                                            
GCCGCGTACCTGGCAAAACATGATATGCAGGTTACCGTAGAGCCGCGCGGCGACCAGGCCGAAGAAACCATTTTGCGAGAAAATCCGGATTTGGTGTTACTCGACATCATGCTACCAGGCAAGGACGGCATGACCATTTGTCGTGATTTACGC  >  NC_000913/1682210‑1682362

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