Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,883	C→A	10.8%	M457I (ATG→ATT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,883	0	C	A	10.8%	234.3 / 9.7	74	M457I (ATG→ATT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base A (2/6);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAC  >  minE/608815‑608953                                                                     |                                                                       aTGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAt                                                                      >  1:1594974/1‑71 (MQ=255)  tGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtt                                                                     >  1:3445172/1‑71 (MQ=255)  tGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtt                                                                     <  1:2717356/71‑1 (MQ=255)    cgcTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtttt                                                                   <  1:1198758/71‑1 (MQ=255)       tCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAg                                                                >  1:87381/1‑71 (MQ=255)        ccATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGt                                                               >  1:3921489/1‑71 (MQ=255)         cATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGt                                                               <  1:2087947/70‑1 (MQ=255)         cATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGt                                                               <  1:3660560/70‑1 (MQ=255)          aTCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                             <  1:750591/71‑1 (MQ=255)           tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtttt                                                                   >  1:3144332/1‑64 (MQ=255)           tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                            <  1:2916003/71‑1 (MQ=255)           tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                            <  1:1473380/71‑1 (MQ=255)           tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                            >  1:3446003/1‑71 (MQ=255)             aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCa                                                                 <  1:2900239/64‑1 (MQ=255)             aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                          <  1:38364/71‑1 (MQ=255)             aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                          <  1:2151196/71‑1 (MQ=255)              cTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGc                                                                        <  1:2597445/56‑1 (MQ=255)              cTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGc                                                                        <  1:1440972/56‑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:2348685/56‑1 (MQ=255)                         ttCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                             <  1:3587387/56‑1 (MQ=255)                          tCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                             <  1:2254112/55‑1 (MQ=255)                          tCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                             <  1:392670/55‑1 (MQ=255)                             ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                          >  1:2616954/1‑55 (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:3283948/1‑70 (MQ=255)                                    gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                    >  1:3545287/1‑70 (MQ=255)                                    gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                    >  1:685629/1‑70 (MQ=255)                                      cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                 <  1:2993798/71‑1 (MQ=255)                                         tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGt                                                        >  1:3864525/1‑45 (MQ=255)                                         tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                               >  1:3720852/1‑70 (MQ=255)                                         tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                               >  1:401538/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:1660998/1‑71 (MQ=255)                                                         cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt              >  1:2065075/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)                                                                     |                                                                       ATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAC  >  minE/608815‑608953

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