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%	171.6 / NA	48	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 (21/27);  total (21/27)

GGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCC  >  NC_000913/631686‑631938                                                                                                                               |                                                                                                                               ggCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCt                                                                                                                    <  1:1431962/139‑1 (MQ=255)       aaGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGc                                                                                                              <  2:679955/139‑1 (MQ=255)       aaGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGc                                                                                                              <  2:796464/139‑1 (MQ=255)            tttAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGt                                                                                                         >  2:30783/1‑139 (MQ=255)                 cccGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAAcc                                                                                                    <  1:1686162/139‑1 (MQ=255)                   cGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                  <  1:208925/139‑1 (MQ=255)                   cGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                  <  1:74956/139‑1 (MQ=255)                                               tGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTg                                                                      >  2:819403/1‑139 (MQ=255)                                                  taataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                   <  2:346611/139‑1 (MQ=255)                                                  taataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                   <  1:975777/139‑1 (MQ=255)                                                    ataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                   <  1:194553/137‑1 (MQ=255)                                                    ataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                   >  2:194553/1‑137 (MQ=255)                                                    ataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACt                                                                 <  2:601288/139‑1 (MQ=255)                                                       aGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGtt                                                              <  1:30783/139‑1 (MQ=255)                                                           ccAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCg                                                          >  2:150148/1‑139 (MQ=255)                                                           ccAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCg                                                          >  2:21892/1‑139 (MQ=255)                                                                  aTGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGa                                                   <  1:819403/139‑1 (MQ=255)                                                                    gATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACg                                                 <  1:1266720/139‑1 (MQ=255)                                                                           aGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                  <  1:1812490/83‑1 (MQ=255)                                                                           aGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                  >  2:1812490/1‑83 (MQ=255)                                                                           aGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                          <  1:936909/139‑1 (MQ=255)                                                                           aGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                          <  1:150148/139‑1 (MQ=255)                                                                             cGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCcaacaa                                                                                                                          >  1:236625/1‑57 (MQ=255)                                                                             cGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCcaacaa                                                                                                                          <  2:236625/57‑1 (MQ=255)                                                                             cGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCg                                        >  2:1381910/1‑139 (MQ=255)                                                                               gCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGc                                      >  1:321607/1‑139 (MQ=255)                                                                               gCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGc                                      >  1:1448669/1‑139 (MQ=255)                                                                                          tCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATg                                                                                                               <  1:300930/55‑1 (MQ=255)                                                                                          tCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATg                                                                                                               >  2:300930/1‑55 (MQ=255)                                                                                          tCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATc                           >  2:50413/1‑139 (MQ=255)                                                                                                       aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                  >  2:520128/1‑55 (MQ=255)                                                                                                       aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                  <  1:520128/55‑1 (MQ=255)                                                                                                       aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGc              >  2:937250/1‑139 (MQ=255)                                                                                                       aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGc              <  1:21892/139‑1 (MQ=255)                                                                                                         tCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCga            <  1:428155/139‑1 (MQ=255)                                                                                                          cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                   <  1:82098/83‑1 (MQ=255)                                                                                                          cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                   >  2:82098/1‑83 (MQ=255)                                                                                                          cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGc                                               <  2:1399850/103‑1 (MQ=255)                                                                                                          cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGc                                               >  1:1399850/1‑103 (MQ=255)                                                                                                            caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaaag         >  2:1722172/1‑139 (MQ=255)                                                                                                             aGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaaaga        <  1:168316/139‑1 (MQ=255)                                                                                                               cTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                          >  1:585844/1‑103 (MQ=255)                                                                                                               cTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                          <  2:585844/103‑1 (MQ=255)                                                                                                                 ggCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAtctt                                                                                               <  2:390476/48‑1 (MQ=255)                                                                                                                 ggCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAtctt                                                                                               >  1:390476/1‑48 (MQ=255)                                                                                                                  gcgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGc   <  2:321607/139‑1 (MQ=255)                                                                                                                   cgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGcc  >  2:1830093/1‑139 (MQ=255)                                                                                                                   cgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTCGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGcc  >  2:321717/1‑139 (MQ=255)                                                                                                                               |                                                                                                                               GGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCC  >  NC_000913/631686‑631938

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: ‑