Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	982,620	C→A	28.5%	G22G (GGC→GGA)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	982,620	0	C	A	28.5%	164.1 / 49.3	84	G22G (GGC→GGA)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base A (11/13);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

CAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTG  >  minE/982551‑982689                                                                      |                                                                      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: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)               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)                 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: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: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:851923/71‑1 (MQ=255)                                                     tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                  <  1:3158471/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)                                                                     gCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                       <  1:1522643/50‑1 (MQ=255)                                                                     gCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctggctg  >  1:1633841/1‑71 (MQ=255)                                                                      |                                                                      CAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTG  >  minE/982551‑982689

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