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

CCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAG  >  NC_000913/3090297‑3090384                                             |                                            ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGc                                       <  2:3902350‑M1/50‑2 (MQ=255)  cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcc                                      <  2:348122‑M1/50‑3 (MQ=255)   aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcc                                      >  2:2935998‑M1/1‑47 (MQ=255)    aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGccta                                    <  1:4155597‑M1/50‑5 (MQ=255)    aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGccta                                    <  1:2448188‑M1/50‑5 (MQ=255)     aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGccta                                    >  2:878356‑M1/1‑45 (MQ=255)     aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGccta                                    >  2:1927830‑M1/1‑45 (MQ=255)     aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGccta                                    >  1:2422594‑M1/1‑45 (MQ=255)      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctaca                                  >  1:3569221‑M1/1‑44 (MQ=255)      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctaca                                  <  2:1492102‑M1/50‑7 (MQ=255)      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctaca                                  >  2:2439245‑M1/1‑44 (MQ=255)      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctaca                                  >  2:644009‑M1/1‑44 (MQ=255)       cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                 <  1:3214995‑M1/50‑8 (MQ=255)       cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                 <  2:2734223‑M1/50‑8 (MQ=255)       cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                 >  1:1897‑M1/1‑43 (MQ=255)       cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                 <  2:4492946‑M1/50‑8 (MQ=255)        aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacact                                <  1:1965851‑M1/50‑9 (MQ=255)        aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacact                                <  1:4189725‑M1/50‑9 (MQ=255)        aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacact                                <  2:4034462‑M1/50‑9 (MQ=255)         gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacacta                               <  1:2553955‑M1/50‑10 (MQ=255)         gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacac                                 >  2:3983319‑M1/1‑41 (MQ=255)          cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactag                              >  1:2340652‑M1/1‑40 (MQ=255)            ttCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcg                            >  1:3911150‑M1/1‑38 (MQ=255)             tCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcg                            <  1:3610553‑M1/49‑13 (MQ=255)              cGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgg                           >  2:3999867‑M1/1‑36 (MQ=255)              cGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGAGcctacactagcgg                           >  1:1559192‑M1/1‑36 (MQ=255)               gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                         >  1:1657940‑M1/1‑35 (MQ=255)               gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                         >  1:1180895‑M1/1‑35 (MQ=255)               gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                         <  1:4164229‑M1/50‑16 (MQ=255)               gCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggg                          <  1:2285237‑M1/49‑15 (MQ=255)               gCTACGTTGCCGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcggga                         >  2:4483216‑M1/1‑35 (MQ=255)                cTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggat                        <  1:2520275‑M1/50‑17 (MQ=255)                 tACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggatt                       >  2:3973808‑M1/1‑33 (MQ=255)                  aCGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattc                      >  2:1447460‑M1/1‑32 (MQ=255)                   cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattct                     <  2:1307332‑M1/50‑20 (MQ=255)                   cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattct                     >  1:4469206‑M1/1‑31 (MQ=255)                   cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattct                     >  1:2154744‑M1/1‑31 (MQ=255)                   cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattc                      <  2:3538856‑M1/49‑19 (MQ=255)                   cGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggatt                       <  1:4277220‑M1/48‑18 (MQ=255)                    gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                    >  1:2605704‑M1/1‑30 (MQ=255)                    gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                    >  2:3379173‑M1/1‑30 (MQ=255)                    gTTGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctt                    >  1:383334‑M1/1‑30 (MQ=255)                     ttGCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttt                   <  2:425257‑M1/50‑22 (MQ=255)                       gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttttt                 >  1:3362479‑M1/1‑27 (MQ=255)                       gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttt                  <  2:1431663‑M1/49‑23 (MQ=255)                       gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttt                  <  1:813068‑M1/49‑23 (MQ=255)                       gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttt                  <  1:2129434‑M1/49‑23 (MQ=255)                       gCTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcgttt                  <  2:688670‑M1/49‑23 (MQ=255)                        cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttg                <  1:3520269‑M1/50‑25 (MQ=255)                        cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttg                <  2:2870325‑M1/50‑25 (MQ=255)                        cTGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattcttttt                 <  1:2910894‑M1/49‑24 (MQ=255)                         tGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt               >  1:828896‑M1/1‑25 (MQ=255)                         tGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt               >  1:2591762‑M1/1‑25 (MQ=255)                         tGGCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt               >  1:4162708‑M1/1‑25 (MQ=255)                          ggCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgt               >  1:1169759‑M1/1‑24 (MQ=255)                           gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta             <  2:1761619‑M1/50‑28 (MQ=255)                           gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta             <  1:1182492‑M1/50‑28 (MQ=255)                           gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta             <  2:2363083‑M1/50‑28 (MQ=255)                           gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtta             <  1:4143074‑M1/50‑28 (MQ=255)                           gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt              <  2:3868375‑M1/49‑27 (MQ=255)                           gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt              >  1:4323380‑M1/1‑23 (MQ=255)                           gCTCGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgtt              >  2:623690‑M1/1‑23 (MQ=255)                              cGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacc          >  2:3943349‑M1/1‑20 (MQ=255)                              cGTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacc          >  1:2922300‑M1/1‑20 (MQ=255)                               gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccc         >  1:3128499‑M1/1‑19 (MQ=255)                               gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccc         >  1:1148559‑M1/1‑19 (MQ=255)                               gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccc         >  2:1098212‑M1/1‑19 (MQ=255)                               gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacca         <  2:4213602‑M1/50‑32 (MQ=255)                               gTTTTAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacc          <  1:159920‑M1/49‑31 (MQ=255)                                ttttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccc        <  1:405314‑M1/50‑33 (MQ=255)                                ttttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccc        >  1:2861200‑M1/1‑18 (MQ=255)                                 tttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccc        >  1:2151165‑M1/1‑17 (MQ=255)                                 tttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccc        >  1:4148242‑M1/1‑17 (MQ=255)                                 tttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccc        >  1:3702093‑M1/1‑17 (MQ=255)                                  ttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccct       <  1:4026803‑M1/49‑34 (MQ=255)                                  ttAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccct       <  2:225040‑M1/49‑34 (MQ=255)                                   tAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccctac     >  2:1430960‑M1/1‑15 (MQ=255)                                   tAACACGGAG‑AAGTGcctacactagcgggattctttttgttaacccctac     >  1:4151199‑M1/1‑15 (MQ=255)                                   tAACACGGAG‑AAGTGcctacactagcgggattctttttgttaaccccta      <  1:1956378‑M1/49‑35 (MQ=255)                                    aaCACGGAG‑AAGTGcctacactagcgggattctttttgttaacccctac     <  2:337420‑M1/49‑36 (MQ=255)                                    aaCACGGAG‑AAGCGcctacactagcgggattctttttgttaacccctacc    >  2:233466‑M1/1‑14 (MQ=255)                                     acacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctaccc   >  1:931830‑M1/1‑13 (MQ=255)                                     acacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctaccc   <  2:4085741‑M1/50‑38 (MQ=255)                                     acacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctaccc   >  1:3136347‑M1/1‑13 (MQ=255)                                     acacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacc    <  2:109994‑M1/49‑37 (MQ=255)                                      cacGGAG‑AAGTGcctacgctagcgggattctttttgttaacccctacccc  <  2:1434846‑M1/50‑39 (MQ=255)                                      cacGGAG‑AAGTGcctacgctagcgggattctttttgttaacccctacccc  <  2:150642‑M1/50‑39 (MQ=255)                                      cacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacccc  <  2:4236771‑M1/50‑39 (MQ=255)                                      cacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctacccc  <  2:737149‑M1/50‑39 (MQ=255)                                      cacGGAG‑AAGTGcctacactagcgggattctttttgttaacccctaccc   >  2:3002847‑M1/1‑12 (MQ=255)                                             |                                            CCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAG  >  NC_000913/3090297‑3090384

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