Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,442,070	A→G	T497A (ACG→GCG)	thiC →	thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,442,070	0	A	G	100.0%	232.2 / NA	67	T497A (ACG→GCG)	thiC	thiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/36);  total (31/36)

ACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAG  >  W3110S.gb/3442002‑3442135                                                                     |                                                                  aCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCg                                                                 >  1:3152843/1‑71 (MQ=255)     cTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             <  1:1507238/71‑1 (MQ=255)       tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:3554079/1‑69 (MQ=255)       tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:2621726/1‑69 (MQ=255)       tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:2329165/1‑69 (MQ=255)       tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:4649158/1‑69 (MQ=255)        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:1660876/1‑68 (MQ=255)        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:2496873/1‑68 (MQ=255)        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:2452539/1‑68 (MQ=255)        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:415071/1‑68 (MQ=255)        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:3379668/1‑68 (MQ=255)        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:3619189/1‑68 (MQ=255)        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             >  1:3554171/1‑68 (MQ=255)          cGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAagag                                                        <  1:94099/71‑1 (MQ=255)           gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGc                                                       <  1:3969478/71‑1 (MQ=255)            tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCa                                                      <  1:3886991/71‑1 (MQ=255)             cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCa                                                      <  1:3135122/70‑1 (MQ=255)                  aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTggg                                                <  1:3615788/71‑1 (MQ=255)                  aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTggg                                                <  1:3117073/71‑1 (MQ=255)                  aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTggg                                                >  1:3812905/1‑71 (MQ=255)                    tGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTc                                              >  1:1434695/1‑71 (MQ=255)                        gcggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGcc                                          >  1:956228/1‑71 (MQ=255)                          ggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa                                                             <  1:713939/50‑1 (MQ=255)                             gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAAt                                      >  1:540980/1‑70 (MQ=255)                             gatgatTGGCTGGTTTGGCGGCGCGATGCTCTGTTACGTAGCGCCAAAagag                                                        <  1:4445127/52‑1 (MQ=39)                             gatgatTGGCTGGTTTAGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAAt                                      >  1:3501713/1‑70 (MQ=255)                                gatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATa                                     >  1:21757/1‑68 (MQ=255)                                gatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATa                                     >  1:3141550/1‑68 (MQ=255)                                 atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGt                                               >  1:3544077/1‑57 (MQ=255)                                  ttGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                <  1:479098/71‑1 (MQ=255)                                  ttGGCTGGTTTGGCTGCGCGATGCTCTCTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                <  1:4090837/71‑1 (MQ=255)                                   tGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                <  1:4594987/70‑1 (MQ=255)                                    tccTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                              <  1:2190995/69‑1 (MQ=255)                                     gCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTgg                                                 <  1:98829/51‑1 (MQ=255)                                     gCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                              <  1:58397/70‑1 (MQ=255)                                     ccTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                              <  1:3202061/69‑1 (MQ=255)                                         gTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCa                                                      <  1:3995750/42‑1 (MQ=255)                                         gTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCa                                                      <  1:2092997/42‑1 (MQ=255)                                         gTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGtt                          >  1:3423124/1‑70 (MQ=255)                                         gTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGtt                          >  1:1331219/1‑70 (MQ=255)                                         gTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGtt                          >  1:3306395/1‑70 (MQ=255)                                         gTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                         <  1:3863546/71‑1 (MQ=255)                                           ttGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATc                                                    >  1:3120884/1‑42 (MQ=255)                                           ttGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATc                                                    >  1:127531/1‑42 (MQ=255)                                                  cgcgATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgggg                 <  1:2816582/70‑1 (MQ=255)                                                  cgcgATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgggg                 <  1:4094809/70‑1 (MQ=255)                                                   gcgATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCt               <  1:3343549/71‑1 (MQ=255)                                                   gcgATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCt               <  1:4486789/71‑1 (MQ=255)                                                   gcgATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCt               <  1:4241395/71‑1 (MQ=255)                                                    cgATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAATATGTTAAGCAGGGGCtt              <  1:645237/71‑1 (MQ=255)                                                    cgATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCtt              <  1:1909540/71‑1 (MQ=255)                                                    cgATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCt               >  1:3730550/1‑70 (MQ=255)                                                    cgATGCTCTCTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCtt              <  1:1245180/71‑1 (MQ=255)                                                       tGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATc           <  1:4014880/71‑1 (MQ=255)                                                        gCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATg                            >  1:1121747/1‑53 (MQ=255)                                                        gCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCa          >  1:4113719/1‑71 (MQ=255)                                                              ttACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCtat     <  1:2324573/70‑1 (MQ=255)                                                              ttACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCta      <  1:152278/69‑1 (MQ=255)                                                               tACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa                     >  1:3273678/1‑53 (MQ=255)                                                               tACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCta      >  1:2470676/1‑68 (MQ=255)                                                                aCGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa                     >  1:3258471/1‑52 (MQ=255)                                                                aCGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATaa   <  1:3670638/70‑1 (MQ=255)                                                                aCGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATaa   <  1:3766633/70‑1 (MQ=255)                                                                aCGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATaa   <  1:1072016/70‑1 (MQ=255)                                                                 cGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAg  <  1:3750958/70‑1 (MQ=255)                                                                 cGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAg  <  1:111760/70‑1 (MQ=255)                                                                 cGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAg  <  1:1078625/70‑1 (MQ=255)                                                                     |                                                                  ACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAG  >  W3110S.gb/3442002‑3442135

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