Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,085,271	A→G	46.1%	S127G (AGC→GGC)	sbcB →	exonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,085,271	0	A	G	46.1%	4.8 / 51.4	39	S127G (AGC→GGC)	sbcB	exonuclease I
Reads supporting (aligned to +/- strand):  ref base A (4/17);  new base G (12/6);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 3.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCC  >  W3110S.gb/2085212‑2085331                                                            |                                                             tCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCAt                                                   >  1:2537418/1‑71 (MQ=255)   gacgacgaAGTCACCCGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGa                                                 >  1:2285729/1‑71 (MQ=255)   gacgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGa                                                 >  1:84198/1‑71 (MQ=255)   gacgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGa                                                 >  1:2468581/1‑71 (MQ=255)    acgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCt                                                           <  1:2950449/60‑1 (MQ=255)    acgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCt                                                           <  1:2871361/60‑1 (MQ=255)    acgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCt                                                           <  1:3280488/60‑1 (MQ=255)    acgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGAt                                                <  1:3715632/71‑1 (MQ=255)      gacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATaa                                              >  1:283800/1‑71 (MQ=255)                    cAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGAtt                                <  1:2370691/71‑1 (MQ=255)                    cAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGAtt                                <  1:3648586/71‑1 (MQ=255)                    cAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGAtt                                <  1:2473687/71‑1 (MQ=255)                    cAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGAtt                                <  1:2366522/71‑1 (MQ=255)                    cAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGAtt                                <  1:2218903/71‑1 (MQ=255)                    cAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGAtt                                <  1:1931606/71‑1 (MQ=255)                    cAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGAtt                                <  1:63337/71‑1 (MQ=255)                    cAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGAtt                                <  1:1407410/71‑1 (MQ=255)                           ttttATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGa                         >  1:1431841/1‑71 (MQ=255)                           ttttATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGa                         >  1:1070056/1‑71 (MQ=255)                            tttATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGAt                        >  1:2364919/1‑71 (MQ=255)                               aTCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGt                      <  1:2022408/70‑1 (MQ=255)                               aTCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGt                      <  1:2764301/70‑1 (MQ=255)                               aTCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGt                      <  1:762062/70‑1 (MQ=255)                                  gTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATg                  >  1:994424/1‑71 (MQ=255)                                  gTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATg                  >  1:3294029/1‑71 (MQ=255)                                  gTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATg                  >  1:860302/1‑71 (MQ=255)                                  gTAATTTCTAAGATCCTTACGCCTGGGGCTGGCAGCATGATAACTcgc                                         >  1:822532/1‑48 (MQ=38)                                   tAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCAt                                                   <  1:2781643/37‑1 (MQ=255)                                   tAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCAt                                                   <  1:94055/37‑1 (MQ=255)                                   tAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCAt                                                   <  1:260215/37‑1 (MQ=255)                                        tCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTa                    >  1:2885914/1‑63 (MQ=255)                                        tCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTa                    >  1:58962/1‑63 (MQ=255)                                        tCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTa                    >  1:1305352/1‑63 (MQ=255)                                        tCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAAc                                             >  1:1522535/1‑38 (MQ=255)                                            cGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGt         <  1:1127213/70‑1 (MQ=255)                                            cGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGt         <  1:2703777/70‑1 (MQ=255)                                             gATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTa       <  1:2699685/71‑1 (MQ=255)                                             gATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTa       <  1:267415/71‑1 (MQ=255)                                                    aCGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGccc  <  1:440173/69‑1 (MQ=255)                                                            |                                                             TCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCC  >  W3110S.gb/2085212‑2085331

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