Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	982,587	T→C	17.1%	S11S (TCT→TCC)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	982,587	0	T	C	17.1%	185.4 / 30.1	82	S11S (TCT→TCC)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base C (7/7);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

TTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGT  >  minE/982518‑982657                                                                      |                                                                       ttCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAAtctc                                                                       <  1:3111044/71‑1 (MQ=255) ttCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAAtctc                                                                       <  1:2498114/71‑1 (MQ=255)  tCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATctct                                                                      >  1:1720894/1‑71 (MQ=255)   ccGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTg                                                                     >  1:606498/1‑71 (MQ=255)    cGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTg                                                                     <  1:3857060/70‑1 (MQ=255)     gCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTgg                                                                    >  1:1395317/1‑70 (MQ=255)         ccccAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGgtgt                                                                 <  1:3888447/69‑1 (MQ=255)              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)                  tAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGgtg                                                                  <  1:760299/59‑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:2852011/56‑1 (MQ=255)                           aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGc                                                            <  1:3906143/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)                                 tCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGgtg                                                                  <  1:1890546/44‑1 (MQ=255)                                 tCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGgtg                                                                  <  1:366003/44‑1 (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:3164697/1‑71 (MQ=255)                                   aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                     >  1:470296/1‑71 (MQ=255)                                    gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACg                                    <  1:1947297/71‑1 (MQ=255)                                     cATGGATACTATCGATCTTGGCAACAACGAAtctc                                                                       >  1:828380/1‑35 (MQ=255)                                     cATGGATACTATCGATCTTGGCAACAACGAAtctc                                                                       >  1:2980205/1‑35 (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: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: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)                                                                      |                                                                       TTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGT  >  minE/982518‑982657

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