Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,162,554	A→G	L39P (CTA→CCA)	ftsP ←	cell division protein required during stress conditions

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,162,554	0	A	G	100.0%	143.5 / NA	40	L39P (CTA→CCA)	ftsP	cell division protein required during stress conditions
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/19);  total (21/19)

GCGCGTGTCCCTGGCGTAAATGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTAGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAGCGCCTGCACAAAGT  >  NC_000913/3162476‑3162625                                                                               |                                                                        gcgcgTGTCCCTGGCGTAAATGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGgc                                                                       >  1:551732/1‑81 (MQ=255)  cgcgTGTCCCTGGCGTAAATGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGgcg                                                                      >  2:551755/1‑81 (MQ=255)   gcgTGTCCCTGGCGTAAATGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGgcgg                                                                     >  2:551770/1‑81 (MQ=255)         cccTGGCGTAAATGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAAc                                                               >  2:551761/1‑81 (MQ=255)          ccTGGCGTAAATGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAAc                                                               <  2:551708/80‑1 (MQ=255)             ggCGTAAATGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCTACGGCGAGATTCAAGTGGCGGCGGAACggg                                                            >  2:551752/1‑80 (MQ=255)               cGTAAATGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAg                                                         >  1:551757/1‑81 (MQ=255)                  aaaTGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCgg                                                      >  1:551765/1‑81 (MQ=255)                   aaTGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGt                                                     >  1:551764/1‑81 (MQ=255)                   aaTGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGt                                                     <  2:551709/81‑1 (MQ=255)                    aTGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGtt                                                    >  2:551760/1‑81 (MQ=255)                    aTGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGtt                                                    <  1:551710/81‑1 (MQ=255)                         cAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGt                                               >  2:551772/1‑81 (MQ=255)                         cAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGt                                               >  1:551766/1‑81 (MQ=255)                             gcgcACGTTGTACAGTCATAAACAGGGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGcc                                           >  2:551763/1‑81 (MQ=255)                              cgcACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGccc                                          <  2:551711/81‑1 (MQ=255)                              cgcACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGccc                                          <  1:551713/81‑1 (MQ=255)                              cgcACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGcc                                           <  2:551712/80‑1 (MQ=255)                                cACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCgg                                        >  2:551753/1‑81 (MQ=255)                                   gTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTg                                     <  1:551714/81‑1 (MQ=255)                                   gTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTg                                     >  2:551750/1‑81 (MQ=255)                                    ttGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTgc                                    >  1:551777/1‑81 (MQ=255)                                    ttGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTgc                                    <  2:551715/81‑1 (MQ=255)                                      gTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTgcgc                                  >  1:551758/1‑81 (MQ=255)                                          aGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGc                              >  1:551767/1‑81 (MQ=255)                                             cATAAACAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCtt                           >  1:551759/1‑81 (MQ=255)                                                 aaCAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAgg                       <  2:551717/81‑1 (MQ=255)                                                 aaCAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAgg                       <  1:551718/81‑1 (MQ=255)                                                 aaCAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAgg                       <  1:551716/81‑1 (MQ=255)                                                  aCAGCGGTTGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAggg                      >  2:551769/1‑81 (MQ=255)                                                         ttGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAgc               <  1:551719/81‑1 (MQ=255)                                                          tGCCCACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAgcg              >  2:551771/1‑81 (MQ=255)                                                             ccACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAGCGCCt           <  2:551721/81‑1 (MQ=255)                                                             ccACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAGCGCCt           <  1:551720/81‑1 (MQ=255)                                                              cACGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAGCGCCTg          <  2:551722/81‑1 (MQ=255)                                                               aCGGCGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAGCGCCTGc         <  1:551723/81‑1 (MQ=255)                                                               aCGGCGAGAGTCGAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGGACAGCGCCTGc         <  2:551724/81‑1 (MQ=255)                                                                   cGAGATTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAGCGCCTGCACaa     <  2:551725/81‑1 (MQ=255)                                                                      gaTTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAGCGCCTGCACAAAGt  <  2:551726/81‑1 (MQ=255)                                                                      gaTTCAAGTGGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAGCGCCTGCACAAAGt  >  2:551775/1‑81 (MQ=255)                                                                               |                                                                        GCGCGTGTCCCTGGCGTAAATGACCAGTGCGCACGTTGTACAGTCATAAACAGCGGTTGCCCACGGCGAGATTCAAGTAGCGGCGGAACGGGTAGCGGTTGCTGTTGCCCGGCTGCGCTGGCCTTCAGGGGAACAGCGCCTGCACAAAGT  >  NC_000913/3162476‑3162625

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