Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,901	A→C	14.1%	V451V (GTT→GTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,901	0	A	C	14.1%	170.2 / 30.1	85	V451V (GTT→GTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (45/28);  new base C (3/9);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 2.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

ATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTT  >  minE/608832‑608971                                                                      |                                                                       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                                        <  1:2478310/70‑1 (MQ=255) atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                                        <  1:1936985/70‑1 (MQ=255)    atAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                                    >  1:1906577/1‑71 (MQ=255)            ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                             >  1:3316256/1‑70 (MQ=255)               cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                         <  1:3893848/71‑1 (MQ=255)               cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                         <  1:3944075/71‑1 (MQ=255)               cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                         >  1:266883/1‑71 (MQ=255)                   gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                      >  1:3545287/1‑70 (MQ=255)                   gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                      >  1:685629/1‑70 (MQ=255)                   gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                      >  1:3283948/1‑70 (MQ=255)                     cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                                   <  1:2993798/71‑1 (MQ=255)                        tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 >  1:401538/1‑70 (MQ=255)                        tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 >  1:3720852/1‑70 (MQ=255)                         gATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGc                                                <  1:3535708/70‑1 (MQ=255)                           tCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                             >  1:934303/1‑71 (MQ=255)                            cGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                            >  1:3949552/1‑71 (MQ=255)                             ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                           >  1:2842195/1‑71 (MQ=255)                             ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                           <  1:798967/71‑1 (MQ=255)                               ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                                    <  1:111933/44‑1 (MQ=255)                               ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                         <  1:3929044/71‑1 (MQ=255)                                 cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                        <  1:2184315/70‑1 (MQ=255)                                 cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                       >  1:479043/1‑71 (MQ=255)                                 cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                       >  1:743499/1‑71 (MQ=255)                                 cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                       >  1:941882/1‑71 (MQ=255)                                   gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                          <  1:2513308/50‑1 (MQ=255)                                   gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                                     <  1:2220937/71‑1 (MQ=255)                                     cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                         >  1:3858894/1‑49 (MQ=255)                                     cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                            >  1:2574848/1‑62 (MQ=255)                                       tCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                  <  1:2889250/70‑1 (MQ=255)                                        cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                >  1:2065075/1‑71 (MQ=255)                                        cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                >  1:1660998/1‑71 (MQ=255)                                        cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                >  1:2574866/1‑71 (MQ=255)                                        cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                >  1:2579454/1‑71 (MQ=255)                                         ggACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCATGtt                               <  1:3209922/71‑1 (MQ=255)                                         ggACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                <  1:1811441/70‑1 (MQ=255)                                          gACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                         >  1:3177881/1‑44 (MQ=255)                                            ccgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGc                            <  1:485992/68‑1 (MQ=255)                                              gaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgc                          <  1:2869311/70‑3 (MQ=255)                                              aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                            <  1:3071286/37‑1 (MQ=255)                                               aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                          >  1:2564407/1‑38 (MQ=255)                                               aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                 >  1:1137796/1‑63 (MQ=255)                                               aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                >  1:2764278/1‑64 (MQ=255)                                               aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                         >  1:1475953/1‑71 (MQ=255)                                               aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                  <  1:2463701/62‑1 (MQ=255)                                               aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                         <  1:1521300/71‑1 (MQ=255)                                                ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                        <  1:2996212/71‑1 (MQ=255)                                                ccTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                <  1:3059050/63‑1 (MQ=255)                                                 cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                       <  1:3665939/71‑1 (MQ=255)                                                    cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                     <  1:2835929/70‑1 (MQ=255)                                                    cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc                    >  1:1781102/1‑71 (MQ=255)                                                      ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                        <  1:833394/49‑1 (MQ=255)                                                      ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                  <  1:921355/71‑1 (MQ=255)                                                       cttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                <  1:3800219/55‑1 (MQ=255)                                                        ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                            <  1:3159143/43‑1 (MQ=255)                                                        ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                    >  1:3053695/1‑51 (MQ=255)                                                        ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                    >  1:1858069/1‑51 (MQ=255)                                                        ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                    >  1:3646875/1‑51 (MQ=255)                                                        ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCga                 <  1:3563169/70‑1 (MQ=255)                                                        ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                  >  1:3240407/1‑69 (MQ=255)                                                         tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                         >  1:3214221/1‑45 (MQ=255)                                                         tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                         >  1:569910/1‑45 (MQ=255)                                                          cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                                  >  1:2858225/1‑35 (MQ=255)                                                          cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                                  >  1:659043/1‑35 (MQ=255)                                                          cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag               <  1:675343/70‑1 (MQ=255)                                                          cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag               >  1:2215457/1‑70 (MQ=255)                                                           aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                               <  1:3144932/53‑1 (MQ=255)                                                           aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa             >  1:922123/1‑71 (MQ=255)                                                           aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa             >  1:3825802/1‑71 (MQ=255)                                                            gTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGttt                              <  1:2342304/53‑1 (MQ=255)                                                             tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                             >  1:2451095/1‑53 (MQ=255)                                                             tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                             >  1:3416380/1‑53 (MQ=255)                                                             tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                          >  1:570867/1‑56 (MQ=255)                                                              ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg          >  1:3185564/1‑71 (MQ=255)                                                               tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCa                                   >  1:1163429/1‑45 (MQ=255)                                                               tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg          >  1:3422021/1‑70 (MQ=255)                                                                ttCGTCAACGCCGCCGTTGATTGCGTACATCATGGTTTTCGcc                                    <  1:1492905/43‑1 (MQ=255)                                                                ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCga                 >  1:1259392/1‑62 (MQ=255)                                                                ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa        >  1:3789051/1‑71 (MQ=255)                                                                ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa        >  1:3152497/1‑71 (MQ=255)                                                                 tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt       >  1:311948/1‑71 (MQ=255)                                                                  cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt       <  1:973447/70‑1 (MQ=255)                                                                  cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt       <  1:1852187/70‑1 (MQ=255)                                                                   gTCCACGCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt      <  1:1468722/70‑1 (MQ=255)                                                                   gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt      <  1:2778101/70‑1 (MQ=255)                                                                    tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                     <  1:2557423/54‑1 (MQ=255)                                                                     cAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                          <  1:1149445/48‑1 (MQ=255)                                                                     cAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt   >  1:206006/1‑71 (MQ=255)                                                                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCa                                   >  1:2346010/1‑38 (MQ=255)                                                                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt  >  1:3542483/1‑71 (MQ=255)                                                                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt  >  1:2708089/1‑71 (MQ=255)                                                                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt  <  1:2912248/71‑1 (MQ=255)                                                                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt   <  1:2048483/70‑1 (MQ=255)                                                                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt   <  1:2327653/70‑1 (MQ=255)                                                                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt   <  1:3233016/70‑1 (MQ=255)                                                                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGG‑TTGCACGCGCAcc                   >  1:489912/1‑53 (MQ=255)                                                                      |                                                                       ATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTT  >  minE/608832‑608971

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