Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,894,681:1	+A	6.8%	coding (531/780 nt)	ygcR ←	predicted flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,894,681	1	.	A	6.8%	314.4 / 22.5	88	coding (531/780 nt)	ygcR	predicted flavoprotein
Reads supporting (aligned to +/- strand):  ref base . (52/30);  new base A (2/4);  total (54/34)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CCATTCGCTGGCGCATTCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAGCGCGTC  >  W3110S.gb/2894611‑2894746                                                                        |                                                                  ccATTCGCTGGCGCATTCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGc                                                                    <  1:1812181/71‑1 (MQ=255)      cGCTGGCGCATTCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGc                                                               >  1:2241440/1‑70 (MQ=255)      cGCTGGCGCATTCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGc                                                               >  1:7173/1‑70 (MQ=255)      cGCTGGCGCATTCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGc                                                               >  1:1361349/1‑70 (MQ=255)      cGCTGGCGCATTCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGc                                                               >  1:1740368/1‑70 (MQ=255)           gcgcATTCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAAccc                                                         <  1:3861540/71‑1 (MQ=255)           gcgcATTCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAAccc                                                         <  1:3707535/71‑1 (MQ=255)                ttCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGCAAAGCGaa                                                            <  1:330361/64‑1 (MQ=255)                ttCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGcaa                                                    >  1:1815704/1‑71 (MQ=255)                  cccGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGcaac                                                   >  1:26275/1‑70 (MQ=255)                   ccGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGCAAAGCGAACCCGGcaac                                                   <  1:2972513/70‑1 (MQ=255)                   ccGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGcaaca                                                  <  1:1186654/70‑1 (MQ=255)                     ggTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAAcc                                                          >  1:3268777/1‑60 (MQ=255)                     ggTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACg                                               >  1:3327201/1‑71 (MQ=255)                        aCTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCaa                                            >  1:4294790/1‑71 (MQ=255)                         cTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCaa                                            >  1:1913025/1‑70 (MQ=255)                         cTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAAc                                           <  1:4210582/71‑1 (MQ=255)                          tGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCg                                                              <  1:243274/51‑1 (MQ=255)                          tGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGcaac                                                   <  1:4617585/62‑1 (MQ=255)                          tGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGcaac                                                   <  1:2330162/62‑1 (MQ=255)                           ggCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGcaa                                                    >  1:2785955/1‑60 (MQ=255)                           ggCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAAcc                                          >  1:4551637/1‑70 (MQ=255)                            gCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCaacaac                                                >  1:3728415/1‑63 (MQ=255)                            gCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAAccc                                         >  1:3155823/1‑70 (MQ=255)                            gCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAAcc                                          <  1:192414/69‑1 (MQ=255)                            gCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCa                                        >  1:3367701/1‑71 (MQ=255)                               aGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGt                                     <  1:3703901/71‑1 (MQ=255)                                ggCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGtt                                    <  1:1505844/71‑1 (MQ=255)                                ggCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGtt                                    <  1:2557912/71‑1 (MQ=255)                                ggCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGtt                                    <  1:3989051/71‑1 (MQ=255)                                 gCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTc                                   <  1:4570432/71‑1 (MQ=255)                                  ccACTTCACCGCACTGACGCACGGCAATCACCGCAGGCAAAGCGAACCCGGCAACAACGCAACCCATGtt                                    >  1:2254157/1‑70 (MQ=255)                                    aCTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGtt                                    >  1:3611473/1‑67 (MQ=255)                                     cTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGcaac                                                   <  1:4271334/51‑1 (MQ=255)                                      ttCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCg                                                              <  1:2908626/39‑1 (MQ=255)                                      ttCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGtt                                    <  1:2244135/65‑1 (MQ=255)                                      ttCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAg                                 >  1:1978084/1‑68 (MQ=255)                                      ttCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAg                                 >  1:1972881/1‑68 (MQ=255)                                       tCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAg                                                                >  1:4315157/1‑36 (MQ=255)                                       tCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTa                               >  1:3204079/1‑69 (MQ=255)                                       tCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACg                             >  1:2146920/1‑71 (MQ=255)                                        cACCGCACTGACGCACTGCAATCACCGCAGTC‑AAGCGAACCCGGCAACAACGc                                              >  1:669832/1‑53 (MQ=255)                                        cACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCAt                                       >  1:2436589/1‑60 (MQ=255)                                        cACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCAt                                       >  1:472371/1‑60 (MQ=255)                                        cACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCa                                  >  1:364045/1‑65 (MQ=255)                                           cGCACTTACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGt                         >  1:4626240/1‑71 (MQ=255)                                           cGCACTGACGCACGGCAATCACCGCAGGCAAAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTg                          <  1:1472445/71‑1 (MQ=255)                                           cGCACTGACGCACGGCAATCACCGCAGGCAAAGCGAACCCGGCAACAACGCAACCCATGTTCAGCACGCTg                          <  1:601283/71‑1 (MQ=255)                                           cGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGt                         >  1:3454131/1‑71 (MQ=255)                                           cGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGt                         <  1:2685737/71‑1 (MQ=255)                                           cGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGt                         <  1:3359716/71‑1 (MQ=255)                                            gCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGcaacaa                                                 >  1:1708508/1‑46 (MQ=255)                                             cACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTc                       <  1:1935481/71‑1 (MQ=255)                                               cTGACGCACGTCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGa                     <  1:3780485/71‑1 (MQ=255)                                                    gCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTc                       <  1:1306555/64‑1 (MQ=255)                                                      aCGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGa               <  1:206190/70‑1 (MQ=255)                                                      aCGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGa               >  1:1901777/1‑70 (MQ=255)                                                       cGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGc                                              <  1:4176182/38‑1 (MQ=255)                                                       cGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGAt              >  1:1745435/1‑70 (MQ=255)                                                        ggCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGc                                              <  1:1552878/37‑1 (MQ=255)                                                        ggCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATa             <  1:527192/70‑1 (MQ=255)                                                         gCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATaaa           >  1:2418311/1‑71 (MQ=255)                                                          cAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATaaa           >  1:2951441/1‑70 (MQ=255)                                                          cAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAAc          <  1:642070/71‑1 (MQ=255)                                                          cAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAAc          <  1:4294105/71‑1 (MQ=255)                                                             tCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTAc                              <  1:1676499/48‑1 (MQ=255)                                                             tCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTAc                              <  1:4240479/48‑1 (MQ=255)                                                              cACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTa                               >  1:165478/1‑46 (MQ=255)                                                              cACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAgcg      >  1:4235877/1‑71 (MQ=255)                                                              cACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAgc       >  1:102920/1‑70 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGt                                >  1:2554729/1‑44 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTa                               >  1:1469162/1‑45 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTa                               >  1:3718978/1‑45 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTa                               >  1:1463798/1‑45 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTa                               >  1:217616/1‑45 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGc                            >  1:4323572/1‑48 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCt                           >  1:279114/1‑49 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTTTTCGAGGGTGATAAACAgcgc     >  1:3626545/1‑71 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATaaa           >  1:724183/1‑65 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAgcgc     >  1:302566/1‑71 (MQ=255)                                                               aCCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAgcgc     >  1:2935980/1‑71 (MQ=255)                                                                  gCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATaaa           >  1:4228792/1‑62 (MQ=255)                                                                  gCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAg        >  1:4293659/1‑65 (MQ=255)                                                                  gCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAg        >  1:1948905/1‑65 (MQ=255)                                                                  gCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAGCGCGTc  >  1:496101/1‑71 (MQ=255)                                                                  gCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGATGTTCGAGGGTGATAAACAg        >  1:3012794/1‑65 (MQ=255)                                                                   cAGGCAAAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAGCGCGTc  >  1:134179/1‑71 (MQ=255)                                                                     ggC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGa                     >  1:1848988/1‑49 (MQ=255)                                                                      gC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAggg                  <  1:3409499/51‑1 (MQ=255)                                                                       c‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAggg                  <  1:4171496/50‑1 (MQ=255)                                                                        |                                                                  CCATTCGCTGGCGCATTCCCGGTACTGGCAAGGCCACTTCACCGCACTGACGCACGGCAATCACCGCAGGC‑AAGCGAACCCGGCAACAACGCAACCCATGTTCAGTACGCTGTTCGAGGGTGATAAACAGCGCGTC  >  W3110S.gb/2894611‑2894746

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