Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,149,522	C→A	Q248H (CAG→CAT)	eptC ←	phosphoethanolamine transferase EptC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,149,522	0	C	A	100.0%	120.6 / NA	36	Q248H (CAG→CAT)	eptC	phosphoethanolamine transferase EptC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/16);  total (20/16)

GTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGCTGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTG  >  NC_000913/4149386‑4149656                                                                                                                                         |                                                                                                                                       gTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATgg                                                                                                                                     >  2:521108/1‑140 (MQ=255) gTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATgg                                                                                                                                     >  1:275678/1‑140 (MQ=255)   tGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGt                                                                                                                                   <  1:79469/140‑1 (MQ=255)             cAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCa                                                                                                                         <  2:194016/140‑1 (MQ=255)                  gtgtACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCAc                                                                                                                    <  2:32610/140‑1 (MQ=255)                  gtgtACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCAc                                                                                                                    <  2:114913/140‑1 (MQ=255)                  gtgtACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCAc                                                                                                                    <  1:387662/140‑1 (MQ=255)                    gtACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCa                                                                                                                  >  2:91777/1‑140 (MQ=255)                          aCGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTa                                                                                                            >  2:382181/1‑140 (MQ=255)                               aCGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTg                                                                                                       >  2:562443/1‑140 (MQ=255)                               aCGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTg                                                                                                       >  2:537795/1‑140 (MQ=255)                                        cGTTATTGAACGCGGTCAGATTCGGATCTGTTTTATGCAGCGCATCCAGCTCCGGCGGGGTTGCACGCGGATAACCGGACGGTGTGATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCCCTAAAGTGCGCGGTTCg                                                                                              <  2:509221/140‑1 (MQ=255)                                                  cacGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGAtt                                                                                    >  1:476423/1‑140 (MQ=255)                                                       tCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCt                                                                               >  2:130620/1‑140 (MQ=255)                                                       tCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCt                                                                               >  1:167588/1‑140 (MQ=255)                                                       tCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCt                                                                               >  2:304170/1‑140 (MQ=255)                                                        cAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCtt                                                                              <  2:414175/140‑1 (MQ=255)                                                           aTTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGa                                                                           >  1:47053/1‑140 (MQ=255)                                                                     gTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAg                                                                 >  2:448207/1‑140 (MQ=255)                                                                     gTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAg                                                                 >  2:504964/1‑140 (MQ=255)                                                                           tGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTggcgg                                                           >  1:516512/1‑140 (MQ=255)                                                                            gCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTA‑CCGATTCATCTTTGAAATTAGCCAGTGGCGGCa                                                         <  1:214070/140‑1 (MQ=255)                                                                                 gcATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGc                                                     <  1:278431/140‑1 (MQ=255)                                                                                  cATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCa                                                    <  1:427611/140‑1 (MQ=255)                                                                                             ggCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTttcatt                                         >  1:307357/1‑140 (MQ=255)                                                                                                             gATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGc                         <  2:516512/140‑1 (MQ=255)                                                                                                             gATAACCGTACAGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGc                         >  2:344998/1‑140 (MQ=255)                                                                                                                        aGACTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGAGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTg              >  2:317092/1‑140 (MQ=255)                                                                                                                          aCTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCt            >  2:364493/1‑140 (MQ=255)                                                                                                                          aCTCATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCt            <  1:130620/140‑1 (MQ=255)                                                                                                                             cATGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGAc         >  1:492388/1‑140 (MQ=255)                                                                                                                              aTGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACg        <  2:275678/140‑1 (MQ=255)                                                                                                                              aTGCGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACg        <  1:562443/140‑1 (MQ=255)                                                                                                                                 cGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATa     <  1:213182/140‑1 (MQ=255)                                                                                                                                 cGTCCGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATa     <  1:387357/140‑1 (MQ=255)                                                                                                                                 cGTCCGCAATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTGCCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATa     <  2:398611/140‑1 (MQ=255)                                                                                                                                    ccGCGATGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTg  >  2:292763/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       GTTGCAGGATTTCAATGGTGTACGGACGAGACGTAACTACGTTATTAAACACGGTCAGATTCGGATCGGTTTTATGCAGCGCATCCAGCTCCGGCGTGGTTTCACGCGGATAACCGTACAGACTCATGCGTCCGCGCTGGGTCGACTCGCCAATCACCAGCACTAAAGTGCGCGGTTCGTTACCCGATTCATCTTTGAAATTAGCCAGTGGCGGCAAGGCATTATTTTCATTCAGTAACTTTGTTAGCGAGTTTAGTTGCTGACGATACTG  >  NC_000913/4149386‑4149656

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