Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	453,654	C→G	D172H (GAT→CAT)	yajG ←	putative lipoprotein YajG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	453,654	0	C	G	100.0%	118.7 / NA	33	D172H (GAT→CAT)	yajG	putative lipoprotein YajG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (17/16);  total (17/16)

AGAGGGCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATCGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAGCT  >  NC_000913/453577‑453726                                                                              |                                                                         agagGGCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGt                                                                       <  2:96469/81‑1 (MQ=255)  gagGGCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTc                                                                      <  2:96470/81‑1 (MQ=255)   agGGCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGt                                                                       >  1:96490/1‑79 (MQ=255)   agGGCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCa                                                                     <  1:96471/81‑1 (MQ=255)     ggCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAg                                                                    >  2:96510/1‑80 (MQ=255)     ggCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGc                                                                   <  2:96472/81‑1 (MQ=255)      gCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCTTGAGACATATCAGCGATGGTATGGGTCAGca                                                                  >  2:96526/1‑81 (MQ=255)         gAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCacac                                                               <  1:96473/81‑1 (MQ=255)          aaaTTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACt                                                              >  2:96515/1‑81 (MQ=255)              tAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGtt                                                           >  1:96523/1‑80 (MQ=255)               aaCGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTaa                                                         >  2:96513/1‑81 (MQ=255)                   cgcgTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGc                                                     <  1:96474/81‑1 (MQ=255)                            cTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCg                                              <  1:96475/79‑1 (MQ=255)                                atgaATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATAttt                                        <  1:96476/81‑1 (MQ=255)                                  gaATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATAttttt                                      >  2:96539/1‑81 (MQ=255)                                   aaTTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTttg                                     <  1:96477/81‑1 (MQ=255)                                        gTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTgg                                 <  2:96478/80‑1 (MQ=255)                                            aTGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGcc                            >  1:96531/1‑81 (MQ=255)                                            aTGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGcc                            >  1:96516/1‑81 (MQ=255)                                               cTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTg                          <  1:96479/80‑1 (MQ=255)                                                  gtgtCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGa                        <  1:96481/79‑1 (MQ=255)                                                  gtgtCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAAc                      >  1:96508/1‑81 (MQ=255)                                                  gtgtCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAAc                      <  2:96480/81‑1 (MQ=255)                                                    gtCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACg                     >  1:96533/1‑80 (MQ=255)                                                    gtCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACg                     >  1:96521/1‑80 (MQ=255)                                                      ccTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCAc                  >  2:96527/1‑81 (MQ=255)                                                           gaCATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTc              >  1:96482/1‑80 (MQ=255)                                                            aCATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCaa            >  2:96518/1‑81 (MQ=255)                                                             cATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAAc           <  1:96483/81‑1 (MQ=255)                                                             cATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAAc           <  2:194135/81‑1 (MQ=255)                                                                   aGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTg       >  2:96517/1‑79 (MQ=255)                                                                    gCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAg    >  2:96524/1‑81 (MQ=255)                                                                      gATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAGCt  <  2:96484/81‑1 (MQ=255)                                                                              |                                                                         AGAGGGCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATCGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAGCT  >  NC_000913/453577‑453726

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