Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,322,128	(C)5→4	coding (1149/1650 nt)	yfaQ ←	YfaQ family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,322,124	0	C	.	100.0%	74.3 / NA	18	coding (1153/1650 nt)	yfaQ	YfaQ family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (9/9);  total (9/9)

CCTGCCCCGCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCG  >  NZ_CP009273/2322001‑2322249                                                                                                                            |                                                                                                                              ccTGCCCCGCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGc                                                                                                               >  2:257591/1‑139 (MQ=255)     ccccGCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGa                                                                                                             <  2:219385/137‑1 (MQ=255)                              cAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                                                                                  <  2:118385/139‑1 (MQ=255)                                                        aCTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCtt                                                        >  2:554579/1‑139 (MQ=255)                                                              aTGATAATGAACCGGACCGCCTGCATAGAGGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGAt                                                  <  2:454948/139‑1 (MQ=255)                                                                        aCCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTg                                        <  1:650269/139‑1 (MQ=255)                                                                           ggATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTg                                     >  2:573264/1‑139 (MQ=255)                                                                                          aTGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAAt                      >  1:316805/1‑139 (MQ=255)                                                                                          aTGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAAt                      >  1:346613/1‑139 (MQ=255)                                                                                          aTGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAAt                      >  2:274337/1‑139 (MQ=255)                                                                                               aTCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCa                 <  1:255570/139‑1 (MQ=255)                                                                                               aTCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCa                 <  2:316805/139‑1 (MQ=255)                                                                                               aTCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCa                 >  2:15960/1‑139 (MQ=255)                                                                                               aTCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCa                 >  1:102450/1‑139 (MQ=255)                                                                                                            gCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCttt    >  2:417983/1‑139 (MQ=255)                                                                                                              ggTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCg  <  1:573264/139‑1 (MQ=255)                                                                                                              ggTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCg  <  1:554579/139‑1 (MQ=255)                                                                                                              ggTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCg  <  2:608346/139‑1 (MQ=255)                                                                                                                            |                                                                                                                              CCTGCCCCGCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCG  >  NZ_CP009273/2322001‑2322249

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