Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,166,540	T→C	R291R (CGT→CGC)	yegQ →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,166,540	0	T	C	100.0%	181.1 / NA	48	R291R (CGT→CGC)	yegQ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/23);  total (25/23)

CGTCCGGGCGAGTATATGACCGCGTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGTACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGAAACTCTGGAAGGTCTGGCGCATCGTGGCTATAC  >  NC_000913/2166406‑2166671                                                                                                                                       |                                                                                                                                    cGTCCGGGCGAGTATATGACCGCGTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCa                                                                                                                                 >  2:373239/1‑139 (MQ=255) cGTCCGGGCGAGTATATGACCGCGTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCa                                                                                                                                 <  2:450796/139‑1 (MQ=255)     cGGGCGAGTATATGACCGCGTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATc                                                                                                                             >  1:901367/1‑139 (MQ=255)      gggCGAGTATATGACCGCGTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCt                                                                                                                            <  2:231430/139‑1 (MQ=255)           aGTATATGACCGCGTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTcta                                                                                                                       <  1:107680/139‑1 (MQ=255)               tatGACCGCGTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTAt                                                                                                                   <  1:4788/139‑1 (MQ=255)                  gACCGCGTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATtgt                                                                                                                >  1:1054107/1‑139 (MQ=255)                     cgcgTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTgca                                                                                                             >  2:429429/1‑139 (MQ=255)                       cgTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTgcacg                                                                                                           >  2:459687/1‑139 (MQ=255)                       cgTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTgcacg                                                                                                           >  2:296572/1‑139 (MQ=255)                         tttGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTgcacgca                                                                                                         <  2:364589/139‑1 (MQ=255)                          ttGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTgcacgcac                                                                                                        <  2:748498/139‑1 (MQ=255)                                  gAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAgg                                                                                                >  1:916031/1‑139 (MQ=255)                                       cGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTAc                                                                                           >  2:575579/1‑139 (MQ=255)                                        ggCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTAcc                                                                                          >  1:164004/1‑139 (MQ=255)                                           aCTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCa                                                                                       <  2:280980/139‑1 (MQ=255)                                              tACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAg                                                                                    <  1:373239/139‑1 (MQ=255)                                                    atGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCg                                                                              >  1:632448/1‑139 (MQ=255)                                                          tCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGAcg                                                                        >  1:922908/1‑139 (MQ=255)                                                          tCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGAcg                                                                        >  1:381128/1‑139 (MQ=255)                                                                           cATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCg                                                       >  2:459761/1‑139 (MQ=255)                                                                           cATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCg                                                       <  2:128721/139‑1 (MQ=255)                                                                             tCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGtt                                                     <  2:922908/139‑1 (MQ=255)                                                                             tCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGtt                                                     <  2:381128/139‑1 (MQ=255)                                                                                     gTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCa                                             <  2:164004/139‑1 (MQ=255)                                                                                       aGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGc                                           <  1:459687/139‑1 (MQ=255)                                                                                           cGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCctgc                                       >  1:1104/1‑139 (MQ=255)                                                                                           cGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCctgc                                       >  1:581300/1‑139 (MQ=255)                                                                                           cGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCctgc                                       >  1:139818/1‑139 (MQ=255)                                                                                           cGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCctgc                                       >  1:382810/1‑139 (MQ=255)                                                                                             ccTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAAt                                                                                                                              <  1:446125/50‑1 (MQ=255)                                                                                             ccTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAAt                                                                                                                              >  2:446125/1‑50 (MQ=255)                                                                                              cTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTgg                                    <  1:459761/139‑1 (MQ=255)                                                                                              cTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTgg                                    >  2:202422/1‑139 (MQ=255)                                                                                               tGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGa                                   >  2:933207/1‑139 (MQ=255)                                                                                                gACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGaa                                  >  2:725450/1‑139 (MQ=255)                                                                                                       aTGGGCGTGCATTAACTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCGCCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGAAACTCTgg                           <  1:351411/139‑1 (MQ=255)                                                                                                        tGGGCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGAAACTCTGGa                          <  1:575579/139‑1 (MQ=255)                                                                                                          ggCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGAAACTCTGGAAg                        >  1:35222/1‑139 (MQ=255)                                                                                                           gCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGAcgccgc                                                                    <  2:714581/94‑1 (MQ=255)                                                                                                           gCGTGCATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGAcgccgc                                                                    >  1:714581/1‑94 (MQ=255)                                                                                                                cATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGAAACTCTGGAAGGTCTgg                  >  1:322539/1‑139 (MQ=255)                                                                                                                cATTCGCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGAAACTCTGGAAGGTCTgg                  <  2:632448/139‑1 (MQ=255)                                                                                                                     gCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTgca                                                                                                             <  1:781743/43‑1 (MQ=255)                                                                                                                     gCTGAAAATCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGAAACTCTGGAAGGTCTGGCGCAt             <  1:725450/139‑1 (MQ=255)                                                                                                                             tCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGtt                                                     >  2:298085/1‑91 (MQ=255)                                                                                                                             tCGAAGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGtt                                                     <  1:298085/91‑1 (MQ=255)                                                                                                                                aaGGTCGCACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGAAACTCTGGAAGGTCTGGCGCATCGTGGCTATAc  <  2:1104/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    CGTCCGGGCGAGTATATGACCGCGTTTGAAGATGAGCACGGCACTTACATCATGAACTCGAAAGATCTGCGCGCCATCGCCCATGTAGAACGCCTGACCAAAATGGGCGTGCATTCGCTGAAAATCGAAGGTCGTACCAAATCTTTCTACTATTGTGCACGCACCGCACAGGTTTACCGCAAAGCTATCGATGACGCCGCTGCGGGAAAACCGTTCGATACCAGCCTGCTGGAAACTCTGGAAGGTCTGGCGCATCGTGGCTATAC  >  NC_000913/2166406‑2166671

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