Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,729,737	T→C	38.7%	V53V (GTT→GTC)	yheV →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,729,737	0	T	C	38.7%	45.3 / 62.4	93	V53V (GTT→GTC)	yheV	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base C (15/21);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GTTGAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTT  >  minE/2729675‑2729797                                                               |                                                             gttgAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCg                                                            >  1:1057366/1‑65 (MQ=255) gttgAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCg                                                            >  1:3413250/1‑65 (MQ=255) gttgAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCg                                                            >  1:2282966/1‑65 (MQ=255) gttgAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCg                                                            >  1:467592/1‑65 (MQ=255) gttgAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCg                                                            >  1:486020/1‑65 (MQ=255) gttgAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCg                                                            >  1:1342676/1‑65 (MQ=255) gttgAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCaaa                                                        <  1:2162029/69‑1 (MQ=255)    gAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAg                                                       >  1:3062688/1‑67 (MQ=255)    gAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAg                                                       >  1:294120/1‑67 (MQ=255)       tgtgTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGa                                                   >  1:2687851/1‑68 (MQ=255)         tgtTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCa                                                <  1:2180694/69‑1 (MQ=255)         tgtTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCa                                                >  1:3001739/1‑69 (MQ=255)         tgtTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCa                                                >  1:492824/1‑69 (MQ=255)         tgtTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCa                                                >  1:1750911/1‑69 (MQ=255)         tgtTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCa                                                >  1:507138/1‑69 (MQ=255)         tgtTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCa                                                >  1:2541159/1‑69 (MQ=255)         tgtTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCa                                                <  1:2650191/69‑1 (MQ=255)                         cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                >  1:272795/1‑69 (MQ=255)                         cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                >  1:879979/1‑69 (MQ=255)                         cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                >  1:2628183/1‑69 (MQ=255)                         cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                >  1:42955/1‑69 (MQ=255)                         cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAttttt                                >  1:3025297/1‑69 (MQ=255)                         cAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGAt                                    >  1:2059247/1‑65 (MQ=255)                          aGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTc                               >  1:136315/1‑69 (MQ=255)                          aGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTc                               >  1:576993/1‑69 (MQ=255)                          aGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTc                               >  1:1267291/1‑69 (MQ=255)                          aGATGAGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTc                               >  1:2803385/1‑69 (MQ=255)                              gCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATc                            >  1:2791897/1‑68 (MQ=255)                               cGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCg                          <  1:3149655/69‑1 (MQ=255)                               cGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCg                          <  1:1044581/69‑1 (MQ=255)                               cGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCg                          <  1:3389714/69‑1 (MQ=255)                               cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                          <  1:1967322/69‑1 (MQ=255)                               cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                          <  1:81187/69‑1 (MQ=255)                               cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                          <  1:208281/69‑1 (MQ=255)                               cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                          <  1:571897/69‑1 (MQ=255)                               cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                          <  1:3191993/69‑1 (MQ=255)                               cGAGAAGCAGACAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCg                          <  1:1702234/69‑1 (MQ=255)                                 agaAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTg                                            <  1:154012/49‑1 (MQ=255)                                 agaAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTg                                            <  1:3032118/49‑1 (MQ=255)                                 agaAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTg                                            <  1:2070410/49‑1 (MQ=255)                                 agaAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTg                                            <  1:2114766/49‑1 (MQ=255)                                  gaAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTg                                            <  1:1351737/48‑1 (MQ=255)                                       aGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACt                      <  1:2997150/65‑1 (MQ=255)                                        gACAAAGAAGCCCGCGATCACGTTCGCAAAGATTAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                 <  1:1676607/69‑1 (MQ=255)                                        gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                 <  1:1322341/69‑1 (MQ=255)                                        gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                 <  1:1621567/69‑1 (MQ=255)                                        gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                 <  1:2687694/69‑1 (MQ=255)                                        gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                 <  1:933207/69‑1 (MQ=255)                                        gACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGa                 <  1:2750307/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:3068546/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:3099564/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:872100/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:2818013/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:326866/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:2601553/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:1356487/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:2446828/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:2420765/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:226562/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:1081273/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:2209708/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:1725706/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGata               >  1:193958/1‑69 (MQ=255)                                          cAAAGAAGCCCGCGATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGat                >  1:416360/1‑68 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:2856202/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:1306684/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:3412933/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:1223865/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:3316423/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:1127117/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:3244751/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:909398/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:1956693/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:3148382/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:2128957/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:213839/69‑1 (MQ=255)                                                      cgATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:2747812/69‑1 (MQ=255)                                                       gATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:278390/68‑1 (MQ=255)                                                       gATCACGTCCGCAAAGATGAGCAAGTGATTGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGttt   <  1:1363913/68‑1 (MQ=255)                                                        aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt  >  1:1314574/1‑68 (MQ=255)                                                        aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt  >  1:1912987/1‑68 (MQ=255)                                                        aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt  >  1:1431896/1‑68 (MQ=255)                                                        aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt  >  1:1663940/1‑68 (MQ=255)                                                        aTCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtttt  >  1:1927134/1‑68 (MQ=255)                                                          cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt    <  1:1563535/64‑1 (MQ=255)                                                          cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt    <  1:1780539/64‑1 (MQ=255)                                                          cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt    <  1:1176238/64‑1 (MQ=255)                                                          cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt    <  1:3243164/64‑1 (MQ=255)                                                          cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt    <  1:2225266/64‑1 (MQ=255)                                                          cACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGtt    <  1:939758/64‑1 (MQ=255)                                                               |                                                             GTTGAATGTGTTAAGTGCGGACATCAGATGCGAGAAGCAGACAAAGAAGCCCGCGATCACGTTCGCAAAGATGAGCAAGTGATCGGGATTTTTCATCCGGACTAGCGATATGCGCCGAGTTTT  >  minE/2729675‑2729797

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