Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,116,901	C→A	A738S (GCG→TCG)	yghJ ←	putative lipoprotein YghJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,116,901	0	C	A	100.0%	86.5 / NA	28	A738S (GCG→TCG)	yghJ	putative lipoprotein YghJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (14/14);  total (14/14)

AAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAG  >  NC_000913/3116778‑3117027                                                                                                                            |                                                                                                                               aaaTGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACCGACACACAAAACAAGACCAGCTCTCTTc                                                                                                                       >  2:177000/1‑133 (MQ=255)   aTGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCt                                                                                                                     >  2:75962/1‑133 (MQ=255)       ccGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAg                                                                                                                 >  2:285420/1‑133 (MQ=255)        cGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGa                                                                                                                <  2:500715/133‑1 (MQ=255)               cgcTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCtca                                                                                                         >  1:361173/1‑133 (MQ=255)                          aaCGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTcc                                                                                              <  1:269693/133‑1 (MQ=255)                                     gATCCCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCAc                                                                                   <  2:272983/133‑1 (MQ=255)                                     gATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCAc                                                                                   <  1:75962/133‑1 (MQ=255)                                          ccTTGCGGATCGTTATTTACTCCCGACTTGTTCAGTGCCATCGACACACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATc                                                                              <  1:419333/133‑1 (MQ=255)                                           cTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATcc                                                                             >  1:476927/1‑133 (MQ=255)                                             tGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTCAAGATTGCTCATCACGGTTTCCATGGTCAGGACCGATCCAc                                                                           >  2:170722/1‑133 (MQ=255)                                             tGCGGATCGTTATTTACTACCCACTTGTTCAGTGCCATCGACACACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAACATTGCGCATCACGATTTCCATGCTCAGAACCGATCCAc                                                                           >  2:193105/1‑133 (MQ=255)                                              gCGGATCGTTGTTTACTACCGACTCGTTGAGTTCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCAcc                                                                          <  2:172624/133‑1 (MQ=255)                                               cGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCt                                                                         >  1:19551/1‑133 (MQ=255)                                                 gATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCttt                                                                       <  2:476927/133‑1 (MQ=255)                                                    cGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGtt                                                                    >  2:490308/1‑133 (MQ=255)                                                      ttATTTACGACCGACGTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCa                                                                  <  1:231868/133‑1 (MQ=255)                                                        aTTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGa                                                                <  1:490308/133‑1 (MQ=255)                                                           tactacCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCa                                                                  <  2:20053/128‑1 (MQ=255)                                                           tactacCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCa                                                                  >  1:20053/1‑128 (MQ=255)                                                           tactacCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAg                                                             <  1:60210/133‑1 (MQ=255)                                                               acCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGa                                                         >  1:302580/1‑133 (MQ=255)                                                                             tGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTcaca                                           >  1:263612/1‑133 (MQ=255)                                                                             tGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTACTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTcaca                                           >  2:112478/1‑133 (MQ=255)                                                                                  tCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTGCAGATAGGCGAACACATCGGTTACATCctg                                      >  1:199794/1‑133 (MQ=255)                                                                                       aGACCTGCGGCATCCAACAGACGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAg                                 <  1:285420/133‑1 (MQ=255)                                                                                                           aCGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTa             >  1:473981/1‑133 (MQ=255)                                                                                                            cGCACAAAACCAGACGAGCTCTCTTCCTTAAGATTACTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTAt            <  1:112478/133‑1 (MQ=255)                                                                                                                      cAGACGAGCTCTCTTGCTTAAGTTTGCTCAGCACGTTTTCCATGATCAGCACCGATCCCCCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAg  <  1:193105/133‑1 (MQ=255)                                                                                                                            |                                                                                                                               AAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAG  >  NC_000913/3116778‑3117027

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