Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	982,605	T→C	22.3%	F17F (TTT→TTC)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	982,605	0	T	C	22.3%	159.9 / 44.3	85	F17F (TTT→TTC)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base C (8/11);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTG  >  minE/982537‑982675                                                                     |                                                                       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)           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:736647/1‑71 (MQ=255)                     ggATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCa                                                  >  1:112889/1‑71 (MQ=255)                       aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAcc                                                <  1:2350193/71‑1 (MQ=255)                         aaaaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGc                                              >  1:3303811/4‑71 (MQ=255)                         aaaaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGc                                              >  1:3427573/4‑71 (MQ=255)                         aaaaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGc                                              >  1:3550938/4‑71 (MQ=255)                         aCTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAcc                                                                 <  1:3636997/52‑1 (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)                             tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAc                                                                  <  1:3543655/47‑1 (MQ=255)                             tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAc                                                                  <  1:456281/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: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: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:1087032/1‑63 (MQ=255)                                                                 tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaaa              >  1:2926376/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)                                                                     |                                                                       AGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTG  >  minE/982537‑982675

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