Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,674,542	C→A	P49P (CCG→CCT)	yphB ←	putative aldose 1‑epimerase YphB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,674,542	0	C	A	100.0%	152.0 / NA	42	P49P (CCG→CCT)	yphB	putative aldose 1‑epimerase YphB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/22);  total (20/22)

AACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGCGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATCGCGCCAGAAGCC  >  NC_000913/2674464‑2674619                                                                               |                                                                              aaCATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAgg                                                                             <  1:476445/81‑1 (MQ=255)       cGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAAc                                                                        <  2:476446/80‑1 (MQ=255)        ggTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACAt                                                                      <  1:476447/81‑1 (MQ=255)         gTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATg                                                                     <  1:476448/81‑1 (MQ=255)            gCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAgg                                                                  >  2:476495/1‑81 (MQ=255)             cAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAgg                                                                  >  2:476499/1‑80 (MQ=255)                  gATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCg                                                             <  2:476449/80‑1 (MQ=255)                    tACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGtt                                                          <  1:476452/81‑1 (MQ=255)                    tACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGtt                                                          <  1:476451/81‑1 (MQ=255)                    tACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGtt                                                          >  1:476510/1‑81 (MQ=255)                    tACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGtt                                                          <  2:476450/81‑1 (MQ=255)                       tCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGc                                                        >  2:476517/1‑80 (MQ=255)                        cACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCca                                                      >  1:476501/1‑81 (MQ=255)                         aCGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGTGGCATCGGTTGCca                                                      >  1:476507/1‑80 (MQ=255)                           gCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGTAAAACATGAGGCATCGGTTGCCacac                                                   <  2:476453/81‑1 (MQ=255)                            ccccTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACAcc                                                  <  2:476454/81‑1 (MQ=255)                                 gCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCttt                                             >  2:476503/1‑81 (MQ=255)                                  ccAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCtttt                                            >  1:476504/1‑81 (MQ=255)                                  ccAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCtttt                                            <  1:476455/81‑1 (MQ=255)                                  ccAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCttt                                             <  1:476456/80‑1 (MQ=255)                                        aaaaCGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTAc                                       >  1:476514/1‑80 (MQ=255)                                         aaaCGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTAc                                       >  2:476513/1‑79 (MQ=255)                                           aCGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAgg                                   <  1:476457/81‑1 (MQ=255)                                             gATTGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGAc                                 >  2:476524/1‑81 (MQ=255)                                               ttGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGt                               <  2:476458/81‑1 (MQ=255)                                               ttGCCGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGt                               <  2:476459/81‑1 (MQ=255)                                                   cGCTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTaa                             >  2:476528/1‑79 (MQ=255)                                                     cTCACCCGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAgc                         >  2:476506/1‑81 (MQ=255)                                                           cGGTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGcgtc                   <  1:476460/81‑1 (MQ=255)                                                             gTTTGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGcgtcgt                 <  1:476461/81‑1 (MQ=255)                                                              tttGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTa                >  1:476519/1‑81 (MQ=255)                                                               ttGCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTAt               >  2:476530/1‑81 (MQ=255)                                                                 gCAAATGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATcg             >  1:476511/1‑81 (MQ=255)                                                                  cAAATGGCACCAGAGGAAAACATTAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATcg             <  2:476462/80‑1 (MQ=255)                                                                   aaaTGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATcgcg           <  1:476464/81‑1 (MQ=255)                                                                   aaaTGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATcgcg           <  1:476463/81‑1 (MQ=255)                                                                    aaTGGCACCAGAGGAAAACATGAGGCATCGGTTGCCAGACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATCgcgc          <  2:476465/81‑1 (MQ=255)                                                                    aaTGGCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATCgcgc          >  2:476531/1‑81 (MQ=255)                                                                       ggCACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATCgcgc          >  2:476508/1‑78 (MQ=255)                                                                         cACCAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATCGCGCCAg       >  2:476529/1‑79 (MQ=255)                                                                           ccAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATCGCGCCAGAAGc   <  2:476466/81‑1 (MQ=255)                                                                            cAGAGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATCGCGCCAGAAGcc  >  1:476542/1‑81 (MQ=255)                                                                               |                                                                              AACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGCAAATGGCACCAGCGGAAAACATGAGGCATCGGTTGCCACACCGCTTTTTTTACCAGGACGTAAGAGCGGCGTCGTATCGCGCCAGAAGCC  >  NC_000913/2674464‑2674619

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