Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,202,619	Δ1 bp	intergenic (‑60/‑140)	iclR ← / → metH	DNA‑binding transcriptional repressor/B12‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,202,618	0	A	.	100.0%	285.4 / NA	57	intergenic (‑59/‑141)	iclR/metH	DNA‑binding transcriptional repressor/B12‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (34/23);  total (34/23)

GGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCAC  >  REL606/4202549‑4202689                                                                      |                                                                        ggTGCGACCATGACGTTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTg                                                                  <  1:258263/76‑1 (MQ=255) ggTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTg                                                                  <  2:1247917/76‑1 (MQ=255)      gACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGaa                                                             >  2:907931/1‑76 (MQ=255)      gACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGaa                                                             >  1:197679/1‑76 (MQ=255)          aTGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGa                                                         >  2:87096/1‑76 (MQ=255)            taCAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATc                                                       >  2:322617/2‑76 (MQ=255)             aCAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATcc                                                      >  2:509651/1‑76 (MQ=255)               aGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGa                                                    <  2:1164184/76‑1 (MQ=255)               aGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGa                                                    >  2:531655/1‑76 (MQ=255)               aGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGa                                                    >  1:1343086/1‑76 (MQ=255)               aGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGa                                                    >  2:565678/1‑76 (MQ=255)                gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGAt                                                   <  2:1214786/76‑1 (MQ=255)                gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGAt                                                   <  2:1525065/76‑1 (MQ=255)                 tctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATg                                                  <  2:595840/76‑1 (MQ=255)                 tctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATg                                                  <  1:23939/76‑1 (MQ=255)                  ctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGa                                                 >  2:232386/1‑76 (MQ=255)                  ctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGa                                                 <  2:1196950/76‑1 (MQ=255)                   tcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAg                                                <  2:335238/76‑1 (MQ=255)                   tcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAg                                                <  1:968075/76‑1 (MQ=255)                   tcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAg                                                <  2:158843/76‑1 (MQ=255)                    ccTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTACTGC‑ATAGTTGCTGAACTGATCCGATGAGt                                               >  1:410845/1‑76 (MQ=255)                    ccTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGt                                               >  2:1284301/1‑76 (MQ=255)                    ccTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGt                                               >  2:357199/1‑76 (MQ=255)                    ccTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGt                                               >  1:777045/1‑76 (MQ=255)                    ccTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGt                                               >  2:735070/1‑76 (MQ=255)                     cTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGtt                                              >  1:870481/1‑76 (MQ=255)                     cTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGtt                                              <  2:1347371/76‑1 (MQ=255)                     cTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGtt                                              <  2:1484872/76‑1 (MQ=255)                     cTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGtt                                              <  2:1406002/76‑1 (MQ=255)                       ttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTaa                                            <  1:181573/76‑1 (MQ=255)                       ttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTaa                                            >  1:803420/1‑76 (MQ=255)                           cTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGtt                                        >  1:963682/1‑76 (MQ=255)                                tCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACa                                   <  1:1115498/76‑1 (MQ=255)                                       aaTCATTTTCATTTTTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCa                            >  1:1007611/1‑76 (MQ=255)                                              ttcattttTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCg                     >  1:113840/1‑76 (MQ=255)                                               tcattttTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgt                    >  1:659599/1‑76 (MQ=255)                                                cattttTATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtg                   <  2:898158/76‑1 (MQ=255)                                                  tttttATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggt                 <  1:763007/76‑1 (MQ=255)                                                  tttttATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggt                 >  1:212623/1‑76 (MQ=255)                                                  tttttATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggt                 <  1:18238/76‑1 (MQ=255)                                                   ttttATTGTTAGCTAATGC‑ATAGTTGCTGAACTGACCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggtg                >  1:1545881/1‑76 (MQ=255)                                                     ttATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTggtggt              >  2:283161/1‑76 (MQ=255)                                                     ttATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTggtggt              <  1:1055023/76‑1 (MQ=255)                                                     ttATTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTggtggt              >  2:1390911/1‑76 (MQ=255)                                                       aTTGTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCg            >  1:1095487/1‑76 (MQ=255)                                                          gTTAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCtt         <  2:101515/76‑1 (MQ=255)                                                            tAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCtttt       >  1:58727/1‑76 (MQ=255)                                                            tAGCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCtttt       >  1:206948/1‑76 (MQ=255)                                                              gCTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAc     >  2:1521815/1‑76 (MQ=255)                                                               cTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc    >  2:1019677/1‑76 (MQ=255)                                                               cTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc    >  1:148557/1‑76 (MQ=255)                                                               cTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc    <  1:1208697/76‑1 (MQ=255)                                                               cTAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc    >  1:1123685/1‑76 (MQ=255)                                                                tAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACca   >  1:962816/1‑76 (MQ=255)                                                                tAATGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACca   >  1:1306980/1‑76 (MQ=255)                                                                 aaTGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcac  <  2:328135/76‑1 (MQ=255)                                                                 aaTGC‑ATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcac  >  2:736438/1‑76 (MQ=255)                                                                      |                                                                        GGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCAC  >  REL606/4202549‑4202689

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