Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,012,553	Δ1 bp	intergenic (+160/+61)	yqeA → / ← yqeB	putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,012,552	0	G	.	100.0%	132.0 / NA	32	intergenic (+159/+62)	yqeA/yqeB	putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (16/16);  total (16/16)

CCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAA  >  NC_000913/3012416‑3012673                                                                                                                                         |                                                                                                                          ccATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAgca                                                                                                                         >  1:727451/1‑137 (MQ=255) ccATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAgca                                                                                                                         >  1:614160/1‑137 (MQ=255)    tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                     >  1:272088/1‑139 (MQ=255)    tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                     >  1:401195/1‑139 (MQ=255)    tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                     >  1:650309/1‑139 (MQ=255)    tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                     >  1:49838/1‑139 (MQ=255)       ggCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTgg                                                                                                                  <  1:58957/139‑1 (MQ=255)                 tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                        <  2:650309/139‑1 (MQ=255)                 tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                        <  2:286907/139‑1 (MQ=255)                            tGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                                                             >  1:717640/1‑139 (MQ=255)                               gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  2:614160/139‑1 (MQ=255)                               gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  2:324123/139‑1 (MQ=255)                               gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  1:381202/139‑1 (MQ=255)                                 tttGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggt                                                                                        >  1:131660/1‑139 (MQ=255)                                                          gcgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCa                                                               >  1:716441/1‑139 (MQ=255)                                                           cgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt                                                              >  1:11857/1‑139 (MQ=255)                                                               ttATCCGTCCTACAAAACATCGCAAATTCAAAAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTTGGTAATACCATCCGGCCATTTTCATTTaa                                                          <  1:310949/139‑1 (MQ=255)                                                                       ccTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          >  1:30226/1‑99 (MQ=255)                                                                       ccTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  2:30226/99‑1 (MQ=255)                                                                                        aTTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCg                                 <  2:716441/139‑1 (MQ=255)                                                                                        aTTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGACGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCg                                 >  1:489929/1‑139 (MQ=255)                                                                                              tAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                           <  2:725505/139‑1 (MQ=255)                                                                                                   tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                           <  1:434862/134‑1 (MQ=255)                                                                                                   tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                           >  2:434862/1‑134 (MQ=255)                                                                                                      aaTGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGa                   <  2:49838/139‑1 (MQ=255)                                                                                                           aCTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCAt              <  2:131660/139‑1 (MQ=255)                                                                                                              tGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCaa           <  2:555097/139‑1 (MQ=255)                                                                                                               gTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAc          >  2:209806/1‑139 (MQ=255)                                                                                                                  ggCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAcatc       >  1:624597/1‑139 (MQ=255)                                                                                                                    ccTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTc                                  <  2:187818/110‑1 (MQ=255)                                                                                                                    ccTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTc                                  >  1:187818/1‑110 (MQ=255)                                                                                                                       gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTATCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa  >  1:475687/1‑139 (MQ=255)                                                                                                                                  cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa  >  1:79671/1‑128 (MQ=255)                                                                                                                                  cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa  <  2:79671/128‑1 (MQ=255)                                                                                                                                         |                                                                                                                          CCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAA  >  NC_000913/3012416‑3012673

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