Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	629,717	C→T	intergenic (+4/‑177)	entH → / → cstA	proofreading thioesterase in enterobactin biosynthesis/carbon starvation protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	629,717	0	C	T	100.0%	89.5 / NA	26	intergenic (+4/‑177)	entH/cstA	proofreading thioesterase in enterobactin biosynthesis/carbon starvation protein A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/16);  total (10/16)

AAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTG  >  NC_000913/629643‑629791                                                                           |                                                                           aaaTCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAg                                                                       >  1:123861/1‑80 (MQ=255)         gtTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGAt                                                               <  1:123816/80‑1 (MQ=255)         gtTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGAt                                                               >  1:123869/1‑80 (MQ=255)            ttCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATgag                                                            <  2:123817/80‑1 (MQ=255)              cGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGAt                                                          <  2:123819/80‑1 (MQ=255)              cGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGAt                                                          <  2:123818/80‑1 (MQ=255)                aTGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTg                                                        >  1:123837/1‑80 (MQ=255)                   aaCAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAg                                                     <  1:123820/80‑1 (MQ=255)                     cAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGaa                                                  <  2:123821/81‑1 (MQ=255)                     cAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGa                                                   >  2:123862/1‑80 (MQ=255)                      aGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGaa                                                  >  1:123866/1‑80 (MQ=255)                        gggCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAAc                                               <  2:123822/81‑1 (MQ=255)                             gcgTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGt                                          >  1:123877/1‑81 (MQ=255)                              cgTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGtt                                         <  2:123823/81‑1 (MQ=255)                               gTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTa                                        <  1:123824/81‑1 (MQ=255)                                 tgctgcACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAAc                                      >  1:123874/1‑81 (MQ=255)                                   ctgcACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACg                                     <  2:123825/80‑1 (MQ=255)                                      cACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGt                                 <  2:123826/81‑1 (MQ=255)                                            tCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATc                             >  2:123868/1‑79 (MQ=255)                                            tCGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGa                           <  1:123827/81‑1 (MQ=255)                                             cGGCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAg                          >  2:123863/1‑81 (MQ=255)                                               gCTGGGTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGtt                        <  2:123828/81‑1 (MQ=255)                                                  gggTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTaa                      >  1:123870/1‑80 (MQ=255)                                                    gTACGGCAGTTTTGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACAt                   <  1:123829/81‑1 (MQ=255)                                                             ttttGGGATGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTgttaagtt          <  2:123830/81‑1 (MQ=255)                                                                     tGAACCTGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTg  <  2:123831/81‑1 (MQ=255)                                                                           |                                                                           AAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTG  >  NC_000913/629643‑629791

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