Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	982,670	A→C	19.1%	N39T (AAT→ACT)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	982,670	0	A	C	19.1%	116.7 / 33.3	68	N39T (AAT→ACT)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base C (3/10);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGG  >  minE/982600‑982739                                                                       |                                                                      gtgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                                       <  1:553774/71‑1 (MQ=255)    tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                                       <  1:1761545/68‑1 (MQ=255)    tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                                                                    <  1:851923/71‑1 (MQ=255)    tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                                                                    <  1:3158471/71‑1 (MQ=255)      tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTg                                                                  >  1:2929104/1‑71 (MQ=255)       cccaaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGc                                                                 >  1:1490063/4‑71 (MQ=255)        cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGc                                                                 <  1:3745069/70‑1 (MQ=255)             cAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTg                                                                     <  1:3667758/61‑1 (MQ=255)              aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCg                                                           <  1:784787/70‑1 (MQ=255)              aGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgc                                                          <  1:52041/71‑1 (MQ=255)                 aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctgg                                                       <  1:2816653/71‑1 (MQ=255)                    gCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctggctg                                                    >  1:1633841/1‑71 (MQ=255)                     aaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg                                                    <  1:3835257/69‑1 (MQ=255)                     aaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGatggttgg                                                   <  1:2702312/70‑1 (MQ=255)                     aaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGatggttgg                                                   <  1:3140303/70‑1 (MQ=255)                      aCGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                  <  1:570782/71‑1 (MQ=255)                       cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGc                                                                 >  1:1324791/1‑55 (MQ=255)                       cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                                 <  1:1362761/71‑2 (MQ=255)                       cGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaa                                                 <  1:2415209/71‑1 (MQ=255)                       cGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                  >  1:2869985/1‑70 (MQ=255)                        gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaa                                                <  1:2718327/71‑1 (MQ=255)                        gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaa                                                <  1:2181387/71‑1 (MQ=255)                        gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                  <  1:1293730/69‑1 (MQ=255)                         ttCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaag                                               <  1:519500/71‑1 (MQ=255)                          tCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaag                                               >  1:777994/1‑70 (MQ=255)                           cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                                                       >  1:27844/1‑45 (MQ=255)                           cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                              >  1:2709370/1‑70 (MQ=255)                           cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                              >  1:2005320/1‑70 (MQ=255)                           cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                              <  1:2544152/70‑1 (MQ=255)                           cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                              >  1:3461625/1‑70 (MQ=255)                             ccGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                  <  1:2051769/64‑1 (MQ=255)                              cGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACt                                          >  1:385486/1‑71 (MQ=255)                              cGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACt                                          >  1:860743/1‑71 (MQ=255)                                cAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                                                            >  1:54342/1‑51 (MQ=255)                                    gACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGc                                                                 >  1:2420514/1‑42 (MQ=38)                                        tataCCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggtt                                                     <  1:2769810/50‑3 (MQ=255)                                        tataCCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                                 <  1:3588852/54‑2 (MQ=255)                                         ataCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                                            >  1:2049951/1‑58 (MQ=255)                                          taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgagtga                               <  1:2102060/70‑1 (MQ=255)                                           aCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgct                                                         <  1:2539005/43‑1 (MQ=255)                                             cAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATg                           <  1:3491763/71‑1 (MQ=255)                                              aaaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtga                                   >  1:3530296/1‑62 (MQ=255)                                              aaaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATgg                          <  1:815703/71‑1 (MQ=255)                                               aaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGa                         >  1:1128136/1‑71 (MQ=255)                                               aaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGa                         >  1:880873/1‑71 (MQ=255)                                               aaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGa                         >  1:2436610/1‑71 (MQ=255)                                                aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgct                                                         >  1:1036833/1‑38 (MQ=255)                                                aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGAt                        <  1:543831/71‑1 (MQ=255)                                                     aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                                            >  1:56866/1‑46 (MQ=255)                                                     aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                     <  1:1946963/53‑1 (MQ=255)                                                     aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGa                   <  1:1398796/71‑1 (MQ=255)                                                     aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCACGTGAGTGATATGGATTTCGa                   >  1:113236/1‑71 (MQ=255)                                                     aaaCATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                                 <  1:2567329/41‑2 (MQ=12)                                                     aaaCATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                                 <  1:3593249/41‑2 (MQ=12)                                                      aaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGa                   >  1:252705/1‑70 (MQ=255)                                                       aCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                                            >  1:3720713/1‑44 (MQ=255)                                                       aCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAgg                                      >  1:977316/1‑50 (MQ=255)                                                        cGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGacac                <  1:1396034/71‑1 (MQ=255)                                                           ttAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                              >  1:2227764/1‑38 (MQ=255)                                                           ttAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACaa              >  1:2633507/1‑70 (MQ=255)                                                           ttAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAAt             >  1:2297196/1‑71 (MQ=255)                                                            tAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAAt             >  1:1844145/1‑70 (MQ=255)                                                             aaaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGacac                >  1:2451702/1‑66 (MQ=255)                                                             aaaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATc            <  1:2428711/70‑1 (MQ=255)                                                             aaaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATc            <  1:3512191/70‑1 (MQ=255)                                                             aaaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATc            <  1:3255836/70‑1 (MQ=255)                                                              aaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGat                              <  1:2273018/51‑1 (MQ=255)                                                              aaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGat                              <  1:2041350/51‑1 (MQ=255)                                                              aaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCAt          >  1:1740261/1‑71 (MQ=255)                                                                  cGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGa       >  1:2472902/1‑70 (MQ=255)                                                                   gAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCg                    <  1:3840614/56‑1 (MQ=255)                                                                    aaaaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                     <  1:1951568/38‑1 (MQ=255)                                                                    aaaaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCa           >  1:3197850/1‑64 (MQ=255)                                                                     aaaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAg   >  1:2054954/1‑71 (MQ=255)                                                                     aaaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAg   >  1:1614808/1‑71 (MQ=255)                                                                      aaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAgg  >  1:3101794/1‑71 (MQ=255)                                                                       |                                                                      GTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGG  >  minE/982600‑982739

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑