Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,103	G→C	12.7%	L717L (CTC→CTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,103	0	G	C	12.7%	192.5 / 15.1	71	L717L (CTC→CTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (37/25);  new base C (2/7);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGG  >  minE/608034‑608172                                                                      |                                                                      gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                                                       <  1:2574956/70‑1 (MQ=255)  aaaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCCTCGag                                                                     <  1:589931/71‑1 (MQ=255)   aaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                                                       >  1:1936877/1‑68 (MQ=255)    aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                                     >  1:566372/1‑69 (MQ=255)    aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                                     >  1:529629/1‑69 (MQ=255)      gTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                                                                   >  1:3640/1‑69 (MQ=255)        aCTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAt                                                               >  1:3378394/1‑71 (MQ=255)         cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                                                              <  1:1265563/69‑1 (MQ=12)         cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                                                              <  1:3645971/69‑1 (MQ=12)           gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                             >  1:2362056/1‑70 (MQ=255)            cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                                                           >  1:2963865/1‑71 (MQ=255)                gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                                       <  1:3744660/71‑1 (MQ=255)                  cAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga                                                        <  1:2729242/68‑2 (MQ=9)                   agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCg                                                         <  1:263141/66‑1 (MQ=9)                   aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                                    <  1:385533/71‑1 (MQ=255)                   aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                                    >  1:396873/1‑71 (MQ=255)                    gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                                   <  1:803100/68‑1 (MQ=9)                    gATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                             <  1:303918/61‑1 (MQ=255)                     aTACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                   >  1:3380607/1‑70 (MQ=255)                       acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                >  1:1149639/1‑71 (MQ=255)                        cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                <  1:251882/70‑1 (MQ=255)                         acGGATGGTCAGCTGCGGAAATTTTTCCGGGTTTTCCATCGCGTCGagcagc                                                                 >  1:25087/1‑52 (MQ=255)                           ggATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttg                                            <  1:1812077/71‑2 (MQ=255)                           ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                             >  1:2612834/1‑70 (MQ=255)                                gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                             <  1:731714/49‑1 (MQ=255)                                 tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                          <  1:160098/51‑1 (MQ=255)                                 tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                      <  1:928150/71‑1 (MQ=255)                                  cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                       >  1:254837/1‑69 (MQ=255)                                   aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                    >  1:1010417/1‑71 (MQ=255)                                   aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                    >  1:3155460/1‑71 (MQ=255)                                   aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                    <  1:3241491/71‑1 (MQ=255)                                    gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                  >  1:2916952/1‑56 (MQ=255)                                    gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                     <  1:1912623/69‑1 (MQ=255)                                     cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                  >  1:40621/1‑71 (MQ=255)                                     cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                  <  1:2538351/71‑1 (MQ=255)                                     cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                  <  1:3583132/71‑1 (MQ=255)                                       taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                >  1:1934025/3‑70 (MQ=255)                                       gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                                >  1:3690568/1‑71 (MQ=255)                                         ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                          >  1:2967272/1‑43 (MQ=255)                                          gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                             <  1:1929460/71‑1 (MQ=255)                                          gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                             >  1:616967/1‑71 (MQ=12)                                             atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGt                                              <  1:1328372/51‑1 (MQ=255)                                             atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                          >  1:564393/1‑71 (MQ=255)                                              tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                         <  1:2149654/71‑1 (MQ=255)                                               ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                        >  1:534188/1‑71 (MQ=255)                                               ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                        >  1:3251582/1‑71 (MQ=255)                                               ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                        >  1:855116/1‑71 (MQ=255)                                                tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                        <  1:3278536/70‑1 (MQ=255)                                                    cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                   <  1:1278827/71‑1 (MQ=255)                                                    cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                   <  1:868636/71‑1 (MQ=255)                                                     gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                   <  1:558239/38‑1 (MQ=255)                                                     gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat                  <  1:1627687/71‑1 (MQ=255)                                                       gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                 >  1:1553947/1‑70 (MQ=255)                                                        ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                    <  1:1953432/50‑1 (MQ=255)                                                        ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                    <  1:237535/50‑1 (MQ=255)                                                        ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga               >  1:1659397/1‑71 (MQ=255)                                                          ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat              >  1:1977977/1‑70 (MQ=255)                                                           tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc            >  1:3689004/1‑71 (MQ=255)                                                            ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTGAGGTGCTGACCACCTTCGATGGATGc            <  1:3419300/70‑1 (MQ=255)                                                            ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc            <  1:1213134/70‑1 (MQ=255)                                                             cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAgg                                   <  1:904983/46‑1 (MQ=255)                                                             cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCtt          >  1:3171065/1‑71 (MQ=255)                                                              aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                      >  1:1568133/1‑58 (MQ=255)                                                              aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                      >  1:1577473/1‑58 (MQ=255)                                                              aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc         >  1:350886/1‑71 (MQ=255)                                                               tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg        >  1:136118/1‑71 (MQ=255)                                                                cgcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                    >  1:1356558/1‑42 (MQ=255)                                                                 gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                 >  1:1872647/1‑44 (MQ=255)                                                                 gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                 >  1:3904503/1‑44 (MQ=255)                                                                 gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                 >  1:2517574/1‑44 (MQ=255)                                                                   gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                             >  1:3597212/1‑46 (MQ=255)                                                                   gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt    >  1:2158204/1‑71 (MQ=255)                                                                    tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                   <  1:2830291/55‑1 (MQ=255)                                                                      gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                               >  1:3094586/1‑41 (MQ=255)                                                                      gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc         <  1:37952/63‑1 (MQ=255)                                                                      gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg  <  1:1361846/70‑1 (MQ=255)                                                                      gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg  >  1:1551415/1‑70 (MQ=255)                                                                      gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg  <  1:1775066/70‑1 (MQ=255)                                                                      gAGCAGCAATTCACGGTTCATCACGTTAACGTTCAg                                    >  1:1227772/1‑36 (MQ=255)                                                                      |                                                                      GAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGG  >  minE/608034‑608172

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