Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,580,786	G→A	7.2%	C216C (TGC→TGT)	ypfG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,580,786	0	G	A	7.2%	207.0 / 3.9	69	C216C (TGC→TGT)	ypfG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (2/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGA  >  W3110S.gb/2580718‑2580850                                                                     |                                                                 atcatcatcAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGag                                                                  <  1:3716262/69‑1 (MQ=255)  tcatcatcaGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGag                                                                  <  1:3496515/68‑1 (MQ=255)    atcatcaGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAgcg                                                             >  1:428378/1‑71 (MQ=255)       atcaGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAgcg                                                             >  1:3028831/1‑68 (MQ=255)       atcaGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCa                                                           <  1:2550042/70‑1 (MQ=255)       atcaGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCa                                                           <  1:738792/70‑1 (MQ=255)        tcaGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAg                                                          <  1:604269/70‑1 (MQ=255)            cgcTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCAt                                                     >  1:1569593/1‑71 (MQ=255)            cgcTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCAt                                                     <  1:2110126/71‑1 (MQ=255)                ttATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCa                                                  <  1:2016139/70‑1 (MQ=255)                 tATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCa                                                      >  1:2274189/1‑65 (MQ=255)                 tATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCa                                                      >  1:2851003/1‑65 (MQ=255)                 tATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCa                                                      >  1:2773608/1‑65 (MQ=255)                 tATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCa                                                      >  1:2682296/1‑65 (MQ=255)                 tATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCa                                                  <  1:1585121/69‑1 (MQ=255)                 tATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATc                                                >  1:2672749/1‑71 (MQ=255)                 tATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATc                                                >  1:1984104/1‑71 (MQ=255)                     atcaGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGcc                                            <  1:2305219/71‑1 (MQ=255)                     atcaGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGc                                             <  1:3847137/70‑1 (MQ=255)                       cggtcagCGCAGTGACATTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAg                                          >  1:249505/3‑71 (MQ=255)                       cggtcagCGCAGTGACATTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAg                                          >  1:3171396/3‑71 (MQ=255)                        agtcagCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGt                                         >  1:1145971/1‑71 (MQ=255)                        agtcagCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGt                                         >  1:2541170/1‑71 (MQ=255)                        agtcagCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGt                                         >  1:3265521/1‑71 (MQ=255)                            agCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTc                                      >  1:2720494/1‑70 (MQ=255)                              cgcAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCa                                    >  1:744315/1‑70 (MQ=255)                               gcAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGAcc                                                       >  1:3031091/1‑50 (MQ=255)                                cAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTc                                                   <  1:1932723/53‑1 (MQ=255)                                cAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTc                                                   <  1:384118/53‑1 (MQ=255)                                  gTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATc                               <  1:1533967/71‑1 (MQ=255)                                     aCATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAg                            >  1:704792/1‑71 (MQ=255)                                       aTTCACATCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCa                          <  1:2896136/71‑1 (MQ=255)                                         tCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCa                                                      >  1:2302941/1‑41 (MQ=255)                                         tCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCa                                                      >  1:3541584/1‑41 (MQ=255)                                            ccTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCg                     <  1:3187915/71‑1 (MQ=255)                                            ccTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCg                     <  1:1148251/71‑1 (MQ=255)                                            ccTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCaaa                        <  1:2956584/68‑1 (MQ=255)                                            ccTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCaaa                        <  1:1049669/68‑1 (MQ=255)                                            ccTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCaaa                        <  1:444698/68‑1 (MQ=255)                                            ccTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCaaa                        <  1:2180877/68‑1 (MQ=255)                                             cTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCg                     <  1:3175516/70‑1 (MQ=255)                                              taacgacgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATc                      >  1:2001764/3‑68 (MQ=255)                                              tCACGACGCAATGGATCAAGCGAGCAGCGCAGACCAt                                                     <  1:1281671/37‑1 (MQ=255)                                              tCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGtt                                        <  1:444001/50‑1 (MQ=255)                                              tCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGtt                                        <  1:2325091/50‑1 (MQ=255)                                              tCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGtt                                        <  1:1990421/50‑1 (MQ=255)                                               cACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCAt                                                 >  1:771501/1‑40 (MQ=255)                                               cACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATc                               >  1:2632876/1‑58 (MQ=255)                                                acgacgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCg                                              >  1:3568807/1‑42 (MQ=255)                                                acgacgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAg                            >  1:884937/1‑60 (MQ=255)                                                acgacgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGtt                   <  1:710860/69‑1 (MQ=255)                                                acgacgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGtt                   >  1:815199/1‑69 (MQ=255)                                                aagacgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCAGCCAGTTTCCATAATCCAGCAAATCGtt                   >  1:1081997/3‑69 (MQ=255)                                                 cgacgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCg                <  1:3647472/71‑1 (MQ=255)                                                  gacgCAATGGATCAAGCGAGCAGCGCAGACCATTc                                                   >  1:2500826/1‑35 (MQ=255)                                                  gacgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGt               <  1:1551457/71‑1 (MQ=255)                                                    cgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGttc             >  1:2942855/1‑71 (MQ=255)                                                    cgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGttc             >  1:2846295/1‑71 (MQ=255)                                                       aaTGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGttcttc          <  1:69005/71‑1 (MQ=255)                                                       aaTGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGttcttc          >  1:228366/1‑71 (MQ=255)                                                       aaTGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGttcttc          >  1:2066451/1‑71 (MQ=255)                                                        aTGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCg         <  1:3753470/71‑1 (MQ=255)                                                         tGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGttcttc          >  1:2813790/1‑69 (MQ=255)                                                         tGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGttcttc          >  1:1736538/1‑69 (MQ=255)                                                             tCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATc                               <  1:2671152/44‑1 (MQ=255)                                                             tCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTga    <  1:1824000/71‑1 (MQ=255)                                                               aaGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTgaga  <  1:3071405/71‑1 (MQ=255)                                                               aaGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTgaga  <  1:2527072/71‑1 (MQ=255)                                                               aaGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTgaga  <  1:3381636/71‑1 (MQ=255)                                                               aaGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTgag   <  1:287423/70‑1 (MQ=255)                                                               aaGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTgag   <  1:3580693/70‑1 (MQ=255)                                                                     |                                                                 ATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGA  >  W3110S.gb/2580718‑2580850

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