Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,857,513	C→T	21.0%	G58E (GGA→GAA)	ydjG ←	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,857,513	0	C	T	21.0%	169.3 / 35.8	76	G58E (GGA→GAA)	ydjG	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (24/36);  new base T (6/10);  total (30/46)
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

CAACCTGTTCACGGGGCAGTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCC  >  W3110S.gb/1857445‑1857582                                                                     |                                                                      cAACCTGTTCACGGGGCAGTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATcc                                                                      <  1:4508513/70‑1 (MQ=255) cAACCTGTTCACGGGGCAGTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATcc                                                                      <  1:3644098/70‑1 (MQ=255)     cTGTTCACGGGGCAGTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGgcg                                                                 >  1:5033742/1‑71 (MQ=255)     cTGTTCACGGGGCAGTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGgcg                                                                 <  1:262052/71‑1 (MQ=255)        ttCACGGGGCAGTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGgcgc                                                                <  1:877925/69‑1 (MQ=255)              gggCAGTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTatcaa                                                        <  1:1945309/71‑1 (MQ=255)               ggCAGTGTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATTCTGgcg                                                                 <  1:2489485/61‑1 (MQ=255)                   gTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATTCTGgcgc                                                                <  1:692713/58‑1 (MQ=255)                   gTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATTCTGgcg                                                                 <  1:4698246/57‑1 (MQ=255)                   gTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTatcaat                                                       >  1:1303265/1‑67 (MQ=255)                     ttttttAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTatcaat                                                       >  1:3231146/1‑65 (MQ=255)                     ttttttAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAt                                                  >  1:1804675/1‑70 (MQ=255)                     ttttttAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAt                                                  >  1:1934632/1‑70 (MQ=255)                      tttttGACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 >  1:570064/1‑70 (MQ=255)                      tttttAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATTCTGGCGCAGTATCAATCAGAtt                                                 >  1:1449941/1‑70 (MQ=255)                      tttttAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATTCTGGCGCAGTATCAATCAGAtt                                                 >  1:952062/1‑70 (MQ=255)                      tttttAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 <  1:4633101/70‑1 (MQ=255)                      tttttAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 >  1:3610150/1‑70 (MQ=255)                      tttttAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 >  1:2393859/1‑70 (MQ=255)                      tttttAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 >  1:2387902/1‑70 (MQ=255)                              gCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCATTATCAATCAGAtt                                                 <  1:4530661/62‑1 (MQ=255)                              gCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 <  1:3641588/62‑1 (MQ=255)                              gCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 <  1:3905124/62‑1 (MQ=255)                               ccTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCa                                                               >  1:96935/1‑47 (MQ=255)                               ccTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 <  1:4189397/61‑1 (MQ=255)                                cTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTaa                                               >  1:2040538/1‑62 (MQ=255)                                cTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAAc                                      <  1:4780629/71‑1 (MQ=255)                                cTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAAc                                      <  1:365526/71‑1 (MQ=255)                                cTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAAc                                      <  1:1107302/71‑1 (MQ=255)                                cTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAAc                                      <  1:3482321/71‑1 (MQ=255)                                cTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAAc                                      <  1:316565/71‑1 (MQ=255)                                cTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAAc                                      <  1:65156/71‑1 (MQ=255)                                  gACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAg                                                              >  1:1974663/1‑45 (MQ=255)                                  gACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAg                                                              >  1:828931/1‑45 (MQ=255)                                   aCCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 <  1:3359349/57‑1 (MQ=255)                                   aCCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACg                                     <  1:1424215/69‑1 (MQ=255)                                     cgacgaTAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAg                                                    <  1:3031798/52‑1 (MQ=255)                                     cgacgaTAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATgg                                 >  1:4052292/1‑71 (MQ=255)                                     cgacgaTAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATgg                                 >  1:2950785/1‑71 (MQ=255)                                     cgacgaTAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATgg                                 >  1:5120035/1‑71 (MQ=255)                                       acgaTAACTTCACTATTGCCAAAGTTATATTCTGGCGCAGTATCAATCAGATTAATGCCACAACGa                                    >  1:5054551/1‑66 (MQ=255)                                         gaTAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCtt                             <  1:3635509/71‑1 (MQ=255)                                            aaCTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAAt                                              >  1:3741432/1‑51 (MQ=255)                                             aCTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAg                         <  1:3814319/71‑1 (MQ=255)                                                 cACTATTGCCAAAGTTATATTCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGaa                       <  1:329145/69‑1 (MQ=255)                                                 cACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGaa                       <  1:2594422/69‑1 (MQ=255)                                                  aCTATTGCCAAAGTTATATTCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATc                     <  1:3642513/70‑1 (MQ=255)                                                  aCTATTGCCAAAGTTATATTCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATc                     <  1:3253352/70‑1 (MQ=255)                                                   cTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCg                    >  1:3379290/1‑70 (MQ=255)                                                     aTTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGa                                                   <  1:498939/37‑1 (MQ=255)                                                      ttGCCAAAGTTATATTCTGGCGCAGTATCAATCAGAtt                                                 <  1:3049618/38‑1 (MQ=255)                                                      ttGCCAAAGTTATATTCTGGCGCAGTATCAATCAGAtt                                                 <  1:2800738/38‑1 (MQ=255)                                                      ttGCCAAAGTTATATCCTGGCGCAGTATCAATCAGAtt                                                 <  1:2115471/38‑1 (MQ=255)                                                      ttGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCcaca                                        <  1:638521/47‑1 (MQ=255)                                                      ttGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCcaca                                        <  1:230551/47‑1 (MQ=255)                                                      ttGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCcaca                                        <  1:4750183/47‑1 (MQ=255)                                                      ttGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGa                                    <  1:3550987/51‑1 (MQ=255)                                                      ttGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGa                                    <  1:4750273/51‑1 (MQ=255)                                                      ttGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGa                                    <  1:876789/51‑1 (MQ=255)                                                      ttGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGa                                    <  1:155547/51‑1 (MQ=255)                                                      ttGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGa                                    <  1:2541716/51‑1 (MQ=255)                                                       tGCCAAAGTTATATTCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTAt                 <  1:5087835/69‑1 (MQ=255)                                                        gCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCaa              >  1:1928739/1‑71 (MQ=255)                                                        gCCAAAGTTATATCCTGGCGAAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCaa              >  1:3708042/1‑71 (MQ=255)                                                         ccAAAGTTATATTCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCaa              >  1:4274158/1‑70 (MQ=255)                                                         ccAAAGTTATATTCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCaa              >  1:4841260/1‑70 (MQ=255)                                                         ccAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAAt                      <  1:112025/62‑1 (MQ=255)                                                         ccAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAAt                      <  1:1800371/62‑1 (MQ=255)                                                         ccAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAAt                      <  1:2779180/62‑1 (MQ=255)                                                          cAAAGTTATATCCTGGCGCAGTATCAATCAGATTaa                                               >  1:2579185/1‑36 (MQ=255)                                                          cAAAGTTATATCCTGGCGCAGTATCAATCAGATTaa                                               >  1:3435006/1‑36 (MQ=255)                                                           aaaGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGAt                                   >  1:4031179/1‑47 (MQ=255)                                                              gTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCg                    <  1:3483726/59‑1 (MQ=255)                                                               ttATATTCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACatat       <  1:160220/71‑1 (MQ=255)                                                                tataTTCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACata        >  1:2286605/1‑69 (MQ=255)                                                                 atatCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATAttt     >  1:940743/1‑71 (MQ=255)                                                                  tatCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATAttt     <  1:84951/70‑1 (MQ=255)                                                                    tCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGcc  <  1:1125330/71‑1 (MQ=255)                                                                     |                                                                      CAACCTGTTCACGGGGCAGTTTTTTTAACGCCTGACCGACGATAACTTCACTATTGCCAAAGTTATATCCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCC  >  W3110S.gb/1857445‑1857582

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