Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,757,387	C→G	intergenic (+277/‑34)	pykF → / → lpp	pyruvate kinase I/murein lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,757,387	0	C	G	100.0%	89.2 / NA	24	intergenic (+277/‑34)	pykF/lpp	pyruvate kinase I/murein lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (5/19);  total (5/19)

GCCTTCCCATCAAAAAAATATTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTACATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGCTACTAAACTGGTACTGGGCGCG  >  NC_000913/1757327‑1757450                                                             |                                                                gCCTTCCCATCAAAAAAATATTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGgaga                                                          <  1:220503/68‑1 (MQ=255)    ttCCCATCAAAAAAATATTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTa                                                       <  1:928096/68‑1 (MQ=255)      cccATCAAAAAAATATTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTAAc                                                     <  2:86757/68‑1 (MQ=255)          tCAAAAAAATATTCTCAACATAAAAAACTTTGTGTAATACTTGTATCGCTAGATGTAGATTAACTCat                                                 <  2:136855/68‑2 (MQ=255)           cAATAAAATATTCTCTACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAAt                                                <  1:53513/68‑1 (MQ=255)           cAAAAAAATATTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAAt                                                <  1:190667/68‑1 (MQ=255)                  atatTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAgg                                         <  2:529685/68‑1 (MQ=255)                  atatTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAgg                                         >  1:472997/1‑68 (MQ=255)                   aatTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAggg                                        <  1:388549/67‑1 (MQ=255)                    atTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGt                                       <  1:137870/68‑1 (MQ=255)                      tctcAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTAAc                                                     <  2:811734/52‑1 (MQ=255)                      tctcAACATAAAAAACTTTGTGTAATACTTGTAACGCTAGATGGAGATTAAc                                                     >  1:811737/1‑52 (MQ=255)                            cATAAAAAACTATGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGCATtaataa                               <  2:919685/68‑1 (MQ=255)                                    aCTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGct                       <  1:964737/68‑1 (MQ=255)                                    aCTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGct                       <  2:299343/68‑1 (MQ=255)                                    aCTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGct                       <  1:278635/68‑1 (MQ=255)                                    aCTTTGTGTAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGct                       <  2:828762/68‑1 (MQ=255)                                          tgtAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGCTACTaat                 <  2:916522/68‑2 (MQ=255)                                          tgtAATACTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGCTACTaaa                 <  1:277104/68‑1 (MQ=255)                                                aCTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGTGTATTAATAATGAAAGCTACTaaa                 <  1:878356/62‑1 (MQ=255)                                                aCTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGCTACTaaa                 >  2:878353/1‑62 (MQ=255)                                                aCTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGCTACTAAACTGGTa           >  1:29023/1‑68 (MQ=255)                                                aCTTGTAACGCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGCTACTAAACTGGTa           >  1:21458/1‑68 (MQ=255)                                                         gCTAGATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGCTACTAAACTGGTACTGGGcgcg  <  1:658336/68‑1 (MQ=255)                                                             |                                                                GCCTTCCCATCAAAAAAATATTCTCAACATAAAAAACTTTGTGTAATACTTGTAACGCTACATGGAGATTAACTCAATCTAGAGGGTATTAATAATGAAAGCTACTAAACTGGTACTGGGCGCG  >  NC_000913/1757327‑1757450

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