Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	3,640,845	0	T→G	40.0%	10.2	31	Y30S (TAC→TCC)	yhjW	predicted metal dependent hydrolase
Reads supporting (aligned to +/- strand):  new base (12/0):  ref base (4/14):  total (16/15)

TAAAATCGTGCGCATAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAG  >  REL606/3640813‑3640880                                 |                                    tAAAATCGTGCGCATAGCTGCCGAAGCGGCGGTaaa                                  <  1:1350217/36‑1 (MQ=255)  aaaaTCGTGCGCAAAGCTGCCGAAGCGGCGGTaaaa                                 <  1:635617/36‑1 (MQ=255)   aaaTCGTGCGCATAGCTGCCGAAGCGGCGGTcaaaa                                >  1:6909557/1‑36 (MQ=255)    aaTCGTGCGCATAGCTGCCGAAGCGGCGGGagaaaa                               >  1:4471053/1‑36 (MQ=21)     aTCGTGCGCATAGCTGCCGAAGCGGCGGTAAAAAAc                              <  1:358114/36‑1 (MQ=255)     aTCGTGCGCATAGCTGCCGAAGCGGCGGGgaaaaaa                              >  1:6719528/1‑35 (MQ=21)      tCGTGGGCATACCTGCCGAAGCGGCGGTAAAAAAcc                             <  1:4212378/36‑1 (MQ=25)       cGTGCGCATAGCTGCCGAAGCGGCGGGaaaaaaaca                            >  1:3485871/1‑33 (MQ=21)         tGCGCATAGCTGCCGAAGCGGCGGGAAAAAACcggc                          >  1:2900446/1‑34 (MQ=21)         tGCGCAGAGCTGCCGAAGCGGCGGTAAAAAACcgcg                          <  1:7470977/36‑1 (MQ=255)          gcgcATAGCTGCCGAAGCGGCGGTAAAAAACCgcgc                         <  1:2221784/36‑1 (MQ=255)            gcATAGCTGCCGAAGCGGCGGGAAAAAACCGCGccc                       >  1:7112704/1‑35 (MQ=21)             cATAGCTGCCGAAGCGGCGGGAAAAAAACACGccct                      >  1:768772/1‑34 (MQ=17)              aTAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCAtt                     <  1:303538/36‑1 (MQ=255)              aTAGCTGCCGAAGCGGCGGGAAAAAACCACGcccat                     >  1:1812463/1‑33 (MQ=18)                aGCTGCCGAAGCGGCGGGAAAAAAACGCGCCAttta                   >  1:6322939/1‑34 (MQ=18)                aGCTGCCGAAGCGGCGGGAAAAAAACCCGCCCTTCa                   >  1:1708889/1‑36 (MQ=17)                 ccTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCAt                  <  1:1621391/35‑1 (MQ=255)                  ctcccGAACCGGCGGTAAAAAACCGCGCCATTCATa                 <  1:3233866/33‑1 (MQ=25)                  cTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATa                 <  1:2664851/36‑1 (MQ=255)                   tGCCGAAGCGGCGGGAAAAAACCGCGCCACTCATaa                >  1:5542233/1‑36 (MQ=21)                     ccGAAGCGGCGGGGAAAAACCGCGCCATTCATaaaa              >  1:2657731/1‑36 (MQ=25)                      cGAAGCGGCGGTAAAAAACCGCGCCATTCATaaaaa             <  1:1616655/36‑1 (MQ=255)                      cGAAGCGGCGGGAAAAAACCGCGCCCTTCATaaaaa             >  1:6965662/1‑36 (MQ=25)                        aaGCGGCGGTCAAAAACCGCGCCATTCATAAAAAgg           >  1:3822159/1‑36 (MQ=255)                          gcggcgGTAAAAAACCGCGCCATTCATAAAAAGGcc         <  1:1797589/36‑1 (MQ=255)                            ggcggTAAAAAACCGCGCCATTCATAAAAAGGCCaa       >  1:7228648/1‑36 (MQ=255)                            ggcggAAAAAAACCGCGCCATTCATAAAAAGGCCaa       <  1:725684/36‑1 (MQ=255)                             gcggTAAAAAACCGCGCCATTCATAAAAAAGCCAat      >  1:1732163/1‑36 (MQ=255)                              gggTAAAAAACCGCGCCATTCATAAAAAGGCCAata     <  1:1934250/35‑1 (MQ=255)                              cagTAAAAAACCGCGCCATTCATAAAAAGGCCAata     <  1:1882928/34‑1 (MQ=255)                                 tAAAAAACCGCGCCATTCATAAAAAGGCCAATATAg  >  1:2499804/1‑36 (MQ=255)                                 |                                    TAAAATCGTGCGCATAGCTGCCGAAGCGGCGGTAAAAAACCGCGCCATTCATAAAAAGGCCAATATAG  >  REL606/3640813‑3640880 Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 10 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG