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%	237.3 / NA	58	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 . (29/29);  total (29/29)

TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTA  >  NC_000913/3012419‑3012684                                                                                                                                      |                                                                                                                                     tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGTCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                                >  2:912200/1‑139 (MQ=255) tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                                >  1:156882/1‑139 (MQ=255) tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                                >  1:170243/1‑139 (MQ=255)   cGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTg                                                                                                                              >  1:86821/1‑139 (MQ=255)    ggCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTgg                                                                                                                             <  2:86821/139‑1 (MQ=255)         aaGATTGTCAAGTAAGTGCGCTTTGTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGttt                                                                                                                        >  2:677646/1‑139 (MQ=255)          aGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTg                                                                                                                       <  2:170243/139‑1 (MQ=255)              tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                                   <  2:515183/139‑1 (MQ=255)               gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcatc                                                                                                                  <  1:677646/139‑1 (MQ=255)                tCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcatca                                                                                                                 >  2:357837/1‑139 (MQ=255)                          gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGTAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGcc                                                                                                       >  1:1033252/1‑139 (MQ=255)                          gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGTCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGcc                                                                                                       <  1:912200/139‑1 (MQ=255)                            gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  1:1035324/139‑1 (MQ=255)                            gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  1:357837/139‑1 (MQ=255)                                           gATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCggg                                                                                                    <  1:187768/125‑1 (MQ=255)                                           gATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCggg                                                                                                    >  2:187768/1‑125 (MQ=255)                                                                   tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  2:299617/1‑100 (MQ=255)                                                                   tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  2:30815/100‑1 (MQ=255)                                                                   tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  1:30815/1‑100 (MQ=255)                                                                   tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  1:299617/100‑1 (MQ=255)                                                                             cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGc                                                                   >  2:841780/1‑124 (MQ=255)                                                                             cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGc                                                                   <  1:841780/124‑1 (MQ=255)                                                                             cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActtc                                                    >  2:597752/1‑139 (MQ=255)                                                                              aTCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  1:1084108/89‑1 (MQ=255)                                                                              aTCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  2:1084108/1‑89 (MQ=255)                                                                                cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  2:525663/87‑1 (MQ=255)                                                                                cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  1:525663/1‑87 (MQ=255)                                                                                cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCttt                                                 >  1:785237/1‑139 (MQ=255)                                                                                     aTTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCg                                            >  2:457252/1‑139 (MQ=255)                                                                                             gATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                                    >  2:606840/1‑139 (MQ=255)                                                                                             gATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                                    <  2:785237/139‑1 (MQ=255)                                                                                                tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  2:1016968/1‑71 (MQ=255)                                                                                                tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  1:1016968/71‑1 (MQ=255)                                                                                                       cACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCa                          <  1:606840/139‑1 (MQ=255)                                                                                                          ttGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  2:716576/1‑61 (MQ=255)                                                                                                          ttGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  1:716576/61‑1 (MQ=255)                                                                                                            gTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  1:472599/59‑1 (MQ=255)                                                                                                            gTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  2:472599/1‑59 (MQ=255)                                                                                                            gTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAc                     <  2:1033252/139‑1 (MQ=255)                                                                                                                  cTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCAtt               <  1:508263/139‑1 (MQ=255)                                                                                                                    gATATGCGTAGTGCTTCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTTGGTAATACCATCCGGCCATTTTCATTTAAGTCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCAGCAACATCATTaa             <  2:664633/139‑1 (MQ=255)                                                                                                                    gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  2:133048/1‑51 (MQ=255)                                                                                                                    gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  2:848081/1‑51 (MQ=255)                                                                                                                    gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  1:133048/51‑1 (MQ=255)                                                                                                                    gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  1:848081/51‑1 (MQ=255)                                                                                                                    gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa             <  1:124194/139‑1 (MQ=255)                                                                                                                    gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa             >  1:1102422/1‑139 (MQ=255)                                                                                                                     aTAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggt                                                                                                   <  1:697899/52‑1 (MQ=255)                                                                                                                     aTAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggt                                                                                                   >  2:697899/1‑52 (MQ=255)                                                                                                                          cGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                                                                        <  1:99997/42‑1 (MQ=18)                                                                                                                          cGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                                                                        >  2:99997/1‑42 (MQ=18)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                                <  1:148964/45‑1 (MQ=25)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                                >  2:148964/1‑45 (MQ=25)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     >  1:950343/1‑40 (MQ=21)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                     <  2:950343/40‑1 (MQ=25)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAca                    >  2:1046545/1‑121 (MQ=255)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAca                    <  1:1046545/121‑1 (MQ=255)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTa  >  2:50045/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                     TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTA  >  NC_000913/3012419‑3012684

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