Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,086,423	G→A	20.1%	A66A (GCG→GCA)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,086,423	0	G	A	20.1%	101.3 / 15.2	69	A66A (GCG→GCA)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (9/5);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

GGCACAGGCGGTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGT  >  minE/2086355‑2086482                                                                     |                                                            ggCACAGGCGGTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGgcg                                                             >  1:2376729/1‑69 (MQ=255) ggCACAGGCGGTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGgcg                                                             >  1:1804245/1‑69 (MQ=255) ggCACAGGCGGTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGgcg                                                             >  1:1739270/1‑69 (MQ=255) ggCACAGGCGGTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGgcg                                                             >  1:1407574/1‑69 (MQ=255) ggCACAGGCGGTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGgcg                                                             >  1:1446367/1‑69 (MQ=255) ggCACAGGCGGTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGgcg                                                             >  1:1616172/1‑69 (MQ=255) ggCACAGGCGGTTTTATAGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGgcg                                                             >  1:252386/1‑69 (MQ=255)      aGGCGGTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGcgcaccg                                                        <  1:2129042/69‑1 (MQ=255)          ggTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACg                                                    >  1:3149895/1‑69 (MQ=255)              ttATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTg                                                 >  1:1449609/1‑68 (MQ=255)              ttATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTg                                                 >  1:655932/1‑68 (MQ=255)              ttATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTg                                                 >  1:778829/1‑68 (MQ=255)              ttATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTg                                                 >  1:901627/1‑68 (MQ=255)                aTCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGAt                                               <  1:2545044/68‑1 (MQ=255)                             aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:736083/69‑1 (MQ=255)                             aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:2356178/69‑1 (MQ=255)                             aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:2695811/69‑1 (MQ=255)                             aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:2088646/69‑1 (MQ=255)                             aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:3341617/69‑1 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGTGcgcaccgcac                                                     >  1:1657921/1‑48 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGcgcaccgca                                                      >  1:1131974/1‑47 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:2672779/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3286326/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:2307985/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3132640/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:2862125/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3380941/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:2287064/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3294415/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1136321/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:428489/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1961304/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:176701/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1629300/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:902186/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1368484/1‑69 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCa                                 >  1:911021/1‑68 (MQ=255)                              tCGTTATCATCTCTATGTTTCACTCACCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3209422/1‑69 (MQ=255)                               cGTTATAATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1129350/1‑68 (MQ=255)                                gTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAg                              <  1:205515/69‑1 (MQ=255)                                gTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAg                              <  1:2695212/69‑1 (MQ=255)                                  tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAgg                             <  1:136365/68‑1 (MQ=255)                                  tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAgg                             <  1:2826920/68‑1 (MQ=255)                                  tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGa                            >  1:3267349/1‑69 (MQ=255)                                  tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGa                            >  1:2733958/1‑69 (MQ=255)                                  tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGa                            >  1:2034952/1‑69 (MQ=255)                                     catcTCTATGTTTCACTCGCCTGCCCGTGGGcgcaccg                                                        <  1:3412092/38‑1 (MQ=255)                                               tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:1609419/67‑1 (MQ=255)                                               tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:2899456/67‑1 (MQ=255)                                               tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:672507/67‑1 (MQ=255)                                               tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:1309249/67‑1 (MQ=255)                                               tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:886676/67‑1 (MQ=255)                                               tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:3280258/67‑1 (MQ=255)                                                ttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAt              >  1:2847203/1‑69 (MQ=255)                                                ttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTCTAt              >  1:75260/1‑69 (MQ=255)                                                 tCACTCGCCTGCCCGTGGGCGCACCGTACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:2361702/69‑1 (MQ=255)                                                 tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:2420371/69‑1 (MQ=255)                                                 tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:1804032/69‑1 (MQ=255)                                                 tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:3033062/69‑1 (MQ=255)                                                 tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:895914/69‑1 (MQ=255)                                                 tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:2718161/69‑1 (MQ=255)                                                      cGCCTGCCCGTGGGCGCACCGCACGCTGATCATgcgct                                       <  1:952597/38‑2 (MQ=255)                                                      cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:340838/38‑1 (MQ=255)                                                      cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:346644/38‑1 (MQ=255)                                                      cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:1380019/38‑1 (MQ=255)                                                      cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:2009319/38‑1 (MQ=255)                                                      cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:2165364/38‑1 (MQ=255)                                                      cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:56226/38‑1 (MQ=255)                                                      cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:700801/38‑1 (MQ=255)                                                              cGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATttccgtttccgt  <  1:2264929/67‑1 (MQ=255)                                                                 gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:787209/1‑58 (MQ=255)                                                                 gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:768848/1‑58 (MQ=255)                                                                 gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:2798603/1‑58 (MQ=255)                                                                 gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:1395387/1‑58 (MQ=255)                                                                 gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:399766/1‑58 (MQ=255)                                                                 gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:2326437/1‑58 (MQ=255)                                                                 gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:295853/1‑58 (MQ=255)                                                                 gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccg         >  1:2560376/1‑57 (MQ=255)                                                                 gggCACACCGCACGCTGATCATGCGTAAGCTAAAAGGACt                          >  1:695683/1‑40 (MQ=25)                                                                  ggCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:3013778/62‑1 (MQ=255)                                                                     gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:747816/54‑1 (MQ=255)                                                                     gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:1551110/54‑1 (MQ=255)                                                                     gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:585069/54‑1 (MQ=255)                                                                     gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:845790/54‑1 (MQ=255)                                                                     gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:562192/54‑1 (MQ=255)                                                                     |                                                            GGCACAGGCGGTTTTATCGCAGAGAAAGATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGT  >  minE/2086355‑2086482

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