Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,412,023	(T)7→8	coding (45/399 nt)	yjfL →	conserved inner membrane protein YjfL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,016	1	.	T	100.0%	64.4 / NA	19	Y13F (TAT→TTT)	yjfL	conserved inner membrane protein YjfL
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (12/7);  total (12/7)

CATCACTTCATATTATTGGTTAATTACAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTA‑TTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTATTCTAAAATTACACCGCACAACGAATGGC  >  NC_000913/4411939‑4412090                                                                               |                                                                           cATCACTTCATATTATTGGTTAATTACAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTA‑ttt                                                                         <  2:888905/81‑1 (MQ=255)     aCTTCATATTATTGGTTAATTACAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTA‑ttttttt                                                                     >  2:888940/1‑81 (MQ=255)         cATATTATTGGTTAATTACAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTA‑tttttttta                                                                   <  1:888907/79‑3 (MQ=255)         cATATTATTGGTTAATTACAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTA‑tttttttta                                                                   <  2:888906/79‑3 (MQ=255)                ttGGTTAATTACAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTgg                                                           >  1:888941/1‑81 (MQ=255)                    ttAATTACAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTAt                                                       <  2:888909/81‑1 (MQ=255)                    ttAATTACAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTAt                                                       <  2:888908/81‑1 (MQ=255)                          aCAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGat                                                 >  2:888937/1‑81 (MQ=255)                           cAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGat                                                 <  1:888910/80‑1 (MQ=255)                              ggACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGattatt                                             >  1:888942/1‑81 (MQ=255)                                aCGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTAtttt                                           <  2:888912/81‑1 (MQ=255)                                aCGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTAtttt                                           >  2:888938/1‑81 (MQ=255)                                aCGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTAtttt                                           <  1:888911/81‑1 (MQ=255)                                           gCATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATc                                 >  2:888952/1‑80 (MQ=255)                                            cATATACTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTa                               >  1:888944/1‑81 (MQ=255)                                                  cTGGATTCACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTATTCTa                          >  1:888955/1‑80 (MQ=255)                                                         cACTTCTTGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTATTCTAAAATTaca                  >  2:888949/1‑81 (MQ=255)                                                             tcttGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTATTCTAAAATTACACCg               >  2:888946/1‑80 (MQ=255)                                                                tGCCTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTATTCTAAAATTACACCGcac            >  2:888947/1‑80 (MQ=255)                                                                  ccTTTAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTATTCTAAAATTACACCGCACAAc         <  1:888913/81‑1 (MQ=255)                                                                    tttAGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGCTTATCTATTCTAAAATTACACCGCACAACg        >  2:888943/1‑80 (MQ=255)                                                                       aGCGCCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTATTCTAAAATTACACCGCACAACGAATg    >  1:888945/1‑81 (MQ=255)                                                                         cgcCTATTTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTATTCTAAAATTACACCGCACAACGAATGGc  <  1:888914/81‑1 (MQ=255)                                                                               |                                                                           CATCACTTCATATTATTGGTTAATTACAAGGACGTTTACAATGCATATACTGGATTCACTTCTTGCCTTTAGCGCCTA‑TTTTTTTATTGGCGTGGCTATGGTGATTATTTTTCTGTTTATCTATTCTAAAATTACACCGCACAACGAATGGC  >  NC_000913/4411939‑4412090

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