Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,116,466	C→T	G883S (GGC→AGC)	yghJ ←	putative lipoprotein YghJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,116,466	0	C	T	100.0%	93.5 / NA	26	G883S (GGC→AGC)	yghJ	putative lipoprotein YghJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/14);  total (12/14)

GTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAA  >  NC_000913/3116398‑3116535                                                                     |                                                                      gtCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTc                                                           >  2:571624/1‑81 (MQ=255)  tCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTcc                                                          <  1:571597/81‑1 (MQ=255)    ggCGTTAAGGCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCa                                                         <  2:571598/80‑1 (MQ=255)    ggCGTTAAGGCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAg                                                        >  2:571626/1‑81 (MQ=255)         tAAGGCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGt                                                   >  1:571631/1‑81 (MQ=255)          aaGGCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGtt                                                  >  2:571630/1‑81 (MQ=255)          aaGGCTTAGTTGCGTATACTGTGGAACATACATGCAACCAGTAACCGGAACACCCGTGCTAGGACGATATTACAGGCAGtt                                                  <  2:951337/81‑1 (MQ=255)           aGGCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTg                                                 <  2:571599/81‑1 (MQ=255)             gCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGAc                                               <  2:571600/81‑1 (MQ=255)                 aGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTc                                            <  2:571601/80‑1 (MQ=255)                 aGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCa                                           >  1:571628/1‑81 (MQ=255)                  gTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCAt                                          <  2:571602/81‑1 (MQ=255)                   ttGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCa                                           >  2:571627/1‑79 (MQ=255)                      cGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTg                                      <  1:571603/81‑1 (MQ=255)                            cTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGc                                <  2:571604/81‑1 (MQ=255)                              gtgGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCg                              >  1:571632/1‑81 (MQ=255)                               tgGAACATACATGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCg                              >  2:571643/1‑80 (MQ=255)                                     atacatGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAg                        >  1:571636/1‑80 (MQ=255)                                     atacatGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGt                       <  1:571605/81‑1 (MQ=255)                                      tacatGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTa                      >  1:571660/1‑81 (MQ=255)                                        catGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTaca                    >  1:571637/1‑81 (MQ=255)                                        catGCCACCAGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTac                     <  1:571606/80‑1 (MQ=255)                                              accaGTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTACACTCttt              <  2:571607/81‑1 (MQ=255)                                                  gTAACCGGAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTACACTCTTTCAGc          >  1:571633/1‑81 (MQ=255)                                                        ggAACCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTACACTCTTTCAGCCCCggg    <  1:571608/81‑1 (MQ=255)                                                          aaCCCCCGTGCTAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTACACTCTTTCACCCCCGGGaa  <  2:571609/81‑1 (MQ=255)                                                                     |                                                                      GTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACATACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAA  >  NC_000913/3116398‑3116535

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