Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	982,608	A→C	22.9%	P18P (CCA→CCC)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	982,608	0	A	C	22.9%	115.8 / 55.4	83	P18P (CCA→CCC)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base C (8/11);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGC  >  minE/982539‑982676                                                                      |                                                                     gggAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCaa                                                                     <  1:3929339/71‑1 (MQ=255)  ggAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCaaa                                                                    >  1:2219350/1‑71 (MQ=255)     aaTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCa                                                                      <  1:2349602/66‑1 (MQ=255)     aaTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCa                                                                      <  1:469843/66‑1 (MQ=255)         aGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGa                                                              <  1:3856790/70‑1 (MQ=255)           gTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCa                                                                      <  1:807170/60‑1 (MQ=255)           gTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGa                                                              >  1:757482/1‑68 (MQ=255)            tCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACgg                                                           >  1:2700050/1‑70 (MQ=255)             cAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGc                                                          >  1:144578/1‑70 (MQ=255)             cAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCa                                                         >  1:1624072/1‑71 (MQ=255)              aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                                        >  1:470296/1‑71 (MQ=255)              aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                                        >  1:3164697/1‑71 (MQ=255)               gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACg                                                       <  1:1947297/71‑1 (MQ=255)                 aTGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACAAGGACGGCACGtt                                                     <  1:3106509/71‑1 (MQ=255)                   ggATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCa                                                   >  1:736647/1‑71 (MQ=255)                   ggATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCa                                                   >  1:112889/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)                       aCTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAcc                                                                  <  1:3636997/52‑1 (MQ=255)                       aCTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                                                >  1:398446/1‑70 (MQ=255)                       aCTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                                                >  1:3766279/1‑70 (MQ=255)                        cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                               <  1:743707/70‑1 (MQ=255)                        cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCa                                              <  1:1231933/71‑1 (MQ=255)                           tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAc                                                                   <  1:456281/47‑1 (MQ=255)                           tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAc                                                                   <  1:3543655/47‑1 (MQ=255)                           tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAg                                                                <  1:3566629/50‑1 (MQ=255)                           tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAg                                                                <  1:1211974/50‑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)                                cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCa                                                                      >  1:3682554/1‑39 (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:293877/71‑1 (MQ=255)                                      aacaacGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCaa                                <  1:1500612/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:3536988/1‑70 (MQ=255)                                            gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCa                           >  1:3532453/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:386150/1‑70 (MQ=255)                                                   tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGtt                    >  1:311942/1‑70 (MQ=255)                                                     gtgtGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                 <  1:2970778/71‑1 (MQ=255)                                                     gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                 <  1:3844028/71‑1 (MQ=255)                                                     gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                 <  1:238927/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:2926376/1‑63 (MQ=255)                                                               tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaaa               >  1:1087032/1‑63 (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)                                                                    ccAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                 >  1:2051826/1‑56 (MQ=255)                                                                     cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGc  <  1:3745069/70‑1 (MQ=255)                                                                      |                                                                     GGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGC  >  minE/982539‑982676

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