Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,085,271	A→G	68.2%	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	68.2%	57.1 / 35.8	44	S127G (AGC→GGC)	sbcB	exonuclease I
Reads supporting (aligned to +/- strand):  ref base A (9/5);  new base G (4/26);  total (13/31)
Fisher's exact test for biased strand distribution p-value = 1.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

TCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTAT  >  W3110S.gb/2085212‑2085327                                                            |                                                         tCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGgcagc                                                 >  1:3025555/1‑69 (MQ=255) tCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGgcagc                                                 >  1:2710424/1‑69 (MQ=255) tCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGgcagc                                                 >  1:3453043/1‑69 (MQ=255)    acgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTAGGCCTGGGGCTGGCAGCAt                                               <  1:3478990/68‑1 (MQ=255)    acgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCAt                                               <  1:875016/68‑1 (MQ=255)    acgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCAt                                               <  1:2151362/68‑1 (MQ=255)      gacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGAt                                            <  1:1298921/69‑1 (MQ=255)       acgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGAt                                            <  1:3068709/68‑1 (MQ=255)         gaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAAc                                         >  1:122807/1‑69 (MQ=255)         gaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAAc                                         >  1:2304407/1‑69 (MQ=255)         gaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAAc                                         >  1:641648/1‑69 (MQ=255)                  cGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTgg                                <  1:984599/69‑1 (MQ=255)                  cGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTgg                                <  1:2123880/69‑1 (MQ=255)                  cGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTgg                                <  1:3660357/69‑1 (MQ=255)                      aCATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGAtt                            <  1:123944/69‑1 (MQ=255)                        aTTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCAt                                               <  1:2349282/48‑1 (MQ=255)                        aTTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCAt                                               <  1:326333/48‑1 (MQ=255)                        aTTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCAt                                               <  1:375993/48‑1 (MQ=255)                        aTTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCAt                                               <  1:752420/48‑1 (MQ=255)                        aTTTTTTATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCAt                                               <  1:931929/48‑1 (MQ=255)                         ttttttATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTAc                         >  1:1980989/1‑69 (MQ=255)                         ttttttATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTAc                         >  1:3400641/1‑69 (MQ=255)                         ttttttATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTAc                         >  1:1052541/1‑69 (MQ=255)                          tttttATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGAtt                            <  1:990057/65‑1 (MQ=255)                          tttttATCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGAtt                            <  1:1556789/65‑1 (MQ=255)                           ttttATCGTAATTTCTACGATACTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGAtt                            <  1:3238067/64‑1 (MQ=255)                               aTCGTAATTTCTACGATGCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGAt                    <  1:1399068/68‑1 (MQ=255)                                tCGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGAt                    <  1:1228388/67‑1 (MQ=255)                                 cGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGt                  <  1:1016565/68‑1 (MQ=255)                                 cGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGt                  <  1:20811/68‑1 (MQ=255)                                 cGTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGt                  <  1:1990545/68‑1 (MQ=255)                                  gTAATTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGt                  <  1:623416/67‑1 (MQ=255)                                    aaTTTCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTAt               <  1:2076353/68‑1 (MQ=255)                                       ttCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGt           <  1:1788999/69‑1 (MQ=255)                                       ttCTACGATCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGt           <  1:2407347/69‑1 (MQ=255)                                       ttCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGt           >  1:97345/1‑69 (MQ=255)                                              aTCCTTACGCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGt     <  1:1407904/68‑1 (MQ=255)                                              aTCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTTGGATTTACTGGATGTTATGCGTGCCTGt     <  1:278597/68‑1 (MQ=255)                                              aTCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGt     <  1:205050/68‑1 (MQ=255)                                              aTCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGt     <  1:1077094/68‑1 (MQ=255)                                                ccTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTAt  <  1:3560135/69‑1 (MQ=255)                                                ccTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTAt  >  1:2404252/1‑69 (MQ=255)                                                      gCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGt     >  1:2676825/1‑60 (MQ=255)                                                      gCCTGGGGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGt     >  1:1669385/1‑60 (MQ=255)                                                            |                                                         TCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTAT  >  W3110S.gb/2085212‑2085327

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