Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	402,056	0	T→G	31.2%	11.7	64	G60G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  new base (19/1):  ref base (17/27):  total (36/28)

CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACT  >  REL606/402021‑402091                                    |                                    ccACATGGCGAAGTACACGCGTTGCGTATGGCGGGt                                     <  1:7173035/36‑1 (MQ=255) ccACATGCCGAAGTACACGCGTTGCGTATGGCGGGt                                     <  1:4920493/36‑1 (MQ=255) ccACATGCCGAAGTACACGCGTTGCGTATGGCGGGt                                     <  1:6007242/36‑1 (MQ=255) ccACATGCCGAACTACACGCGTTGCGTATGGCGGGt                                     <  1:3469740/36‑1 (MQ=255)  cacaTGCCGAAGTACACGCGTTGCGTATGGCGGGTg                                    <  1:1115001/36‑1 (MQ=255)  cacaTGCCGAAGTACACGCGTTGCGTATGGCGGGTg                                    <  1:6306183/36‑1 (MQ=255)   ccaTGCCGAAGTACACGCGTTGCGTATGGCGGGTGa                                   <  1:4115811/35‑1 (MQ=255)   acaTGCCGAAGTACACGCGTTGCGTATGGCGGGTgg                                   >  1:3638531/1‑35 (MQ=255)    caTGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                  <  1:7212926/36‑1 (MQ=255)     aTGCCGTAGTACACGGGTTGCGTATGGCGGGTGaaa                                 <  1:6222936/36‑1 (MQ=25)     aTGCCGAAGTACACGCGTTGCGTATGGCGGGTGgaa                                 >  1:4258095/1‑36 (MQ=255)     aTGCCGAAGTACACGCGGTGCGTATGGCGGGTGaaa                                 <  1:5255871/36‑1 (MQ=255)      tGCCGAAGTACACGCGTTGCGTATGGCGGGTGaaaa                                <  1:3261147/36‑1 (MQ=255)      tGCCGAAGTACACGCGTTGCGTATGGCGGGGGgaaa                                >  1:5532704/1‑36 (MQ=255)       gCCGAAGTACACGCGTTGCGTATGGCGGGGGaaaaa                               >  1:2290762/1‑36 (MQ=255)        ccGAAGTACACGCTTTGCGTATGGCGGGTGAAAAAg                              <  1:2466786/36‑1 (MQ=255)        ccGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAg                              <  1:3818254/36‑1 (MQ=255)        ccGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAg                              >  1:4793634/1‑36 (MQ=25)        ccGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAg                              >  1:7189266/1‑36 (MQ=25)         cGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGc                             <  1:4929883/36‑1 (MQ=255)          gAAGTACACTCGTTGCGTATGGCGGGTGAAAAAGcc                            <  1:5973268/36‑1 (MQ=255)          gAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGcc                            <  1:3542241/36‑1 (MQ=255)          gAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGcc                            <  1:6304752/36‑1 (MQ=255)          gAAGTACACGCGTTGCGTATGGCGGGGGaaaaaaac                            >  1:5861669/1‑33 (MQ=25)          gAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGcc                            >  1:3359844/1‑36 (MQ=25)          gAAGTACACGCGTTGCGTATGGCGGGGGAAAAAgac                            >  1:4539829/1‑34 (MQ=255)           aaGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCa                           <  1:3221843/36‑1 (MQ=255)            aGTACACGCGTTGCGTATGGCGGGGGAAAAAACCaa                          >  1:862659/1‑36 (MQ=25)             gttcacGCGTTGCGTATGGCGGGTGAAAAAGCCaaa                         <  1:3021180/33‑1 (MQ=255)             gTACACGCGTTGCGTATGGCGGGGGGAAAAGCCcaa                         >  1:7406082/1‑36 (MQ=21)             gTACACGCGTTGCGTATGGCGGGGGGAAAAGCCaaa                         >  1:4307217/1‑36 (MQ=25)             gTACACGCGTTGCGTATGGCGGGGGAAAAAGCCaaa                         <  1:4884849/36‑1 (MQ=255)              tACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAg                        <  1:7352248/36‑1 (MQ=255)              tACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAg                        <  1:1478127/36‑1 (MQ=255)              tACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAg                        <  1:2018526/36‑1 (MQ=255)              tACACGCGTTGCGTATGGCGGGGGAAAAAGGCAAAg                        >  1:6806824/1‑36 (MQ=25)               acacGCGTTGCGTATGGCGGGGGGAAAAGTCAAAgg                       >  1:6834850/1‑36 (MQ=21)               acacGCGTTGCGTATGGCGGGGGAAAAAGACaaacg                       >  1:749264/1‑34 (MQ=21)                cacGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGt                      >  1:2465596/1‑36 (MQ=255)                cacGCGTAGCGTATGGCGGGTGAAAAAGCCAAAGGt                      <  1:6689140/36‑1 (MQ=255)                 acGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTg                     <  1:7223559/36‑1 (MQ=255)                 acGCGTTGCGTATGGCGGGTGAAAAAGACAAAgggg                     >  1:6424304/1‑34 (MQ=25)                 acGCGTTGCGTATGGCGGGGGGAAAAGGCAAAgggg                     >  1:6761567/1‑34 (MQ=18)                 acGCGTTGCGTATGGCGGGGGAAAAAACCCACGGTg                     >  1:3751393/1‑36 (MQ=18)                    cgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGa                  <  1:3702732/36‑1 (MQ=255)                    cgTTGCGTATGGCGGGGGGAAAAGGCAAAGGGGCgg                  >  1:4491030/1‑35 (MQ=18)                    cgTTGCGTATGGCGGGGGAAAAAGCCAAAGGCGCgg                  >  1:469555/1‑35 (MQ=21)                     gTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAc                 >  1:2556950/1‑36 (MQ=255)                     gTTGCGTATGGCGGGGGGGAAAGCCAAAGGTGGGAc                 >  1:6075486/1‑36 (MQ=18)                      ttGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAcc                <  1:2078410/36‑1 (MQ=255)                      ttGCGTATGGCGGGTGAAAAAGCCAAAGGGGGGAcc                >  1:2272928/1‑36 (MQ=21)                      tggCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAcc                <  1:407756/34‑1 (MQ=255)                       tGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCg               >  1:2637011/1‑36 (MQ=255)                       tGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCg               <  1:5382581/36‑1 (MQ=255)                       tGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGAACg               >  1:4848432/1‑36 (MQ=21)                         cGTATGGCGGGTGAAAAAACCAAAGGTGCGACCGcc             >  1:4907024/1‑36 (MQ=255)                          gTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCt            >  1:2276194/1‑36 (MQ=255)                          gTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCt            <  1:5014708/36‑1 (MQ=255)                          gTATGGCGGGTGAAAAAGCCAAAAGTGCGAACGCCt            >  1:3210910/1‑36 (MQ=21)                           tagggCGGGTGAAAAAGCCAAAGGTGCGACCGCCTa           <  1:5446779/33‑1 (MQ=255)                           tATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTa           <  1:4629880/36‑1 (MQ=255)                            aTGGCGGGTGAAAAAGCCAAAAGTGCGACCGCCTac          >  1:223237/1‑35 (MQ=25)                             tGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATg         >  1:2715832/1‑36 (MQ=255)                             tGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATg         >  1:2386002/1‑36 (MQ=255)                             tGGCGGGTGAAAAAGCCAAAGGAGGCACCCCCTATg         >  1:981093/1‑36 (MQ=17)                                cGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTca      >  1:4578459/1‑36 (MQ=255)                                  ggTGAAAAAGCCAAAGGTGCGACCGCCCATGCcaca    >  1:2274198/1‑36 (MQ=25)                                   gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCacac   >  1:915935/1‑36 (MQ=255)                                    tGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACt  >  1:107980/1‑36 (MQ=255)                                    |                                    CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACT  >  REL606/402021‑402091 Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG