Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,074,318	G→T	38.7%	A39E (GCG→GAG)	ycdM ←	predicted monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,074,318	0	G	T	38.7%	93.3 / 52.7	62	A39E (GCG→GAG)	ycdM	predicted monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (22/16);  new base T (12/12);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATC  >  W3110S.gb/1074248‑1074380                                                                       |                                                               atcgaaatGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTgcggcg                                                                <  1:982218/71‑1 (MQ=255) atcgaaatGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTgcggcg                                                                <  1:4826247/71‑1 (MQ=255)   cgaaatGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGctcg                                                              >  1:3436572/1‑71 (MQ=255)   cgaaatGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGctcg                                                              >  1:4632241/1‑71 (MQ=255)          gTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGt                                                             >  1:3116905/1‑65 (MQ=255)          gTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTc                                                        >  1:2680341/1‑70 (MQ=255)          gTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTc                                                        >  1:5073025/1‑70 (MQ=255)          gTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCg                                                       <  1:4508617/71‑1 (MQ=255)           tAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGa                                                      <  1:181026/71‑1 (MQ=255)            aGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGa                                                      >  1:2320480/1‑70 (MQ=255)                cTCCGCTTTTTGCACAATGGCTTTATTGAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAAt                                                   >  1:589258/1‑69 (MQ=255)                cTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAAt                                                   >  1:4364764/1‑69 (MQ=255)                cTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGGGGGTCGAAATGa                                                 <  1:943315/71‑1 (MQ=255)                  ccGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGagc                                               >  1:1346357/1‑71 (MQ=255)                  ccGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGagc                                               >  1:2521147/1‑71 (MQ=255)                   cGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGagc                                               <  1:4946192/70‑1 (MQ=255)                   cGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGag                                                <  1:1019115/69‑1 (MQ=255)                    gCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTg                                                            >  1:5031182/1‑56 (MQ=255)                    gCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTg                                                            >  1:1145322/1‑56 (MQ=255)                    gCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTg                                                            >  1:1730045/1‑56 (MQ=255)                        tttGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCg                                         >  1:5033582/1‑71 (MQ=255)                         ttGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGagccagcc                                          >  1:4465570/1‑69 (MQ=255)                         ttGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGagccagcc                                          >  1:1573377/1‑69 (MQ=255)                          tGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCgt                                        <  1:902070/70‑1 (MQ=255)                          tGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgtt                                       >  1:2583669/1‑71 (MQ=255)                           gCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCgtt                                       <  1:4570968/70‑1 (MQ=255)                           gCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgttg                                      >  1:3739549/1‑71 (MQ=255)                           gCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgtt                                       <  1:2253078/70‑1 (MQ=255)                                aTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGc                                  >  1:245697/1‑70 (MQ=255)                                aTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTgg                                                           <  1:2862040/45‑1 (MQ=255)                                aTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGcc                                 <  1:68840/71‑1 (MQ=255)                                aTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGcc                                 <  1:4393947/71‑1 (MQ=255)                                aTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGcc                                 <  1:913655/71‑1 (MQ=255)                                 tGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGcc                                 >  1:387364/1‑70 (MQ=255)                                 tGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGcc                                 >  1:2088290/1‑70 (MQ=255)                                       tATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGagccagcc                                          <  1:2775726/55‑1 (MQ=255)                                       tATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGagccagcc                                          <  1:4483225/55‑1 (MQ=255)                                        aTTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGGCAATAGGt                          <  1:1839756/70‑1 (MQ=255)                                        aTTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGt                          <  1:5097853/70‑1 (MQ=255)                                         ttCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTAc                        >  1:2101535/1‑71 (MQ=255)                                           cAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGagcca                                             >  1:2021927/1‑48 (MQ=255)                                           cAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACg                       >  1:4254950/1‑70 (MQ=255)                                              ttCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATa                   >  1:157813/1‑71 (MQ=255)                                              ttCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATa                   >  1:2698028/1‑71 (MQ=255)                                               tCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgttg                                      >  1:2115420/1‑51 (MQ=255)                                               tCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgttg                                      >  1:2947382/1‑51 (MQ=255)                                                  aaGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGcc               <  1:3289973/71‑1 (MQ=255)                                                       cGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATAc                  >  1:874244/1‑63 (MQ=255)                                                       cGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATAc                  >  1:4711407/1‑63 (MQ=255)                                                        ggCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAAttt          <  1:3373097/70‑1 (MQ=255)                                                         gCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgttgtt                                    >  1:1755745/1‑43 (MQ=255)                                                          cATGTACTGCGGCTCGTGGGTCTAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtc        <  1:1757799/70‑1 (MQ=255)                                                          cATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCTTTGTTGCCAATAGGTACGAATACGCCAATTTtc        <  1:372974/70‑1 (MQ=255)                                                          cATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtc        <  1:2468082/70‑1 (MQ=255)                                                          cATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtca       <  1:500766/71‑1 (MQ=255)                                                          cATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtca       >  1:3195781/1‑71 (MQ=255)                                                          cATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtca       >  1:2750421/1‑71 (MQ=255)                                                           aTGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtcat      <  1:766789/71‑1 (MQ=255)                                                            tGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtc        <  1:2796581/68‑1 (MQ=255)                                                            tGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtc        <  1:704490/68‑1 (MQ=255)                                                            tGTACTGCGGCGCGTGGGGCGAAATGAGCCAGCCgttg                                      <  1:2581323/38‑1 (MQ=255)                                                              tACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAg                            >  1:3342016/1‑46 (MQ=255)                                                              tACTGCGGCGAGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtcatcat   >  1:5074514/1‑71 (MQ=255)                                                                cTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtcatcatc  <  1:1268880/70‑1 (MQ=255)                                                                cTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtcatcatc  >  1:138412/1‑70 (MQ=255)                                                                cTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtcatcatc  <  1:1297747/70‑1 (MQ=255)                                                                  gcggctCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGaa                     <  1:4881522/49‑1 (MQ=255)                                                                   cggcgCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGc                >  1:164231/1‑53 (MQ=255)                                                                   cggcgCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGc                >  1:1254077/1‑53 (MQ=255)                                                                       |                                                               ATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATC  >  W3110S.gb/1074248‑1074380

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