Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	982,617	C→T	27.5%	D21D (GAC→GAT)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	982,617	0	C	T	27.5%	130.6 / 57.1	84	D21D (GAC→GAT)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base T (11/12);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGG  >  minE/982550‑982686                                                                    |                                                                      tCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACgg                                                                     >  1:2700050/1‑70 (MQ=255)  cAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGc                                                                    >  1:144578/1‑70 (MQ=255)  cAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCa                                                                   >  1:1624072/1‑71 (MQ=255)   aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                                                  >  1:3164697/1‑71 (MQ=255)   aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                                                  >  1:470296/1‑71 (MQ=255)    gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACg                                                                 <  1:1947297/71‑1 (MQ=255)      aTGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACAAGGACGGCACGtt                                                               <  1:3106509/71‑1 (MQ=255)        ggATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCa                                                             >  1:112889/1‑71 (MQ=255)        ggATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCa                                                             >  1:736647/1‑71 (MQ=255)          aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAcc                                                           <  1:2350193/71‑1 (MQ=255)            aaaaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGc                                                         >  1:3550938/4‑71 (MQ=255)            aaaaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGc                                                         >  1:3427573/4‑71 (MQ=255)            aaaaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGc                                                         >  1:3303811/4‑71 (MQ=255)            aCTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                                                          >  1:3766279/1‑70 (MQ=255)            aCTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                                                          >  1:398446/1‑70 (MQ=255)             cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                         <  1:743707/70‑1 (MQ=255)             cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCa                                                        <  1:1231933/71‑1 (MQ=255)                tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAcc                                                           <  1:932822/65‑1 (MQ=255)                tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCa                                                        <  1:2412675/68‑1 (MQ=255)                tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                      >  1:3322947/1‑70 (MQ=255)                 cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATg                                                     <  1:52514/70‑1 (MQ=255)                  gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAc                                                                       >  1:1808230/1‑51 (MQ=255)                  gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                       >  1:2157967/1‑67 (MQ=255)                  gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                   >  1:3867343/1‑71 (MQ=255)                   aTCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                  >  1:2787413/1‑71 (MQ=255)                     ctcggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGACATGACGta                                                >  1:1891871/4‑71 (MQ=255)                     cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                >  1:2786558/1‑71 (MQ=255)                        ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                  <  1:1468420/66‑1 (MQ=255)                         gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                            <  1:1824892/71‑1 (MQ=255)                           aacaacGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCaa                                          <  1:1500612/71‑1 (MQ=255)                           aacaacGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCaa                                          <  1:293877/71‑1 (MQ=255)                            acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                         <  1:2874776/55‑1 (MQ=255)                            acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                          >  1:226656/1‑70 (MQ=255)                            acaacGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGACGGAACGTTCACTGCCATGACGTATACCaaa                                         <  1:1574190/71‑1 (MQ=255)                             caacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                          >  1:2300351/1‑69 (MQ=255)                             caacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAAGTTCa                                                             >  1:2322837/1‑50 (MQ=255)                              aacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                                                          >  1:2506737/1‑52 (MQ=255)                               acGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATg                                                     <  1:3672106/56‑1 (MQ=255)                               acGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGc                                      <  1:2943475/71‑1 (MQ=255)                               acGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAAGGCACGTTCACCGCAATGACGTATACCAAAAGc                                      <  1:499196/71‑1 (MQ=255)                               acGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGc                                      <  1:2543305/71‑1 (MQ=255)                                cGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCGCTGCCATGACGTATACCAAAAGCa                                     <  1:744226/71‑1 (MQ=255)                                cGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCa                                           >  1:2377970/1‑65 (MQ=255)                                cGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAg                                       >  1:2063079/1‑69 (MQ=255)                                cGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGc                                      >  1:667690/1‑70 (MQ=255)                                cGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCa                                     <  1:296415/71‑1 (MQ=255)                                 gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                    <  1:3690503/71‑1 (MQ=255)                                 gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCa                                     >  1:3532453/1‑70 (MQ=255)                                 gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCa                                     >  1:3536988/1‑70 (MQ=255)                                   aTCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGc                                      <  1:39852/67‑1 (MQ=255)                                      tctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGt                               >  1:1360865/1‑71 (MQ=255)                                       ctGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACAtt                              >  1:1404636/1‑71 (MQ=255)                                       ctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                  >  1:816828/1‑67 (MQ=255)                                       ctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                                   >  1:2360745/1‑66 (MQ=255)                                        tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGtt                              >  1:311942/1‑70 (MQ=255)                                        tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGtt                              >  1:386150/1‑70 (MQ=255)                                          gtgtGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                           <  1:2970778/71‑1 (MQ=255)                                          gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                           <  1:238927/71‑1 (MQ=255)                                          gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                           <  1:3844028/71‑1 (MQ=255)                                             tACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCa                                     <  1:176364/58‑1 (MQ=255)                                               cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                            >  1:3916346/1‑65 (MQ=255)                                               cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAAcc                       <  1:3676827/70‑1 (MQ=255)                                               cGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                           <  1:2735242/66‑1 (MQ=255)                                                ggCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGa                     <  1:41107/71‑1 (MQ=255)                                                 gCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGa                     >  1:2761011/1‑70 (MQ=255)                                                 gCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACAttt                             <  1:1738321/62‑1 (MQ=255)                                                  cGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaa                          <  1:23565/64‑1 (MQ=255)                                                   gtgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                  <  1:553774/71‑1 (MQ=255)                                                    tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaaa                         >  1:1087032/1‑63 (MQ=255)                                                    tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaaa                         >  1:2926376/1‑63 (MQ=255)                                                      tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                  <  1:1761545/68‑1 (MQ=255)                                                      tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg               <  1:3158471/71‑1 (MQ=255)                                                      tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg               <  1:851923/71‑1 (MQ=255)                                                        tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTg             >  1:2929104/1‑71 (MQ=255)                                                         cccaaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGc            >  1:1490063/4‑71 (MQ=255)                                                         ccAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                           >  1:2051826/1‑56 (MQ=255)                                                          cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGc            <  1:3745069/70‑1 (MQ=255)                                                           caaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGaaa                   >  1:3023933/2‑62 (MQ=255)                                                               cAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTg                <  1:3667758/61‑1 (MQ=255)                                                                aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCg      <  1:784787/70‑1 (MQ=255)                                                                aGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgc     <  1:52041/71‑1 (MQ=255)                                                                 ggATGGAACGTTCACTGGCATGACGTATACCAAAAGCAAAACATTTaaaa                         >  1:2531996/1‑50 (MQ=255)                                                                  gACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                            <  1:572479/46‑1 (MQ=255)                                                                   aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGttt                             <  1:671834/44‑1 (MQ=255)                                                                   aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctgg  <  1:2816653/71‑1 (MQ=255)                                                                    |                                                                      TCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGG  >  minE/982550‑982686

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