Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,886	T→C	11.2%	K456K (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,886	0	T	C	11.2%	189.7 / 13.2	71	K456K (AAA→AAG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base C (2/6);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCG  >  minE/608816‑608955                                                                       |                                                                      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:3660560/70‑1 (MQ=255)        cATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGt                                                                 <  1:2087947/70‑1 (MQ=255)         aTCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                               <  1:750591/71‑1 (MQ=255)          tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtttt                                                                     >  1:3144332/1‑64 (MQ=255)          tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                              <  1:1473380/71‑1 (MQ=255)          tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                              >  1:3446003/1‑71 (MQ=255)          tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                              <  1:2916003/71‑1 (MQ=255)            aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCa                                                                   <  1:2900239/64‑1 (MQ=255)            aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                            <  1:2151196/71‑1 (MQ=255)            aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                            <  1:38364/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: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:685629/1‑70 (MQ=255)                                   gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                      >  1:3545287/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:743499/1‑71 (MQ=255)                                                 cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                       >  1:479043/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: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)                                                                       |                                                                      TGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCG  >  minE/608816‑608955

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