Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,123,764	G→A	A231A (GCC→GCT)	cytR ←	DNA‑binding transcriptional repressor CytR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,123,764	0	G	A	100.0%	155.9 / NA	44	A231A (GCC→GCT)	cytR	DNA‑binding transcriptional repressor CytR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/29);  total (15/29)

GCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAA  >  NC_000913/4123687‑4123837                                                                              |                                                                          gcCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCt                                                                        >  1:767692/1‑81 (MQ=255)  ccATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCtt                                                                       <  1:767635/81‑1 (MQ=255)     tCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCg                                                                     >  1:767712/1‑80 (MQ=255)      cacaTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCg                                                                     >  1:767680/1‑79 (MQ=255)        caTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAgg                                                                 <  1:767636/81‑1 (MQ=255)        caTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAgg                                                                 >  2:767700/1‑81 (MQ=255)           cGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTg                                                               >  2:767714/1‑80 (MQ=255)            gCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTg                                                               >  1:767716/1‑79 (MQ=255)                  gCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGc                                                        >  1:767717/1‑80 (MQ=255)                   cAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCa                                                      >  2:767718/1‑81 (MQ=255)                   cAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCa                                                      <  1:767637/81‑1 (MQ=255)                    agaagaCAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCAc                                                     >  2:767705/1‑81 (MQ=255)                    agaagaCAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCAc                                                     <  2:767638/81‑1 (MQ=255)                    agaagaCAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCAc                                                     <  2:767637/81‑1 (MQ=255)                            gcagTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATg                                             <  1:767640/81‑1 (MQ=255)                            gcagTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATg                                             <  1:767641/81‑1 (MQ=255)                             cagTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGt                                            <  2:767641/81‑1 (MQ=255)                              agTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTa                                           >  1:767714/1‑81 (MQ=255)                               gTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTAt                                          >  1:767720/1‑81 (MQ=255)                                  ggCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGc                                       <  2:767642/81‑1 (MQ=255)                                  ggCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGc                                       <  2:767643/81‑1 (MQ=255)                                      gTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGAt                                   <  1:767645/81‑1 (MQ=255)                                         gtgGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCa                                 <  1:767646/80‑1 (MQ=255)                                            gAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAAcc                              <  1:767647/80‑1 (MQ=255)                                            gAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCa                             <  1:767649/81‑1 (MQ=255)                                            gAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCa                             <  2:767647/81‑1 (MQ=255)                                             aaGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCAt                            <  2:767649/81‑1 (MQ=255)                                             aaGATCAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCAt                            <  1:767651/81‑1 (MQ=255)                                                  cAAGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGc                       >  1:767710/1‑81 (MQ=255)                                                   aaGCAGCTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGcc                      <  2:767652/81‑1 (MQ=255)                                                   aaGCAGCTGCTGCATCGCTTTGCTTCCAGCATCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGc                       <  2:767651/80‑1 (MQ=255)                                                     gcagcTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCg                     >  1:767732/1‑80 (MQ=255)                                                      cagcTGCTGCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCa                   <  2:767653/81‑1 (MQ=255)                                                         ctgctgCATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCg                <  1:767655/81‑1 (MQ=255)                                                           gctgcATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGAc              >  2:767717/1‑81 (MQ=255)                                                            ctgcATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGa               <  2:767655/79‑1 (MQ=255)                                                            ctgcATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGAc              <  1:767657/80‑1 (MQ=255)                                                               cATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCa           <  2:767658/80‑1 (MQ=255)                                                               cATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAg          <  1:767658/81‑1 (MQ=255)                                                                   gCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAgcgc       <  1:767660/80‑1 (MQ=255)                                                                    cTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCt     <  2:767660/81‑1 (MQ=255)                                                                       tGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGaa  >  2:767712/1‑81 (MQ=255)                                                                       tGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGa   <  1:767662/80‑1 (MQ=255)                                                                        gCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGaa  <  2:767662/80‑1 (MQ=255)                                                                              |                                                                          GCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAA  >  NC_000913/4123687‑4123837

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