Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,729,758	C→T	49.4%	I60I (ATC→ATT)	yheV →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,729,758	0	C	T	49.4%	‑0.3 / 76.5	85	I60I (ATC→ATT)	yheV	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/21);  new base T (15/27);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

CAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGCTGCAGAA  >  minE/2729699‑2729826                                                            |                                                                     cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                                             >  1:2628183/1‑69 (MQ=255) cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                                             >  1:272795/1‑69 (MQ=255) cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                                             >  1:879979/1‑69 (MQ=255) cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                                             >  1:3025297/1‑69 (MQ=255) cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                                             >  1:42955/1‑69 (MQ=255) cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAt                                                                 >  1:2059247/1‑65 (MQ=255)  aGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTc                                                            >  1:1267291/1‑69 (MQ=255)  aGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTc                                                            >  1:576993/1‑69 (MQ=255)  aGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTc                                                            >  1:136315/1‑69 (MQ=255)  aGATGAGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTc                                                            >  1:2803385/1‑69 (MQ=255)      gCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATc                                                         >  1:2791897/1‑68 (MQ=255)       cGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCg                                                       <  1:1044581/69‑1 (MQ=255)       cGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCg                                                       <  1:3149655/69‑1 (MQ=255)       cGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCg                                                       <  1:3389714/69‑1 (MQ=255)       cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                                                       <  1:3191993/69‑1 (MQ=255)       cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                                                       <  1:1702234/69‑1 (MQ=255)       cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                                                       <  1:208281/69‑1 (MQ=255)       cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                                                       <  1:571897/69‑1 (MQ=255)       cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                                                       <  1:1967322/69‑1 (MQ=255)       cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                                                       <  1:81187/69‑1 (MQ=255)               aGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACt                                                   <  1:2997150/65‑1 (MQ=255)                gACAAAGAAGCCCGCGATCACGTTCGCAAAGATTAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                                              <  1:1676607/69‑1 (MQ=255)                gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                                              <  1:933207/69‑1 (MQ=255)                gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                                              <  1:1621567/69‑1 (MQ=255)                gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                                              <  1:2750307/69‑1 (MQ=255)                gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                                              <  1:1322341/69‑1 (MQ=255)                gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                                              <  1:2687694/69‑1 (MQ=255)                 aCAAAGAAGCCCTCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                                              <  1:3376710/68‑1 (MQ=255)                 aCAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                                              <  1:2482214/68‑1 (MQ=255)                 aCAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                                              <  1:2688978/68‑1 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:1081273/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:2601553/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:2446828/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:872100/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:2818013/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:3068546/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:3099564/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:326866/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:2420765/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:226562/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:1725706/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:1356487/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:2209708/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata                                            >  1:193958/1‑69 (MQ=255)                  cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGat                                             >  1:416360/1‑68 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:3244751/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:3316423/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:909398/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:3148382/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:1127117/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:3412933/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:2856202/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:1956693/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:2128957/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:2747812/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:213839/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:1306684/69‑1 (MQ=255)                              cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:1223865/69‑1 (MQ=255)                               gATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:278390/68‑1 (MQ=255)                               gATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:1363913/68‑1 (MQ=255)                                aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt                               >  1:1927134/1‑68 (MQ=255)                                aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt                               >  1:1912987/1‑68 (MQ=255)                                aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt                               >  1:1431896/1‑68 (MQ=255)                                aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt                               >  1:1314574/1‑68 (MQ=255)                                aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt                               >  1:1663940/1‑68 (MQ=255)                                  cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt                                 <  1:1176238/64‑1 (MQ=255)                                  cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt                                 <  1:1563535/64‑1 (MQ=255)                                  cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt                                 <  1:3243164/64‑1 (MQ=255)                                  cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt                                 <  1:1780539/64‑1 (MQ=255)                                  cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt                                 <  1:939758/64‑1 (MQ=255)                                  cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt                                 <  1:2225266/64‑1 (MQ=255)                                       ccGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:2533631/59‑1 (MQ=255)                                        cGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTa                      <  1:3337602/69‑1 (MQ=255)                                        cGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTa                      <  1:3433038/69‑1 (MQ=255)                                        cGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTa                      >  1:610988/1‑69 (MQ=255)                                            aaGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt                                <  1:2905734/55‑1 (MQ=255)                                            aaGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTg                   >  1:733514/1‑68 (MQ=255)                                                 gAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTaca             >  1:1880111/1‑69 (MQ=255)                                                 gAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTaca             >  1:2261647/1‑69 (MQ=255)                                                     aaGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGc         <  1:1174256/69‑1 (MQ=255)                                                     aaGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGc         <  1:2433036/69‑1 (MQ=255)                                                     aaGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGc         <  1:94159/69‑1 (MQ=255)                                                     aaGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGc         <  1:3011296/69‑1 (MQ=255)                                                        tGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGc         >  1:478848/1‑66 (MQ=255)                                                        tGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGc         >  1:3233310/1‑66 (MQ=255)                                                            cGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGCTGCAGaa  >  1:2032430/1‑69 (MQ=255)                                                            cGGGATATTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGCTGCAGaa  >  1:2454/1‑69 (MQ=255)                                                            |                                                                     CAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTTTTTAAGCTAGTGAGTACACGGCTGCAGAA  >  minE/2729699‑2729826

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