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

ACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAG  >  NC_000913/3090296‑3090384                                              |                                            aCCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTg                                        >  1:4056649/1‑50 (MQ=255)  ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGc                                       <  1:1127388‑M1/50‑2 (MQ=255)  ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGc                                       >  1:1458233‑M1/1‑49 (MQ=255)   cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcc                                      >  1:2727028‑M1/1‑48 (MQ=255)   cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcc                                      <  1:1750973‑M1/50‑3 (MQ=255)    aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcc                                      >  1:108004‑M1/1‑47 (MQ=255)    aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcc                                      >  1:4399886‑M1/1‑47 (MQ=255)     aaaCAGCTTTCACTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGccta                                    <  1:452889‑M1/50‑5 (MQ=255)      aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctac                                   <  1:3664988‑M1/50‑6 (MQ=255)      aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctac                                   <  1:4658811‑M1/50‑6 (MQ=255)      aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctac                                   <  1:3390152‑M1/50‑6 (MQ=255)      aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGccta                                    >  1:1678129‑M1/1‑45 (MQ=255)       aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctaca                                  >  1:266912‑M1/1‑44 (MQ=255)       aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctaca                                  <  1:3364360‑M1/50‑7 (MQ=255)        cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                 >  1:2409744‑M1/1‑43 (MQ=255)        cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                 >  1:586574‑M1/1‑43 (MQ=255)         aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                 >  1:2541503‑M1/1‑42 (MQ=255)         aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                 >  1:1869895‑M1/1‑42 (MQ=255)          gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacacta                               <  1:2903365‑M1/50‑10 (MQ=255)           cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactag                              <  1:4156108‑M1/50‑11 (MQ=255)             ttCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcg                            >  1:5245664‑M1/1‑38 (MQ=255)             ttCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcg                            >  1:4045368‑M1/1‑38 (MQ=255)              tCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcg                            <  1:4017334‑M1/49‑13 (MQ=255)              tCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagc                             <  1:4361151‑M1/48‑12 (MQ=255)               cGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctactctagcgg                           <  1:6961573‑M1/49‑14 (MQ=255)               cGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgg                           >  1:1508282‑M1/1‑36 (MQ=255)               cGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgg                           >  1:1616809‑M1/1‑36 (MQ=255)               cGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgg                           <  1:7584453‑M1/49‑14 (MQ=255)               cGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgg                           >  1:1973159‑M1/1‑36 (MQ=255)               cGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgg                           <  1:3129204‑M1/49‑14 (MQ=255)               cGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgg                           >  1:3529528‑M1/1‑36 (MQ=255)                gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                         >  1:2449427‑M1/1‑35 (MQ=255)                gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                         >  1:2431595‑M1/1‑35 (MQ=255)                gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                         <  1:3188943‑M1/50‑16 (MQ=255)                gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                         >  1:1895784‑M1/1‑35 (MQ=255)                gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcccacactagcgggt                         >  1:6579562‑M1/1‑35 (MQ=255)                 cTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggat                        <  1:35610‑M1/50‑17 (MQ=255)                 cTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggat                        >  1:5799137‑M1/1‑34 (MQ=255)                  tACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggatt                       >  1:880597‑M1/1‑33 (MQ=255)                   aCGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattc                      >  1:4668699‑M1/1‑32 (MQ=255)                    cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattct                     >  1:6059400‑M1/1‑31 (MQ=255)                    cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattc                      <  1:4443943‑M1/49‑19 (MQ=255)                    cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattc                      <  1:99256‑M1/49‑19 (MQ=255)                    cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattc                      <  1:1758127‑M1/49‑19 (MQ=255)                    cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggatt                       <  1:5521940‑M1/48‑18 (MQ=255)                    cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggatt                       <  1:2693422‑M1/48‑18 (MQ=255)                     gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                    >  1:3762801‑M1/1‑30 (MQ=255)                     gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                    >  1:6183876‑M1/1‑30 (MQ=255)                     gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                    >  1:5061401‑M1/1‑30 (MQ=255)                      ttGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttt                   <  1:6452521‑M1/50‑22 (MQ=255)                      ttGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttt                   <  1:2115862‑M1/50‑22 (MQ=255)                      ttGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttt                   <  1:2397470‑M1/50‑22 (MQ=255)                       tGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttt                  >  1:3693447‑M1/1‑28 (MQ=255)                        gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttttt                 >  1:3765121‑M1/1‑27 (MQ=255)                        gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttttt                 >  1:448306‑M1/1‑27 (MQ=255)                         cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttg                <  1:5044765‑M1/50‑25 (MQ=255)                         cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttg                <  1:2046138‑M1/50‑25 (MQ=255)                         cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttttt                 <  1:3539455‑M1/49‑24 (MQ=255)                          ttgctCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt               >  1:2722706‑M1/3‑25 (MQ=255)                          tGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt               >  1:6920844‑M1/1‑25 (MQ=255)                          tGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt               >  1:7574174‑M1/1‑25 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta             <  1:6289081‑M1/50‑28 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta             <  1:339104‑M1/50‑28 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta             <  1:2121567‑M1/50‑28 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt              >  1:2491610‑M1/1‑23 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt              <  1:6040690‑M1/49‑27 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt              >  1:2642665‑M1/1‑23 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt              >  1:7223118‑M1/1‑23 (MQ=255)                            gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt               <  1:590581‑M1/48‑26 (MQ=255)                               cGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacc          >  1:4970856‑M1/1‑20 (MQ=255)                                gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccc         >  1:181194‑M1/1‑19 (MQ=255)                                gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccc         >  1:6103036‑M1/1‑19 (MQ=255)                                gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccc         >  1:3748638‑M1/1‑19 (MQ=255)                                 ttttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccc        <  1:4367381‑M1/50‑33 (MQ=255)                                   ttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccccta      >  1:6011719‑M1/1‑16 (MQ=255)                                   ttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccct       <  1:5195878‑M1/49‑34 (MQ=255)                                    tAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccctac     >  1:6556170‑M1/1‑15 (MQ=255)                                    tAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccctac     >  1:4323959‑M1/1‑15 (MQ=255)                                    tAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccccta      <  1:1550015‑M1/49‑35 (MQ=255)                                     aaCACGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacc    >  1:1879101‑M1/1‑14 (MQ=255)                                     aaCACGGAG‑AAGTGcctacactagcgggattctttttgttaacccctac     <  1:5650662‑M1/49‑36 (MQ=255)                                      acacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctaccc   >  1:5792209‑M1/1‑13 (MQ=255)                                       cacGGAG‑AAGTGcctacgctagcgggattctttttgttaacccctacccc  <  1:577641‑M1/50‑39 (MQ=255)                                       cacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacccc  <  1:107503‑M1/50‑39 (MQ=255)                                              |                                            ACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAG  >  NC_000913/3090296‑3090384

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: ‑