Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,202,568	C→T	intergenic (‑9/‑191)	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,568	0	C	T	100.0%	145.4 / NA	40	intergenic (‑9/‑191)	iclR/metH	DNA‑binding transcriptional repressor/B12‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/22);  total (18/22)

CCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATCC  >  REL606/4202495‑4202637                                                                          |                                                                      ccAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGATGGGGGCGGCCCTGACAGtcttct                                                                     >  1:1503482/1‑76 (MQ=255)  cAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTcttctt                                                                    <  2:1190911/76‑1 (MQ=255)   aGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCttt                                                                   <  2:1698039/76‑1 (MQ=255)    gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCtttt                                                                  <  1:1812961/76‑1 (MQ=255)     tCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCttttt                                                                 <  2:549695/76‑1 (MQ=255)     tCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGGCCATGCCAGGCtttttttt                                                                 >  2:862504/1‑76 (MQ=255)      cGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGGCCCAGGCAGGCTTTTTTTTc                                                                >  2:1872512/1‑76 (MQ=255)       gCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGAAGGGGGGGACCATGACAGGCTTCttttttt                                                               >  1:1379967/1‑74 (MQ=255)         tGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTTGCGGGGATGGGGGGGGCCATGACAGTCTTCTTTTTTTGt                                                             >  1:437231/1‑76 (MQ=255)          ggTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCATGACAGTCTTCTTTTTTTGTa                                                            >  2:1112985/1‑76 (MQ=255)            tGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATc                                                          <  2:203949/76‑1 (MQ=255)             gCGGTGGCAACGNCGNGTTTTCTGCCGCGTTTCGNGGGAATAGGTGCGACCATGACAGTCTTCTTTTTCTGTATCg                                                         <  1:1657382/76‑1 (MQ=255)                 tGGCAGCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGa                                                     <  1:959288/76‑1 (MQ=255)                 tGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCATGACAGTCTTCTTTTTTTTTTTCGTGGa                                                     >  2:367408/1‑76 (MQ=255)                  ggCAGCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGaa                                                    <  1:1016091/76‑1 (MQ=255)                   gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGaaa                                                   <  2:1572896/76‑1 (MQ=255)                      aCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAAtca                                                <  2:1808813/76‑1 (MQ=255)                       cggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAAtcat                                               <  2:1450571/76‑1 (MQ=255)                         gcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAAtcattt                                             >  2:1873129/1‑76 (MQ=255)                          cggGTTTTCTGCCGCGTTTCGCGTGAANGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAAtcatttt                                            >  2:3981/1‑76 (MQ=255)                           gggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAAtcattttc                                           <  1:717810/76‑1 (MQ=255)                            ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAAtcattttca                                          <  2:369841/76‑1 (MQ=255)                              ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAAtcattttcatt                                        >  1:1153013/1‑76 (MQ=255)                               tttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAAtcattttcattt                                       <  2:466308/76‑1 (MQ=255)                                 tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCAttttt                                     <  2:996525/76‑1 (MQ=255)                                   tGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAt                                   <  2:698632/76‑1 (MQ=255)                                    gCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtt                                  >  2:454982/1‑76 (MQ=255)                                    gCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtt                                  >  2:1287771/1‑76 (MQ=255)                                    gCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtt                                  <  1:629894/76‑1 (MQ=255)                                       gcgTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGtt                               <  2:1444459/76‑1 (MQ=255)                                        cgTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTNTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTa                              <  1:1994840/76‑1 (MQ=255)                                         gTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAg                             >  2:1148232/1‑76 (MQ=255)                                          tttCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGc                            >  2:942312/1‑76 (MQ=255)                                            tCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTa                          <  2:356380/76‑1 (MQ=255)                                             cgcgGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTaa                         <  2:351108/76‑1 (MQ=255)                                               cgGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATg                       >  2:1193331/1‑76 (MQ=255)                                                  gAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCaa                    <  1:151894/76‑1 (MQ=255)                                                       ggTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGtt               <  1:362634/76‑1 (MQ=255)                                                          gCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCt            >  2:418936/1‑76 (MQ=255)                                                          gCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCt            >  1:333668/1‑76 (MQ=255)                                                            gACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGa          >  1:1903728/1‑76 (MQ=255)                                                             aCCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTGTTATTGTTAGCTAATGCAATAGTTGCTGaa         >  1:44246/1‑76 (MQ=255)                                                             aCCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCAGTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGaa         >  1:1711935/1‑76 (MQ=255)                                                               cATNACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACt       >  1:1839288/1‑76 (MQ=255)                                                               cATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACt       >  1:929117/1‑76 (MQ=255)                                                                 tGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGa     >  1:1300592/1‑76 (MQ=255)                                                                    cAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATcc  >  2:670127/1‑76 (MQ=255)                                                                    cAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATcc  >  1:1505555/1‑76 (MQ=255)                                                                          |                                                                      CCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATCC  >  REL606/4202495‑4202637

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