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%	96.0 / NA	23	coding (1153/1650 nt)	yfaQ	YfaQ family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (13/10);  total (13/10)

TAACGCTCTCCCTGCCCCGCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCG  >  NZ_CP009273/2321991‑2322249                                                                                                                                      |                                                                                                                              tAACGCTCTCCCTGCCCCGCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGccccgt                                                                                                                          >  1:19385/1‑137 (MQ=255)   aCGCTCTCCCTGCCCCGCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGccccgtgg                                                                                                                        >  2:501727/1‑135 (MQ=255)                 ccGCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCa                                                                                                         <  1:501727/139‑1 (MQ=255)                   gCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAg                                                                                                        >  1:229943/1‑138 (MQ=255)                   gCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAg                                                                                                        >  2:346449/1‑138 (MQ=255)                        ggCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAacac                                                                                                  >  2:440226/1‑139 (MQ=255)                        ggCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAaca                                                                                                   <  2:347661/138‑1 (MQ=255)                         gCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACg                                                                                                 <  1:277713/139‑1 (MQ=255)                            gTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAACAGAAACACGCTg                                                                                              >  2:714/1‑139 (MQ=255)                                   gccgccAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTgg                                                                                         <  1:390491/137‑1 (MQ=255)                                    ccgccAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGgcgc                                                                                      >  2:78990/1‑139 (MQ=255)                                      gccAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctg                                                                                    <  1:536489/139‑1 (MQ=255)                                                  tCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATg                                                                        >  2:99676/1‑139 (MQ=255)                                                                                     ggATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTg                                     >  1:613088/1‑139 (MQ=255)                                                                                      gATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCg                                                                                                              <  2:392240/65‑1 (MQ=255)                                                                                      gATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCg                                                                                                              >  1:392240/1‑65 (MQ=255)                                                                                      gATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTg                                     >  1:530897/1‑138 (MQ=255)                                                                                               cATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTg                           >  1:620343/1‑139 (MQ=255)                                                                                                   gATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACGCGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGACTTCCCGATTGGCGGTTTGCTGTAAAAAAGTGCGa                        <  1:112736/138‑1 (MQ=255)                                                                                                         aTCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCa                 >  1:348584/1‑139 (MQ=255)                                                                                                         aTCCTGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCa                 <  1:346449/139‑1 (MQ=255)                                                                                                             tGCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGt             >  2:484703/1‑139 (MQ=255)                                                                                                              gCGTCCAGGCGGTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGt             >  2:187640/1‑138 (MQ=255)                                                                                                                        ggTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCg  <  1:78990/139‑1 (MQ=255)                                                                                                                        ggTCATCGTTCGCG‑CCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCg  <  2:613088/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                              TAACGCTCTCCCTGCCCCGCCTGGGCCGTGGCTTGCCGCCAGGAAAGTGTCCCTTCGGTAGCGCGACTGCCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCG  >  NZ_CP009273/2321991‑2322249

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