Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	187,404	T→C	13.3%	E416G (GAA→GGA)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	187,404	0	T	C	13.3%	197.2 / 19.0	75	E416G (GAA→GGA)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (41/24);  new base C (3/7);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 8.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGT  >  W3110S.gb/187338‑187473                                                                   |                                                                      aTTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCgg                                                                   <  1:3391436/71‑1 (MQ=255)    ggCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGaa                                                                 >  1:619169/1‑70 (MQ=255)    ggCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGaa                                                                 >  1:1202154/1‑70 (MQ=255)    ggCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGaa                                                                 >  1:145041/1‑70 (MQ=255)    ggCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGaa                                                                 >  1:4403809/1‑70 (MQ=255)    ggCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGaa                                                                 >  1:4124134/1‑70 (MQ=255)    ggCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGaa                                                                 <  1:1903043/70‑1 (MQ=255)    ggCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGaa                                                                 <  1:700831/70‑1 (MQ=255)    ggCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGaa                                                                 >  1:557661/1‑70 (MQ=255)    ggCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGaa                                                                 >  1:3316020/1‑70 (MQ=255)       aaTAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAatat                                                              >  1:1614279/1‑70 (MQ=255)             ccGCGCCGCACCGCTCCGGGGTGACTCAGAATGCTCAAAAACAGTTTTCGTGCTCCCgg                                                                   <  1:1684009/59‑1 (MQ=255)             ccGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTCCCgg                                                                   <  1:3433297/59‑1 (MQ=255)               gcgcCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAata                                                               >  1:5105544/1‑61 (MQ=255)               gcgcCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCg                                                     <  1:2426090/71‑1 (MQ=255)                  ccgcACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAatat                                                              <  1:4873976/59‑1 (MQ=255)                   cgcACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCgg                                                                   <  1:1963188/53‑1 (MQ=255)                   cgcACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGgttg                                                 >  1:2553003/1‑71 (MQ=255)                    gcACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGgttgt                                                >  1:3372847/1‑71 (MQ=255)                       ccGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGttgttg                                              >  1:4330679/1‑70 (MQ=255)                        cGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCa                                            <  1:3853878/71‑1 (MQ=255)                         gCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGgttg                                                 >  1:2048912/1‑65 (MQ=255)                         gCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGgttg                                                 >  1:595884/1‑65 (MQ=255)                          cTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAat                                                                >  1:4972482/1‑49 (MQ=255)                               gggTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTac                                                           <  1:4773962/49‑1 (MQ=255)                                ggTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTCCCgg                                                                   <  1:2736920/40‑1 (MQ=255)                                   gACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACgg                                 >  1:2882104/1‑71 (MQ=255)                                    aCGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGTGACggg                                <  1:662714/71‑1 (MQ=255)                                    aCGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACggg                                >  1:618107/1‑71 (MQ=255)                                     cGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATgtc                                      <  1:2123060/64‑3 (MQ=255)                                        aGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCg                                                     >  1:3958200/1‑46 (MQ=255)                                        aGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATg                            >  1:3138605/1‑71 (MQ=255)                                          aaTGCTCAAAAACAGTTTTCGTGCTCCCGGAATATTACACAg                                                       <  1:2950086/42‑1 (MQ=255)                                          aaTGCTCAAAAACAGTTTTCGTGCTCCCGGAATATTACACAg                                                       <  1:866668/42‑1 (MQ=255)                                           aTGCTCAAAAACAGTTTTCGTGCTCCCGGAATATTACACAGCGGttgttg                                              >  1:5037628/1‑50 (MQ=255)                                           aTGCTCAAAAACAGTTTTCGTGCTCCCGGAATATTACACAGCGGttgttg                                              >  1:4389066/1‑50 (MQ=255)                                            tGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTa                                                            <  1:1804372/35‑1 (MQ=255)                                            tGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTa                                                            <  1:2160012/35‑1 (MQ=255)                                               tCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGa                                          >  1:1960611/1‑50 (MQ=255)                                               tCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGa                                          >  1:2936246/1‑50 (MQ=255)                                                cAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACt                    >  1:3626184/1‑71 (MQ=255)                                                cAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACt                    >  1:923071/1‑71 (MQ=255)                                                cAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACt                    >  1:1725123/1‑71 (MQ=255)                                                 aaaaaCAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAAc                     <  1:2244659/69‑1 (MQ=255)                                                 aaaaaCAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTg                   <  1:4211604/71‑1 (MQ=255)                                                    aaCAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGgcg                >  1:1029632/1‑71 (MQ=255)                                                    aaCAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGgcg                <  1:3020818/71‑1 (MQ=255)                                                    aaCAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGgcg                >  1:2983421/1‑71 (MQ=255)                                                    aaCAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGgcg                >  1:4559988/1‑71 (MQ=255)                                                     aCAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGa                                    >  1:3147721/1‑50 (MQ=255)                                                      cAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCa              <  1:3623989/71‑1 (MQ=255)                                                       aGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAg             >  1:3832253/1‑71 (MQ=255)                                                         ttttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACgg                                 <  1:1918717/49‑1 (MQ=255)                                                         ttttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTgg                  <  1:3234511/64‑1 (MQ=255)                                                         ttttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTg                   <  1:2789651/63‑1 (MQ=255)                                                         ttttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTg                   <  1:4434925/63‑1 (MQ=255)                                                         ttttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTg                   <  1:3472055/63‑1 (MQ=255)                                                         ttttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTg                   <  1:1951616/63‑1 (MQ=255)                                                         ttttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAg             >  1:5061097/1‑69 (MQ=255)                                                         ttttCGTGCTCCCGGAATATTACACAGCGGTTGTTGCa                                            <  1:486072/38‑1 (MQ=255)                                                         ttttCGTGCTCCCGGAATATTACACAGCGGTGGTTGCa                                            <  1:9453/38‑1 (MQ=37)                                                          tttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCgt          >  1:1910856/1‑71 (MQ=255)                                                           ttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGc                               >  1:2954748/1‑49 (MQ=255)                                                           ttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGgc                 >  1:2744165/1‑63 (MQ=255)                                                           ttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGgc                 >  1:2225631/1‑63 (MQ=255)                                                           ttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGgc                 >  1:747321/1‑63 (MQ=255)                                                           ttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCgt          <  1:4978021/70‑1 (MQ=255)                                                           ttCGTGCATCCGGAATATTACACAGCGGTTGTTGCAGATgg                                       >  1:3269874/1‑41 (MQ=255)                                                            tCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCgtgg        >  1:1733429/1‑71 (MQ=255)                                                            tCGTGCTCCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCgtgg        >  1:701860/1‑71 (MQ=255)                                                             cGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCgtg         <  1:4144256/69‑1 (MQ=255)                                                              gTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCgtggtg      >  1:4597589/1‑71 (MQ=255)                                                              gTGATTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCgtggtg      >  1:1777336/1‑71 (MQ=255)                                                                 cTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAg   >  1:2796133/1‑71 (MQ=255)                                                                 cTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAg   >  1:1035176/1‑71 (MQ=255)                                                                   tCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGt  <  1:4477279/70‑1 (MQ=255)                                                                   |                                                                      ATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGT  >  W3110S.gb/187338‑187473

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: ‑