Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,340,671	G→A	21.9%	T347T (ACC→ACT)	adiC ←	arginine:agmatin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,340,671	0	G	A	21.9%	175.4 / 44.9	82	T347T (ACC→ACT)	adiC	arginine:agmatin
Reads supporting (aligned to +/- strand):  ref base G (29/35);  new base A (12/6);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

GTGTAAAGATATGGCACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAATCAGCCCCGC  >  W3110S.gb/4340609‑4340739                                                               |                                                                     gtgtAAAGATATGGCACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTAGTCGCGtt                                                              <  1:517256/71‑1 (MQ=255)    tAAAGATATGGCACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTg                                                            <  1:3302374/70‑1 (MQ=255)    tAAAGATATGGCACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTg                                                            <  1:293099/70‑1 (MQ=255)    tAAAGATATGGCACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTg                                                            <  1:2186694/70‑1 (MQ=255)      aaGATATGGCACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTAGTCGCGTTTGGt                                                          >  1:385435/1‑70 (MQ=255)         atatGGCACCAGTGTAAAGATGACCGAGACGTAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGaaa                                                      >  1:1253805/1‑71 (MQ=255)            tGGCACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAAtgc                                                   <  1:1949604/71‑1 (MQ=255)             ggCACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAAtgct                                                  >  1:2480620/1‑71 (MQ=255)               cACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAAtgctgc                                                <  1:338995/71‑1 (MQ=255)                   aGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAg                                            <  1:681067/71‑1 (MQ=255)                       tAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTg                                                                      >  1:4645925/1‑41 (MQ=255)                       tAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAAtgc                                                   <  1:3440017/60‑1 (MQ=255)                       tAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGctgctg                                              >  1:2065808/1‑65 (MQ=255)                       tAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTgg                                        <  1:4303371/71‑1 (MQ=255)                       tAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAATGCTGCTGAg                                            >  1:3561955/1‑67 (MQ=255)                       tAAAGATGACCGAGACGGAAGAAACCAGACC‑AACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGa                                       <  1:1632150/71‑1 (MQ=255)                         aaGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGt                                                               >  1:2995469/1‑46 (MQ=255)                         aaGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGt                                                               >  1:2068916/1‑46 (MQ=255)                         aaGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGt                                                               >  1:2122323/1‑46 (MQ=255)                         aaGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGaa                                      <  1:4582073/71‑1 (MQ=255)                            aTGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGaa                                      <  1:2048160/68‑1 (MQ=255)                            aTGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGaa                                      <  1:561238/68‑1 (MQ=255)                            aTGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGa                                    <  1:2384340/70‑1 (MQ=255)                             tGACCGAGACGGAAGAAACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATg                                  >  1:564275/1‑71 (MQ=255)                               aCCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGa                                                        <  1:3164706/47‑1 (MQ=255)                               aCCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATgg                                 <  1:1118774/70‑1 (MQ=255)                               aCCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATgg                                 <  1:2330573/70‑1 (MQ=255)                                ccGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGTGTTTGGTGAAATGCTGCTGAGCTGGAAGATg                                  >  1:4501408/1‑68 (MQ=255)                                ccGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATg                                  >  1:731856/1‑68 (MQ=255)                                ccGAGACGGAAGAAACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAAtgct                                                  <  1:331627/52‑1 (MQ=255)                                    gaCGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATgg                                 >  1:1371615/1‑65 (MQ=255)                                    gaCGGAAGAAACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGtcatc                            >  1:1534848/1‑70 (MQ=255)                                      cGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGa                                    <  1:3792565/60‑1 (MQ=255)                                      cGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGa                                    <  1:4311725/60‑1 (MQ=255)                                       ggAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGa                                                        <  1:189383/39‑1 (MQ=255)                                       ggAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGaa                                      <  1:3061664/57‑1 (MQ=255)                                       ggAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGaa                                      <  1:603805/57‑1 (MQ=255)                                       ggAAGAAACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGa                                       <  1:3433909/56‑1 (MQ=255)                                        gaagaaACCAGACCGAACTCTTTGGTCGCGTTTGGTGa                                                        <  1:1453570/38‑1 (MQ=255)                                        gaagaaACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaaaa                        >  1:1556777/1‑70 (MQ=255)                                          agaaACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaaaa                        >  1:2099270/1‑68 (MQ=255)                                          agaaACCAGACCGAACTCTTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaaaa                        >  1:3378674/1‑68 (MQ=255)                                           gaaACCAGACCGAACTCTTTGGTCGCGTTTGGTGaa                                                       <  1:1501204/36‑1 (MQ=255)                                            aaaCCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATAc                     <  1:907777/69‑1 (MQ=255)                                              aCCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGaa                                      >  1:1030322/1‑50 (MQ=255)                                              aCCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGa                                    >  1:432230/1‑52 (MQ=255)                                              aCCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGa                                    >  1:552628/1‑52 (MQ=255)                                              aCCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGa                                    >  1:170834/1‑52 (MQ=255)                                              aCCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACcgac                 >  1:3622830/1‑71 (MQ=255)                                              aCCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACcga                  >  1:4505927/1‑70 (MQ=255)                                              aCCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACcg                   >  1:129958/1‑69 (MQ=255)                                              aCCAGACCGAACTCTTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACcgac                 >  1:1100696/1‑71 (MQ=255)                                               ccAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTgg                                        <  1:2005249/47‑1 (MQ=255)                                               ccAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACcgacg                <  1:2782886/71‑1 (MQ=255)                                                cAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATg                                  >  1:3711029/1‑52 (MQ=255)                                                 aGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGtcatca                           >  1:4575626/1‑58 (MQ=255)                                                 aGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGAt              <  1:2370131/71‑1 (MQ=255)                                                 aGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGAt              <  1:806457/71‑1 (MQ=255)                                                   aCCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaa                          <  1:929443/57‑1 (MQ=255)                                                   aCCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaa                          <  1:1433161/57‑1 (MQ=255)                                                   aCCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaa                          <  1:204692/57‑1 (MQ=255)                                                     cGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGa                                    >  1:3355736/1‑45 (MQ=255)                                                     cGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaaaa                        >  1:202650/1‑57 (MQ=255)                                                     cGAACTCTTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAATc          <  1:1533502/71‑1 (MQ=255)                                                       aaCTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaa                          <  1:3409823/53‑1 (MQ=255)                                                       aaCTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAAtta         >  1:2259148/1‑68 (MQ=255)                                                        aCTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAATCAGc       >  1:3989084/1‑71 (MQ=255)                                                         ctctTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAATCAg        >  1:3769662/1‑69 (MQ=255)                                                         ctctTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAATCAGc       >  1:2372807/1‑70 (MQ=255)                                                           ctTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATaa            >  1:3264429/1‑63 (MQ=255)                                                           ctTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATaa            >  1:4082589/1‑63 (MQ=255)                                                           ctTTAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATaa            >  1:2505803/1‑63 (MQ=255)                                                            tttGGTCGCGTTTGGTTAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAATCAGCCCCg   >  1:2776868/1‑71 (MQ=255)                                                            tttGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGaa                                      <  1:166113/36‑1 (MQ=255)                                                            tttGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaa                          <  1:4334217/48‑1 (MQ=255)                                                            tttGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCaa                          <  1:4511072/48‑1 (MQ=255)                                                            tttGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACcgacga               >  1:4711325/1‑59 (MQ=255)                                                            tttGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATaa            >  1:2787861/1‑62 (MQ=255)                                                            tttGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATaa            <  1:3219163/62‑1 (MQ=255)                                                            tttAGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGa                                    >  1:4293732/1‑38 (MQ=255)                                                              tGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATg                                  >  1:4769479/1‑38 (MQ=255)                                                              tGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAATCAGCCCCGc  >  1:1762391/1‑70 (MQ=255)                                                              tGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAATCAGCCCCGc  <  1:4055571/70‑1 (MQ=255)                                                               |                                                                     GTGTAAAGATATGGCACCAGTGTAAAGATGACCGAGACGGAAGAAACCAGACCGAACTCTTTGGTCGCGTTTGGTGAAATGCTGCTGAGCTGGAAGATGGTCATCAAAATACCGACGATAATCAGCCCCGC  >  W3110S.gb/4340609‑4340739

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