Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	982,595	A→G	18.8%	Y14C (TAC→TGC)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	982,595	0	A	G	18.8%	180.4 / 35.8	85	Y14C (TAC→TGC)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base G (8/8);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAA  >  minE/982531‑982661                                                                 |                                                                   gCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGc                                                                >  1:215630/1‑69 (MQ=255) gCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCg                                                               >  1:740402/1‑70 (MQ=255)  cGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgtg                                                             >  1:708109/1‑71 (MQ=255)   gTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgtgt                                                            <  1:3417394/71‑1 (MQ=255)       atgggAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTc                                                         <  1:1102382/68‑1 (MQ=255)         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)              aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGgtgt                                                              <  1:17921/58‑1 (MQ=255)              aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGc                                                                <  1:3906143/56‑1 (MQ=255)              aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGc                                                                <  1:2852011/56‑1 (MQ=255)              aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                                              <  1:3854722/58‑1 (MQ=255)                gAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgtgt                                                            >  1:2254021/1‑58 (MQ=255)                 aGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGttt                                                          <  1:3331692/59‑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: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)                               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:1211974/50‑1 (MQ=255)                                   tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAg                                                 <  1:3566629/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)                                       tCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTc                                                         >  1:3256351/1‑38 (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: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)                                                                 |                                                                   GCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAA  >  minE/982531‑982661

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