Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	410,389	G→T	R140L (CGC→CTC)	thiI →	thiamine biosynthesis protein ThiI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	410,389	0	G	T	100.0%	187.3 / NA	48	R140L (CGC→CTC)	thiI	thiamine biosynthesis protein ThiI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/26);  total (22/26)

AACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCGCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAG  >  REL606/410317‑410460                                                                         |                                                                        aaCATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGt                                                                      >  2:686813/1‑76 (MQ=255)  aCATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTg                                                                     <  2:1551137/76‑1 (MQ=255)   cATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGa                                                                    <  1:1339607/76‑1 (MQ=255)   cATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGa                                                                    <  2:1645156/76‑1 (MQ=255)    aTGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGaa                                                                   <  1:1795228/76‑1 (MQ=255)     tGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAg                                                                  >  1:174969/1‑76 (MQ=255)       aTTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCt                                                                <  1:128125/76‑1 (MQ=255)        ttttAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTg                                                               <  1:319357/76‑1 (MQ=255)        ttttAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTg                                                               >  1:1861368/1‑76 (MQ=255)        ttttAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTg                                                               <  1:1763264/76‑1 (MQ=255)          ttAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGAc                                                             >  1:1560512/1‑76 (MQ=255)              cTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAAt                                                         <  1:613382/76‑1 (MQ=255)              cTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAAt                                                         >  2:820197/1‑76 (MQ=255)                cGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATcc                                                       >  2:1794587/1‑76 (MQ=255)                 gATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCg                                                      <  1:1094061/76‑1 (MQ=255)                  aTTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCgg                                                     >  1:1441494/1‑76 (MQ=255)                   ttGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGa                                                    >  1:1664259/1‑76 (MQ=255)                     gATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAg                                                  <  2:521857/76‑1 (MQ=255)                       tgtgGAACGTTACGTCGGCGGCGGTTTAAATTAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATTCGGAAGTg                                                >  2:1477825/1‑76 (MQ=255)                        gtgGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGa                                               <  2:83272/76‑1 (MQ=255)                         tgGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGAc                                              >  1:1105698/1‑76 (MQ=255)                         tgGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGAc                                              <  1:986362/76‑1 (MQ=255)                              cGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTcc                                         <  2:1396163/76‑1 (MQ=255)                                ttACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCAt                                       >  2:573379/1‑76 (MQ=255)                                 tACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATc                                      <  1:1728837/76‑1 (MQ=255)                                     tCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCtgga                                  >  2:610224/1‑76 (MQ=255)                                     tCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCtgga                                  <  1:1016996/76‑1 (MQ=255)                                      cggcggcggTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCtggaa                                 <  2:802738/76‑1 (MQ=255)                                      cggcggcggTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCtggaa                                 <  1:941843/76‑1 (MQ=255)                                        gcggcggTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCtggaagt                               >  1:1779363/1‑76 (MQ=255)                                          ggcggTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCtggaagtgg                             <  1:1375243/76‑1 (MQ=255)                                           gcggTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCtggaagtgga                            <  1:1375119/76‑1 (MQ=255)                                             ggTTTAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCtggaagtggaag                          >  2:1727029/1‑76 (MQ=255)                                               tttAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGAc                        >  1:913235/1‑76 (MQ=255)                                               tttAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGAc                        >  1:1632425/1‑76 (MQ=255)                                                ttAAATCAGCTTATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACg                       >  1:125986/1‑76 (MQ=255)                                                ttAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACg                       >  2:1244087/1‑76 (MQ=255)                                                 tAAATCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGa                      >  1:601411/1‑76 (MQ=255)                                                  aaaTCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGAt                     <  2:437134/76‑1 (MQ=255)                                                    aTCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCg                   >  1:793930/1‑76 (MQ=255)                                                    aTCAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCg                   <  2:482100/76‑1 (MQ=255)                                                      cAGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGtc                 >  1:235304/1‑76 (MQ=255)                                                       aGCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGtct                <  2:439895/76‑1 (MQ=255)                                                        gCATATTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGtctc               <  1:317410/76‑1 (MQ=255)                                                          atatTGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCct             <  2:1845208/76‑1 (MQ=255)                                                              tGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCctgctg         >  2:1405251/1‑76 (MQ=255)                                                              tGAGTCCGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCctgctg         <  1:1700117/76‑1 (MQ=255)                                                                   ccGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTaa    <  2:62731/76‑1 (MQ=255)                                                                   ccGCGCTCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTaa    >  2:920397/1‑76 (MQ=255)                                                                    cgcgctcgTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTaaa   >  1:244322/1‑76 (MQ=255)                                                                     gcgctcgTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAg  >  2:1506034/1‑76 (MQ=255)                                                                     gcgctcgTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAg  >  1:190567/1‑76 (MQ=255)                                                                         |                                                                        AACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAGTCCGCGCGCGTGAAGCTGACCAATCCGGAAGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAG  >  REL606/410317‑410460

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