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%	196.9 / NA	40	intergenic (+89/‑6)	yggI/endA	conserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (17/23);  total (17/23)

ACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGA  >  NC_000913/3090296‑3090383                                              |                                           aCCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTg                                       >  1:262975/1‑50 (MQ=255)  ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGc                                      <  1:1320325‑M1/50‑2 (MQ=255)   cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcc                                     >  1:415964‑M1/1‑48 (MQ=255)    aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcct                                    <  1:1341308‑M1/50‑4 (MQ=255)    aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcc                                     >  1:2679679‑M1/1‑47 (MQ=255)     aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGccta                                   <  1:3077712‑M1/50‑5 (MQ=255)      aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctac                                  <  1:3394662‑M1/50‑6 (MQ=255)       aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctaca                                 <  1:642816‑M1/50‑7 (MQ=255)         aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacact                               <  1:1051296‑M1/50‑9 (MQ=255)         aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacact                               <  1:744965‑M1/50‑9 (MQ=255)          gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGCGcctacac                                >  1:1132564‑M1/1‑41 (MQ=255)           cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactag                             >  1:546636‑M1/1‑40 (MQ=255)           cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactag                             <  1:2250199‑M1/50‑11 (MQ=255)           cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactag                             >  1:3387129‑M1/1‑40 (MQ=255)              tCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcg                           <  1:2273603‑M1/49‑13 (MQ=255)                gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                        >  1:1258621‑M1/1‑35 (MQ=255)                gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                        >  1:95379‑M1/1‑35 (MQ=255)                gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggg                         <  1:3050630‑M1/49‑15 (MQ=255)                    cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggatt                      <  1:2668759‑M1/48‑18 (MQ=255)                     gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                   <  1:271934‑M1/50‑21 (MQ=255)                     gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                   <  1:3822056‑M1/50‑21 (MQ=255)                     gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                   <  1:3540507‑M1/50‑21 (MQ=255)                      ttGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttt                  >  1:3725528‑M1/1‑29 (MQ=255)                        gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttt                 <  1:1285663‑M1/49‑23 (MQ=255)                        gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttt                 <  1:1820820‑M1/49‑23 (MQ=255)                         cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttg               >  1:1166289‑M1/1‑26 (MQ=255)                         cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttg               <  1:406124‑M1/50‑25 (MQ=255)                         cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttg               >  1:4119998‑M1/1‑26 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt             >  1:974735‑M1/1‑23 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt             >  1:312671‑M1/1‑23 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt             <  1:1768922‑M1/49‑27 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt              <  1:356587‑M1/48‑26 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt              <  1:45828‑M1/48‑26 (MQ=255)                             cTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaa           >  1:1010968‑M1/1‑22 (MQ=255)                                gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacc         <  1:2561045‑M1/49‑31 (MQ=255)                                gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaac          <  1:1482643‑M1/48‑30 (MQ=255)                                  tttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccc       >  1:4062636‑M1/1‑17 (MQ=255)                                    tAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccccta     <  1:3294436‑M1/49‑35 (MQ=255)                                     aaCACGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacc   >  1:1873966‑M1/1‑14 (MQ=255)                                      acacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctaccc  >  1:2639066‑M1/1‑13 (MQ=255)                                              |                                           ACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGA  >  NC_000913/3090296‑3090383

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