Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	67,858	G→A	13.9%	I160I (ATC→ATT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	67,858	0	G	A	13.9%	106.0 / 8.0	43	I160I (ATC→ATT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (21/16);  new base A (2/4);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCA  >  W3110S.gb/67792‑67919                                                                   |                                                              gCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCTCCTGACGCATCCAGGAGCCGATAc                                                           >  1:1938918/1‑71 (MQ=255) gCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATAc                                                           >  1:2198467/1‑71 (MQ=255) gCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCAATa                                                            <  1:1250136/70‑1 (MQ=255) gCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCAATa                                                            <  1:153482/70‑1 (MQ=255) gCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCAATa                                                            <  1:947759/70‑1 (MQ=255) gCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCAATa                                                            <  1:820583/70‑1 (MQ=255)      aTCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCATGACGCATCCAGGAGCCGATACGCTCa                                                      <  1:365214/71‑1 (MQ=255)             gACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCtt                                               <  1:61879/71‑1 (MQ=255)              aCTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTc                                                       <  1:2006343/62‑1 (MQ=255)              aCTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCa                                                      >  1:3024564/1‑63 (MQ=255)                  tCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACg                                                          >  1:3162402/1‑55 (MQ=255)                   cAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTAt                                         >  1:2621381/1‑71 (MQ=255)                       tGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTg                                     <  1:3363579/71‑1 (MQ=255)                       tGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTg                                     >  1:1011636/1‑71 (MQ=255)                        gACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGc                                    >  1:1696737/1‑71 (MQ=255)                           gggTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAg                                 <  1:1251286/71‑1 (MQ=255)                              tATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGa                              <  1:1226208/71‑1 (MQ=255)                                tCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGa                                                              <  1:71884/37‑1 (MQ=255)                                tCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGa                                                              <  1:3257177/37‑1 (MQ=255)                                  cTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTa                                          <  1:1049339/55‑1 (MQ=255)                                  cTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTa                                          <  1:2006342/55‑1 (MQ=255)                                  cTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCgg                          >  1:477171/1‑71 (MQ=255)                                    gTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGt                         >  1:1970349/1‑70 (MQ=255)                                     tttAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGt                         <  1:1626799/69‑1 (MQ=255)                                           gaCCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACgg                 <  1:3084213/71‑1 (MQ=255)                                              cGCCTGACGCATCCAGGAGCCAATACGCTCATGTGCTTGTTTATCCTGCCAGTGACCGGTAACAACGGCAt              >  1:1816948/1‑71 (MQ=255)                                              cGCCTGACGCATCCAGGAGCCAATACGCTCATGTGCTTGTTTATCCTGCCAGTGACCGGTAACAACGGCAt              >  1:2553575/1‑71 (MQ=255)                                                 cTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATTCTGCCAGTGACCGGTAACCACGGCATGtt           >  1:3014726/1‑71 (MQ=255)                                                 cTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCa                    >  1:1558262/1‑62 (MQ=255)                                                 cTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCa                    >  1:3862301/1‑62 (MQ=255)                                                 cTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGtt           >  1:349425/1‑71 (MQ=255)                                                 cTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGtt           >  1:2260120/1‑71 (MQ=255)                                                 cTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGtt           >  1:2871853/1‑71 (MQ=255)                                                 cTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGtt           >  1:262239/1‑71 (MQ=255)                                                 cTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGtt           >  1:942693/1‑71 (MQ=255)                                                   gACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCa               >  1:468444/1‑65 (MQ=255)                                                   gACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGc         >  1:256348/1‑71 (MQ=255)                                                     cGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTg       <  1:1860712/71‑1 (MQ=255)                                                        aTCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACg    <  1:719739/71‑1 (MQ=255)                                                         tCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGt                         <  1:860266/49‑1 (MQ=255)                                                           cAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCa                    >  1:2394142/1‑52 (MQ=255)                                                           cAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGa      >  1:3631735/1‑66 (MQ=255)                                                           cAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCa  <  1:3685715/70‑1 (MQ=255)                                                                   |                                                              GCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCA  >  W3110S.gb/67792‑67919

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