Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,996,898	A→T	50.0%	M196K (ATG→AAG)	ygeQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,996,898	0	A	T	50.0%	‑3.5 / 140.9	88	M196K (ATG→AAG)	ygeQ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/14);  new base T (23/21);  total (53/35)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGACCTAAAGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATC  >  W3110S.gb/2996829‑2996959                                                                      |                                                              tAAAGACCTAAAGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCAt                                                              <  1:1126738/71‑1 (MQ=255)   aaGACCTAAAGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCa                                                            >  1:4265900/1‑71 (MQ=255)    aGACCTAAAGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCa                                                            >  1:2436874/1‑70 (MQ=255)      aCCTAAAGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACgg                                                         >  1:1456051/1‑71 (MQ=255)          aaaGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCa                                                            >  1:2844537/1‑64 (MQ=255)          aaaGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCa                                                            >  1:4640379/1‑64 (MQ=255)          aaaGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACgggg                                                       >  1:1425009/1‑69 (MQ=255)           aaGCATCTTTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAg                                                    >  1:3142374/1‑71 (MQ=255)           aaGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACggggg                                                      <  1:3992824/69‑1 (MQ=255)           aaGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAg                                                    >  1:2319528/1‑71 (MQ=255)           aaGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACggggg                                                      <  1:899617/69‑1 (MQ=255)           aaGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAg                                                    >  1:2058180/1‑71 (MQ=255)            aGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGt                                                   >  1:335149/1‑71 (MQ=255)                tcatTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACg                                                          <  1:1402317/60‑1 (MQ=255)                  atTTTTTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTg                                              >  1:477887/1‑70 (MQ=255)                  atTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACggg                                                        >  1:3614996/1‑60 (MQ=255)                  atTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCa                                                            >  1:2780674/1‑56 (MQ=255)                  atTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACgg                                                         >  1:4700358/1‑59 (MQ=255)                  atTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACgg                                                         >  1:1652363/1‑59 (MQ=255)                  atTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACgg                                                         >  1:4088570/1‑59 (MQ=255)                  atTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTg                                              >  1:3568928/1‑70 (MQ=255)                  atTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCATGCACATCTGCCATCACgg                                                         >  1:126262/1‑59 (MQ=255)                   ttttGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTc                                            <  1:4501170/71‑1 (MQ=255)                     ttGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACgggg                                                       >  1:2036220/1‑58 (MQ=255)                      tGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCaaa                                         >  1:3582266/1‑71 (MQ=255)                      tGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCaaa                                         >  1:1182554/1‑71 (MQ=255)                      tGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCaa                                          <  1:4027895/70‑1 (MQ=255)                      tGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCaa                                          <  1:4796014/70‑1 (MQ=255)                      tGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCaa                                          <  1:850592/70‑1 (MQ=255)                      tGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCaa                                          <  1:1013338/70‑1 (MQ=255)                       gTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAAt                                        >  1:911603/1‑71 (MQ=255)                           aTGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCgg                                     <  1:1122042/70‑1 (MQ=255)                            tGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACTGGGGAGTTGCTGTCAAATCGGa                                    <  1:417075/70‑1 (MQ=255)                            tGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGa                                    <  1:4433488/70‑1 (MQ=255)                            tGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGat                                   <  1:1945254/71‑1 (MQ=255)                              gTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGatat                                 >  1:724439/1‑71 (MQ=255)                                tCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCAc                                                           >  1:686685/1‑43 (MQ=255)                                tCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGatat                                 >  1:4903374/1‑69 (MQ=255)                                tCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATg                                >  1:779679/1‑70 (MQ=255)                                 ccACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGca                              >  1:1461804/1‑71 (MQ=255)                                 ccACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGca                              >  1:1457221/1‑71 (MQ=255)                                  cacaGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAAt                                        >  1:1668962/1‑60 (MQ=255)                                  cacaGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAAt                                        >  1:2491570/1‑60 (MQ=255)                                    caGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGCAGTTGCTGTCAAATCGGATATGcacat                           <  1:1667198/71‑2 (MQ=255)                                      gCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTg                                                 <  1:1395186/47‑1 (MQ=255)                                      gCTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGata                                  >  1:4114588/1‑62 (MQ=255)                                       cTGCGTCTTGCCTTCTTCCTGCACATCTGCCTTCAc                                                           <  1:3490029/36‑1 (MQ=37)                                       cTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGa                                    >  1:3943878/1‑59 (MQ=255)                                       cTGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcat                        >  1:5125615/1‑71 (MQ=255)                                        tGCGGCTTGCCTTCTTCCTGCACATCTGCCTTCACggg                                                        >  1:313628/1‑38 (MQ=255)                                         gCGTCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAg                           <  1:843155/66‑1 (MQ=255)                                         gCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGat                                   >  1:2912106/1‑58 (MQ=255)                                         gCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAg                           <  1:261399/66‑1 (MQ=255)                                         gCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcatt                       >  1:257260/1‑70 (MQ=255)                                          cGGCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcatt                       <  1:2719606/69‑1 (MQ=255)                                          cGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAg                           <  1:3850915/65‑1 (MQ=255)                                           ggCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcattcat                    <  1:2245143/71‑1 (MQ=255)                                           ggCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTg                                              >  1:1129979/1‑45 (MQ=255)                                            gattGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTg                                                 <  1:509988/39‑1 (MQ=38)                                            gCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcattcat                    <  1:1990761/70‑1 (MQ=255)                                            gCTTGCCTTCTTCCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcattcat                    <  1:4841262/70‑1 (MQ=255)                                            gCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcattc                      <  1:1201914/68‑1 (MQ=255)                                              ttGCCTTCTTCCTGCACATCTGCCATCACGGTGGAGTTGCTGTCAAATCGGATATGCACAGcattca                     >  1:820632/1‑67 (MQ=255)                                              ttGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcattca                     >  1:357917/1‑67 (MQ=255)                                              ttGCCTTCTTCCTGCACATATGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCAttt                  >  1:2920466/1‑70 (MQ=255)                                                gCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGc                <  1:5128151/70‑1 (MQ=255)                                                   ttcttcCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCgg                                     <  1:1728432/46‑1 (MQ=255)                                                   ttcttcCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGa                                    <  1:2562228/47‑1 (MQ=255)                                                   ttcttcCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGa                                    <  1:1369458/47‑1 (MQ=255)                                                   ttcttcCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcattc                      <  1:679841/61‑1 (MQ=255)                                                    tcttcCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACgg            >  1:342892/1‑70 (MQ=255)                                                    tcttcCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCa               <  1:3336905/67‑1 (MQ=255)                                                     cttcCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGca                              >  1:3426271/1‑51 (MQ=255)                                                     cttcCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcattcat                    >  1:4468704/1‑61 (MQ=255)                                                     cttcCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACggg           >  1:4552504/1‑70 (MQ=255)                                                      ttcCTGCACATCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGaa         >  1:3643555/1‑71 (MQ=255)                                                      ttcCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGcac                             >  1:675446/1‑51 (MQ=255)                                                        ccTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAAtt       >  1:1636166/1‑71 (MQ=255)                                                          tGCACATCTGCCATCACGGGGGAGTTGCTGTCATATCGGa                                    <  1:8123/40‑1 (MQ=255)                                                             acaTCTGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATtccatc  <  1:3277862/71‑1 (MQ=255)                                                             acaTCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATtccatc  <  1:1908974/71‑1 (MQ=255)                                                              caTCTGCCTTCACGGGGGAGTTGCTGTCAAATCgg                                     <  1:4927449/35‑1 (MQ=255)                                                                  tGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcattcat                    >  1:474601/1‑48 (MQ=255)                                                                  tGCCTTCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGcattcat                    >  1:4626456/1‑48 (MQ=255)                                                                    ccATCACGGGGGAGTTGCTGTCAAATCGGATATGca                              >  1:4743081/1‑36 (MQ=255)                                                                    ccATCACGGGGGAGTTGCTGTCAAATCGGATATGca                              >  1:2573666/1‑36 (MQ=255)                                                                    ccATCACGGGGGAGTTGCTGTCAAATCGGATATGca                              >  1:4603793/1‑36 (MQ=255)                                                                    ccATCACGGGGGAGTTGCTGTCAAATCGGATATGca                              >  1:2716261/1‑36 (MQ=255)                                                                      |                                                              TAAAGACCTAAAGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATC  >  W3110S.gb/2996829‑2996959

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