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%	112.7 / NA	32	D172H (GAT→CAT)	yajG	putative lipoprotein YajG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (17/15);  total (17/15)

AGGGCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATCGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAGCTGGC  >  NC_000913/453579‑453729                                                                            |                                                                            agGGCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCa                                                                        <  1:87814/81‑1 (MQ=255)    gCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGca                                                                     <  2:87815/81‑1 (MQ=255)    gCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGca                                                                     >  1:87861/1‑81 (MQ=255)     cAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGcac                                                                    <  1:87816/81‑1 (MQ=255)      aGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGcaca                                                                   >  2:87857/1‑81 (MQ=255)       gAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCacac                                                                  >  2:87866/1‑81 (MQ=255)          aTTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGt                                                               >  1:87845/1‑81 (MQ=255)          aTTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGt                                                               >  1:87875/1‑81 (MQ=255)                  gcgTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCa                                                       >  2:87853/1‑81 (MQ=255)                         gCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGa                                                >  1:87854/1‑81 (MQ=255)                          cTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGc                                                  <  1:87818/78‑1 (MQ=255)                          cTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGa                                                <  1:87817/80‑1 (MQ=255)                             gatgaATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGatat                                             <  2:87819/80‑1 (MQ=255)                                gaATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATAttttt                                         >  2:87881/1‑81 (MQ=255)                                gaATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATAttttt                                         >  1:87859/1‑81 (MQ=255)                                   ttCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTttgtt                                      <  2:87820/81‑1 (MQ=255)                                      gTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTttgttg                                     >  2:87877/1‑79 (MQ=255)                                       tGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAg                                  <  1:87821/81‑1 (MQ=255)                                          aTGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGcc                               >  1:87871/1‑81 (MQ=255)                                          aTGCTGGTGTCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGcc                               >  1:87863/1‑81 (MQ=255)                                                gtgtCCTGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAAc                         <  1:87822/81‑1 (MQ=255)                                                      tGAGACATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCt                   >  2:87872/1‑81 (MQ=255)                                                        agaCATATCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCt                   <  1:87823/79‑1 (MQ=255)                                                               tCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGtt           <  2:87824/80‑1 (MQ=255)                                                               tCAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGtt           <  1:87825/80‑1 (MQ=255)                                                                cAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTg          >  2:87869/1‑80 (MQ=255)                                                                cAGCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTg          >  2:87870/1‑80 (MQ=255)                                                                  gCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAg       >  1:87873/1‑81 (MQ=255)                                                                  gCGATGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAg       >  1:87865/1‑81 (MQ=255)                                                                     aTGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAGCTg    <  2:87826/81‑1 (MQ=255)                                                                      tGGTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAGCTgg   <  2:87827/81‑1 (MQ=255)                                                                        gTATGGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAGCTGGc  <  1:87828/80‑1 (MQ=255)                                                                            |                                                                            AGGGCAGAAATTAACGCGCGTTCTGCTTGATGAATTCGTGGATGCTGGTGTCCTGAGACATATCAGCGATGGTATCGGTCAGCACACTGTTAACCGCATCGGCGATATTTTTGTTGGAGGCCTGGAACGCACCTTCAACGTTGTAGCTGGC  >  NC_000913/453579‑453729

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: ‑