Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,632,686	A→G	34.8%	V14A (GTT→GCT)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,632,686	0	A	G	34.8%	91.6 / 95.0	92	V14A (GTT→GCT)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base G (18/14);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

TTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAA  >  W3110S.gb/2632618‑2632753                                                                     |                                                                    ttttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGcg                                                                   <  1:1395514/71‑1 (MQ=255)   ttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGcgtc                                                                 >  1:1670367/1‑71 (MQ=255)   ttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGcgtc                                                                 >  1:1994383/1‑71 (MQ=255)      gTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtcgt                                                               >  1:3465019/1‑70 (MQ=255)       tCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCa                                                             <  1:2629684/71‑1 (MQ=255)         aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaa                                                            <  1:762318/70‑1 (MQ=255)         aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaa                                                            <  1:3552684/70‑1 (MQ=255)         aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaa                                                            <  1:4532446/70‑1 (MQ=255)         aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaa                                                            <  1:3058789/70‑1 (MQ=255)          gCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGcg                                                                   <  1:5126857/62‑1 (MQ=255)            tGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGt                                                        >  1:4959155/1‑71 (MQ=255)             gggTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACg                                                         <  1:102420/69‑1 (MQ=255)             gggTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACg                                                         <  1:1458116/69‑1 (MQ=255)             gggTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACg                                                         <  1:2919940/69‑1 (MQ=255)               gTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCag                                                     <  1:529149/71‑1 (MQ=255)               gTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTc                                                       >  1:220131/1‑69 (MQ=255)               gTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTc                                                       >  1:2573007/1‑69 (MQ=255)               gTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTc                                                       >  1:1435796/1‑69 (MQ=255)                     aGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCaaa                                                           <  1:3537125/59‑1 (MQ=255)                     aGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCaaa                                                           <  1:2836078/59‑1 (MQ=255)                     aGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCaaa                                                           <  1:2828580/59‑1 (MQ=255)                     aGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttc                                               >  1:1026620/1‑71 (MQ=255)                     aGGTCAGCAGTATTCCGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttc                                               >  1:3370334/1‑71 (MQ=255)                           gcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGc                                         <  1:1289539/71‑1 (MQ=255)                            cagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGc                                         <  1:40101/70‑1 (MQ=255)                            cagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGc                                         >  1:3615167/1‑70 (MQ=255)                            cagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGc                                         >  1:3486862/1‑70 (MQ=255)                             agTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGa                                       <  1:1604179/71‑1 (MQ=255)                             agTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGa                                       <  1:3519494/71‑1 (MQ=255)                             agTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGa                                       <  1:4613683/71‑1 (MQ=255)                              gTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTa                                           <  1:541334/66‑1 (MQ=255)                               tATTCGGCAGAACGGTAGAGTTAGCAGGAACGAGGAGAACTTCGTCa                                                             <  1:1795506/47‑1 (MQ=39)                                aTTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATa                                     >  1:4800864/1‑70 (MQ=255)                                aTTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATAc                                    >  1:3996270/1‑71 (MQ=255)                                aTTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATAc                                    >  1:4292767/1‑71 (MQ=255)                                aTTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATAc                                    >  1:4922713/1‑71 (MQ=255)                                 ttCTGCAGAACGGTAGAGTGAGCAGGAACTAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   >  1:447339/1‑71 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   >  1:3787758/1‑71 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCttt                                            >  1:2520246/1‑62 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   <  1:199284/71‑1 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   <  1:4261907/71‑1 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   <  1:4556608/71‑1 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   >  1:4561503/1‑71 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   >  1:1218131/1‑71 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   <  1:489211/71‑1 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   <  1:2259548/71‑1 (MQ=255)                                 ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   >  1:4935618/1‑71 (MQ=255)                                   cGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaaa                                                           >  1:1394587/1‑45 (MQ=255)                                    ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                                <  1:4228152/71‑1 (MQ=255)                                    ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                                <  1:4414506/71‑1 (MQ=255)                                    ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                                <  1:2227918/71‑1 (MQ=255)                                    ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                                <  1:4215029/71‑1 (MQ=255)                                    ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                                <  1:2989170/71‑1 (MQ=255)                                    ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                                <  1:4902383/71‑1 (MQ=255)                                    ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                                <  1:3157769/71‑1 (MQ=255)                                        gAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCaaa                                                           <  1:1634648/40‑1 (MQ=255)                                          aCGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCa                                                      >  1:334423/1‑43 (MQ=255)                                           cGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATggg                          >  1:3285778/1‑70 (MQ=255)                                             gTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                                <  1:3998243/62‑1 (MQ=255)                                              tAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAat                      <  1:679177/71‑1 (MQ=255)                                              tAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAat                      <  1:4095527/71‑1 (MQ=255)                                                gagTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAatat                    >  1:847622/1‑71 (MQ=255)                                                 agTGAGCAGGAACGAGGAGAGGGTCGTCAAACGTCa                                                      >  1:3320918/1‑36 (MQ=38)                                                 agTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCa                                                      >  1:3700992/1‑36 (MQ=255)                                                 agTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCa                                                      >  1:4416887/1‑36 (MQ=255)                                                  gTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATAtct                  >  1:2991524/1‑71 (MQ=255)                                                  gTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATAtct                  >  1:1375843/1‑71 (MQ=255)                                                  gTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAata                     >  1:2634164/1‑68 (MQ=255)                                                  gTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAata                     >  1:2276535/1‑68 (MQ=255)                                                  gTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAata                     >  1:4056313/1‑68 (MQ=255)                                                     aGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGa               >  1:957455/1‑71 (MQ=255)                                                       cAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGAc              >  1:3501350/1‑70 (MQ=255)                                                        aGGAACGAGGAGAGCGTCGTCAGACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCa            <  1:272087/71‑1 (MQ=255)                                                         tgAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCag           <  1:4128634/70‑1 (MQ=255)                                                         ggAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCag           <  1:884237/71‑1 (MQ=255)                                                          gAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCaga          <  1:3605238/71‑1 (MQ=255)                                                          gAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCaga          <  1:2076751/71‑1 (MQ=255)                                                           aaCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATAtctc                 >  1:851440/1‑63 (MQ=255)                                                            aCGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAata                     >  1:925578/1‑58 (MQ=255)                                                            aCGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGt        >  1:3279987/1‑71 (MQ=255)                                                            aCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATTGGCAATATCTCGACCAGAGt        >  1:1654120/1‑71 (MQ=255)                                                             cGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                                   <  1:4499770/43‑1 (MQ=255)                                                             cGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTg       <  1:2621418/71‑1 (MQ=255)                                                              gaggagAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTgg      >  1:1923933/1‑71 (MQ=255)                                                              gaggagAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTgg      >  1:1917934/1‑71 (MQ=255)                                                               aggagAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCa                              >  1:4627886/1‑46 (MQ=255)                                                                ggagAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCa            >  1:3067964/1‑63 (MQ=255)                                                                 gagaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCa                              >  1:290136/1‑44 (MQ=255)                                                                 gagaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCa                              >  1:5141434/1‑44 (MQ=255)                                                                 gagaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCa                              >  1:1120010/1‑44 (MQ=255)                                                                   gaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTtaa  >  1:1639320/1‑70 (MQ=255)                                                                   gaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTtaa  >  1:4336721/1‑70 (MQ=255)                                                                     |                                                                    TTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAA  >  W3110S.gb/2632618‑2632753

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