Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,090,341	Δ1 bp	intergenic (+89/‑6)	yggI → / → endA	conserved protein YggI/DNA‑specific endonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,090,341	0	T	.	100.0%	149.4 / NA	31	intergenic (+89/‑6)	yggI/endA	conserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (13/18);  total (13/18)

CAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAGC  >  NC_000913/3090298‑3090385                                            |                                             cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcc                                       <  1:1183629‑M1/50‑3 (MQ=255)    aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGccta                                     >  1:1913277‑M1/1‑45 (MQ=255)      cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                  >  1:1865159‑M1/1‑43 (MQ=255)       aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                  >  1:1336511‑M1/1‑42 (MQ=255)           ttCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcg                             >  1:1124583‑M1/1‑38 (MQ=255)           ttCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcg                             >  1:202875‑M1/1‑38 (MQ=255)              gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                          >  1:1511848‑M1/1‑35 (MQ=255)               cTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggat                         <  1:20500‑M1/50‑17 (MQ=255)                tACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggatt                        >  1:1260981‑M1/1‑33 (MQ=255)                 aCGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattc                       >  1:1480530‑M1/1‑32 (MQ=255)                   gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                     <  1:1042813‑M1/50‑21 (MQ=255)                    ttGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttt                    >  1:1544527‑M1/1‑29 (MQ=255)                      gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttttt                  >  1:679274‑M1/1‑27 (MQ=255)                      gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttt                   <  1:974142‑M1/49‑23 (MQ=255)                       cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttttt                  <  1:1992066‑M1/49‑24 (MQ=255)                       cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttttt                  <  1:69237‑M1/49‑24 (MQ=255)                          gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta              <  1:1588398‑M1/50‑28 (MQ=255)                          gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta              <  1:641802‑M1/50‑28 (MQ=255)                          gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta              <  1:2109256‑M1/50‑28 (MQ=255)                          gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt               <  1:385195‑M1/49‑27 (MQ=255)                               ttttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccc         >  1:1076157‑M1/1‑18 (MQ=255)                                tttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccc         >  1:1787023‑M1/1‑17 (MQ=255)                                 ttAACACGGAG‑AAGTGcctacgctagcgggattctttttgttaacccct        <  1:2119301‑M1/49‑34 (MQ=255)                                 ttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccct        <  1:2208914‑M1/49‑34 (MQ=255)                                  tAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccccta       <  1:2031307‑M1/49‑35 (MQ=255)                                   aaCACGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacc     >  1:408566‑M1/1‑14 (MQ=255)                                    acacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacc     <  1:1325966‑M1/49‑37 (MQ=255)                                    acacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacc     <  1:1328159‑M1/49‑37 (MQ=255)                                     cacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacccc   <  1:2202855‑M1/50‑39 (MQ=255)                                     cacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacccc   <  1:1275652‑M1/50‑39 (MQ=255)                                      acGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacccca  <  1:24276‑M1/50‑40 (MQ=255)                                            |                                             CAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAGC  >  NC_000913/3090298‑3090385

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