Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,249	T→C	10.6%	A335A (GCA→GCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,249	0	T	C	10.6%	196.4 / 12.1	75	A335A (GCA→GCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base C (3/5);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACG  >  Exported/609181‑609314                                                                     |                                                                  aGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                                                  <  1:1815458/70‑1 (MQ=255) aGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGc                                                                 <  1:3223435/71‑1 (MQ=255)     aaCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAg                                                             <  1:5120824/71‑1 (MQ=255)      aCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGa                                                            <  1:1618778/71‑1 (MQ=255)       cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                                                >  1:2858751/1‑66 (MQ=255)       cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGAt                                                           <  1:1468218/71‑1 (MQ=255)        ggAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                                                >  1:4980736/1‑65 (MQ=255)           aGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                                                        <  1:2476638/70‑1 (MQ=255)            ggagTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGt                                                      <  1:2964678/68‑1 (MQ=255)                tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTggg                                                       >  1:4588029/1‑66 (MQ=255)                tttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                                   <  1:151642/70‑1 (MQ=255)                 ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                                   >  1:1515688/1‑69 (MQ=255)                 ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                                   >  1:2138134/1‑69 (MQ=255)                  tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCg                                                         >  1:2000489/1‑62 (MQ=255)                  tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                                                  >  1:2712682/1‑69 (MQ=255)                  tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                                                  >  1:2474570/1‑69 (MQ=255)                   ttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                                                <  1:1269378/70‑3 (MQ=255)                     ggTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCagag                                              >  1:1969865/1‑70 (MQ=255)                      gtcaCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGaa                                            <  1:4104932/68‑1 (MQ=255)                      gTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                                             <  1:2512797/70‑1 (MQ=255)                         aCCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                                         >  1:845085/1‑71 (MQ=255)                         aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGcc                                                  <  1:765651/62‑1 (MQ=255)                         aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAg                                         <  1:3508966/71‑1 (MQ=255)                           cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGtt                                       >  1:4221212/1‑71 (MQ=255)                           cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                        >  1:3660695/1‑70 (MQ=255)                           cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                        >  1:2647042/1‑70 (MQ=255)                           cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                                                >  1:2796722/1‑60 (MQ=255)                              ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                                    <  1:1293650/71‑1 (MQ=255)                                tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                  <  1:3889893/71‑1 (MQ=255)                                tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                  <  1:2079290/71‑1 (MQ=255)                                  cGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                                                  <  1:3613668/37‑1 (MQ=255)                                  cGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                                                  <  1:3607041/37‑1 (MQ=255)                                    aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt                              >  1:1698968/1‑71 (MQ=255)                                    aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAt                               >  1:624818/1‑70 (MQ=255)                                     ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                             >  1:572776/1‑71 (MQ=255)                                     ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                             >  1:4124624/1‑71 (MQ=255)                                       gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                                                        <  1:2518763/42‑1 (MQ=255)                                       gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                           <  1:850734/71‑1 (MQ=255)                                       gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                           <  1:1691527/71‑1 (MQ=255)                                        tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                          <  1:3307660/55‑1 (MQ=255)                                        tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                          <  1:4986670/55‑1 (MQ=255)                                        tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                          <  1:3899602/55‑1 (MQ=255)                                        tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                          <  1:4247526/55‑1 (MQ=255)                                        tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                           >  1:4205076/1‑70 (MQ=255)                                           aGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg                        >  1:4514107/1‑70 (MQ=255)                                               ccATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTGCGGAGTACg                    >  1:3484677/1‑70 (MQ=255)                                                cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                           <  1:4291754/46‑1 (MQ=255)                                                cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                           <  1:3358088/46‑1 (MQ=255)                                                cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                           <  1:2868013/46‑1 (MQ=255)                                                 aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                           <  1:4092524/45‑1 (MQ=255)                                                 aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                 >  1:4629968/1‑71 (MQ=255)                                                  tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                           >  1:2804630/1‑44 (MQ=255)                                                  tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                                         >  1:797551/1‑46 (MQ=255)                                                  tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                <  1:1712763/71‑1 (MQ=255)                                                  tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                <  1:3820930/71‑1 (MQ=255)                                                    ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa               <  1:3361070/70‑1 (MQ=255)                                                    ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa               <  1:821195/70‑1 (MQ=255)                                                      accGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                           >  1:378396/1‑56 (MQ=255)                                                       ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                             >  1:5616054/1‑53 (MQ=255)                                                        cGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                                             >  1:2719212/1‑36 (MQ=255)                                                        cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGaa                                            >  1:4196245/1‑37 (MQ=255)                                                        cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg           <  1:1657633/70‑1 (MQ=255)                                                        cGATAGATTCGGTTGCCCAGATCGGGTAGCCAGAGaa                                            >  1:2396641/1‑37 (MQ=255)                                                         gATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATc                                   >  1:503729/1‑45 (MQ=255)                                                          aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTc                                      <  1:5490046/41‑1 (MQ=255)                                                          aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc        >  1:4952182/1‑71 (MQ=255)                                                           tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                             <  1:4114301/49‑1 (MQ=255)                                                           tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa          <  1:2988351/68‑1 (MQ=255)                                                            aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGc            >  1:1460356/1‑65 (MQ=255)                                                            aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc       >  1:5282517/1‑70 (MQ=255)                                                            aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc        >  1:4756260/1‑69 (MQ=255)                                                            aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa      >  1:399609/1‑71 (MQ=255)                                                             gATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa      <  1:4702818/70‑1 (MQ=255)                                                             acttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCa      >  1:5649674/3‑70 (MQ=255)                                                                 cGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACg  >  1:4007674/1‑70 (MQ=255)                                                                     |                                                                  AGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACG  >  Exported/609181‑609314

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