Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	986,248	C→T	38.6%	I163I (ATC→ATT)	ynfE →	oxidoreductase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	986,248	0	C	T	38.6%	48.6 / 56.1	88	I163I (ATC→ATT)	ynfE	oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base T (12/22);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

GTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTA  >  minE/986181‑986313                                                                    |                                                                  gTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCg                                                                  <  1:1040113/69‑1 (MQ=255) gTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCg                                                                  <  1:1653103/69‑1 (MQ=255)   aGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCg                                                                  >  1:2662966/1‑67 (MQ=255)   aGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGtt                                                                <  1:1063936/69‑1 (MQ=255)   aGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGtt                                                                <  1:1686844/69‑1 (MQ=255)   aGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGtt                                                                <  1:1168732/69‑1 (MQ=255)    gCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTg                                                               >  1:130778/1‑69 (MQ=255)    gCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTg                                                               >  1:2740792/1‑69 (MQ=255)    gCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTg                                                               >  1:2837210/1‑69 (MQ=255)                        cAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:293796/69‑1 (MQ=255)                        cAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:1273453/69‑1 (MQ=255)                        cAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:3164640/69‑1 (MQ=255)                        cAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:2182829/69‑1 (MQ=255)                        cAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:640371/69‑1 (MQ=255)                         aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:2851348/68‑1 (MQ=255)                         aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:2672283/68‑1 (MQ=255)                         aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:534732/68‑1 (MQ=255)                         aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:1427811/68‑1 (MQ=255)                            atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCg                                        >  1:2615386/1‑68 (MQ=255)                            atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGc                                       >  1:264010/1‑69 (MQ=255)                            atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGc                                       >  1:2427687/1‑69 (MQ=255)                            atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGc                                       >  1:2138223/1‑69 (MQ=255)                            atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGc                                       >  1:982644/1‑69 (MQ=255)                            atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGc                                       >  1:2059266/1‑69 (MQ=255)                             tGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:690351/64‑1 (MQ=255)                                         tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:3125454/52‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggctg                        <  1:2484088/67‑3 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1392613/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:3402034/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:386132/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:3077072/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1142105/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:2894328/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:2841499/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:607665/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1533385/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:2439650/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:700753/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1885281/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:2290297/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1905572/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1071952/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1779022/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1969549/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1897522/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:1693056/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:989995/67‑1 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcgg                        >  1:1115810/1‑67 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGt                       >  1:114489/1‑68 (MQ=255)                                             tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGt                       >  1:1349599/1‑68 (MQ=255)                                               tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCa                     <  1:588295/68‑1 (MQ=255)                                               tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCa                     <  1:3278377/68‑1 (MQ=255)                                               tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCa                     <  1:456043/68‑1 (MQ=255)                                               tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCa                     <  1:2986539/68‑1 (MQ=255)                                                 ttCAGTACTCTTCGGGGATCGTTGGCGGCAATGTGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAAc                  >  1:3338556/1‑69 (MQ=255)                                                 ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGcc                                      >  1:1036038/1‑49 (MQ=255)                                                 ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGGGGTCAAAc                  >  1:3228586/1‑69 (MQ=255)                                                 ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAAc                  >  1:2659349/1‑69 (MQ=255)                                                 ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAAc                  >  1:550656/1‑69 (MQ=255)                                                 ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAAc                  >  1:1577636/1‑69 (MQ=255)                                                 ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAAc                  >  1:3181749/1‑69 (MQ=255)                                                 ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGGCAAAc                  >  1:948631/1‑69 (MQ=255)                                                 ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGAGGTCAAAc                  >  1:1199870/1‑69 (MQ=255)                                                  tCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                           <  1:448959/43‑1 (MQ=255)                                                      tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtg             >  1:2915213/1‑69 (MQ=255)                                                      tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtg             >  1:2951468/1‑69 (MQ=255)                                                      tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtg             >  1:692686/1‑69 (MQ=255)                                                      tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCt              >  1:3087930/1‑68 (MQ=255)                                                      tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCt              >  1:2122670/1‑68 (MQ=255)                                                          ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatga         <  1:851754/69‑1 (MQ=255)                                                          ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatga         <  1:998854/69‑1 (MQ=255)                                                          ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatga         <  1:155164/69‑1 (MQ=255)                                                          ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatga         <  1:1954604/69‑1 (MQ=255)                                                          ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatga         <  1:2706781/69‑1 (MQ=255)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:2931800/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:2927977/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:12330/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:1197702/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:3246379/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:2355338/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:805392/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:2296002/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:2243434/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:2234238/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCa                                     >  1:1657924/1‑40 (MQ=38)                                                           ttCGGGGATTGTTGGCGGAAATATGACCCGCTCTTCTCCa                                     >  1:1420492/1‑40 (MQ=37)                                                            tCGGGGATTGTTGGCGGCAATATGACCCGCTCTTCTCCATCAGCCTcggcgg                        <  1:2386024/52‑1 (MQ=39)                                                                 gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTa  >  1:3114787/1‑69 (MQ=255)                                                                    |                                                                  GTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTA  >  minE/986181‑986313

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