Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	631,812	T→C	F640S (TTC→TCC)	cstA →	carbon starvation protein involved in peptide utilization, APC peptide transporter family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	631,812	0	T	C	100.0%	114.4 / NA	33	F640S (TTC→TCC)	cstA	carbon starvation protein involved in peptide utilization, APC peptide transporter family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/15);  total (18/15)

CGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGC  >  NC_000913/631678‑631937                                                                                                                                       |                                                                                                                              cGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCcaac                                                                                                                           <  1:906182/139‑1 (MQ=255)                    tttAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTACTAAGTTCCAGGCAATGATCGACAGCGGCAATTTTCCGTCGCAGTATACTGAGTCACAGATGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGt                                                                                                        >  1:933766/1‑139 (MQ=255)                    tttAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGt                                                                                                        <  1:556908/139‑1 (MQ=255)                                              ccTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 <  2:745557/120‑1 (MQ=255)                                              ccTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 >  1:745557/1‑120 (MQ=255)                                                       tGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 >  1:1250670/1‑111 (MQ=255)                                                       tGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 <  2:1250670/111‑1 (MQ=255)                                                                 ttCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGTGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTtct                                                           >  1:234645/1‑139 (MQ=255)                                                                 ttCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTtct                                                           >  1:612163/1‑139 (MQ=255)                                                                   ccAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 <  2:898565/99‑1 (MQ=255)                                                                   ccAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 >  1:898565/1‑99 (MQ=255)                                                                     cgacggagATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 <  2:854120/90‑1 (MQ=255)                                                                     cgacggagATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 >  1:854120/8‑97 (MQ=255)                                                                      ggCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGAtt                                                      <  2:933766/139‑1 (MQ=255)                                                                                       gCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 >  1:1409978/1‑79 (MQ=255)                                                                                       gCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 <  2:1409978/79‑1 (MQ=255)                                                                                       gCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGc                                     >  1:748242/1‑139 (MQ=255)                                                                                             ttCCGTCGCAGTATACTGAGTCACAGCAGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTATTTATGGTGGTCGTGGTGGTTCTGGCACTGATCTCGATTAAGACGGCACTTGCGGCATTGaa                               >  2:182995/1‑139 (MQ=255)                                                                                                         ataCTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGcc                   <  1:182995/139‑1 (MQ=255)                                                                                                            cTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGt                                                                                                        >  1:602207/1‑51 (MQ=255)                                                                                                            cTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGt                                                                                                        <  2:602207/51‑1 (MQ=255)                                                                                                               aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 <  1:66512/55‑1 (MQ=255)                                                                                                               aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                 >  2:66512/1‑55 (MQ=255)                                                                                                                gTCACAGCTGGCGCAACT‑GTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGAc                                                 >  1:437861/1‑101 (MQ=255)                                                                                                                gTCACAGCTGGCGCAACT‑GTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGAc                                                 <  2:437861/101‑1 (MQ=255)                                                                                                                    caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGc                                     >  2:525818/1‑110 (MQ=255)                                                                                                                    caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGc                                     <  1:525818/110‑1 (MQ=255)                                                                                                                    caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaaag        >  2:726043/1‑139 (MQ=255)                                                                                                                      gCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaaagaa      >  2:61294/1‑139 (MQ=255)                                                                                                                      gCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaaagaa      >  1:1290932/1‑139 (MQ=255)                                                                                                                         ggtgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACg   <  2:234645/136‑1 (MQ=255)                                                                                                                         ggCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACg   >  2:876622/1‑139 (MQ=255)                                                                                                                          gcgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGc  <  1:61294/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                              CGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGC  >  NC_000913/631678‑631937

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