Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,892	T→C	12.0%	K454K (AAA→AAG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,892	0	T	C	12.0%	194.8 / 15.7	75	K454K (AAA→AAG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (38/28);  new base C (2/7);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 7.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

CCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAA  >  minE/608822‑608960                                                                       |                                                                     ccATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGt                                                                      >  1:3921489/1‑71 (MQ=255)  cATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGt                                                                      <  1:3660560/70‑1 (MQ=255)  cATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGt                                                                      <  1:2087947/70‑1 (MQ=255)   aTCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                    <  1:750591/71‑1 (MQ=255)    tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                                   <  1:2916003/71‑1 (MQ=255)    tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                                   >  1:3446003/1‑71 (MQ=255)    tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                                   <  1:1473380/71‑1 (MQ=255)      aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                 <  1:38364/71‑1 (MQ=255)      aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                 <  1:2151196/71‑1 (MQ=255)         tcatcaTAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACAAACCTGCATTTTCAGtt                                                                     >  1:2104412/1‑64 (MQ=255)         tcatcaTAGTTCAGGACATCGCCTTTGATAGGTTCAGACTTCGGACCAACCTGCATTTTCAGtt                                                                     >  1:1984154/1‑64 (MQ=255)          catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                              >  1:2737071/1‑70 (MQ=255)           atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGt                                                               <  1:119104/68‑1 (MQ=255)           atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                             <  1:1936985/70‑1 (MQ=255)           atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                             <  1:2478310/70‑1 (MQ=255)              atAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                         >  1:1906577/1‑71 (MQ=255)                  ttCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                    <  1:3587387/56‑1 (MQ=255)                  ttCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                    <  1:2348685/56‑1 (MQ=255)                   tCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                    <  1:392670/55‑1 (MQ=255)                   tCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                    <  1:2254112/55‑1 (MQ=255)                      ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                 >  1:2616954/1‑55 (MQ=255)                      ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                  >  1:3316256/1‑70 (MQ=255)                         cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                              <  1:3944075/71‑1 (MQ=255)                         cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                              >  1:266883/1‑71 (MQ=255)                         cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                              <  1:3893848/71‑1 (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)                                  tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGt                                                               >  1:3864525/1‑45 (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:941882/1‑71 (MQ=255)                                           cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                            >  1:479043/1‑71 (MQ=255)                                           cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                            >  1:743499/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:2579454/1‑71 (MQ=255)                                                  cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                     >  1:2574866/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:1858069/1‑51 (MQ=255)                                                                  ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                         >  1:3053695/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:569910/1‑45 (MQ=255)                                                                   tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                              >  1:3214221/1‑45 (MQ=255)                                                                    cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                       >  1:659043/1‑35 (MQ=255)                                                                    cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                       >  1:2858225/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:3825802/1‑71 (MQ=255)                                                                     aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa  >  1:922123/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)                                                                       |                                                                     CCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAA  >  minE/608822‑608960

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