Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,234	A→G	19.4%	G340G (GGT→GGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,234	0	A	G	19.4%	151.3 / 28.5	72	G340G (GGT→GGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base G (8/6);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

TGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGA  >  Exported/609164‑609301                                                                       |                                                                    tGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATacca                                                                     >  1:1355781/1‑71 (MQ=255)   gtgtGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATcccgcc                                                                   >  1:3131886/1‑71 (MQ=255)   gtgtACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATaccac                                                                    >  1:1664215/1‑70 (MQ=255)      tACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGAt                                                                <  1:2213120/71‑1 (MQ=255)        cAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGa                                                                 >  1:1558318/1‑68 (MQ=255)        cAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATa                                                               <  1:4373535/70‑1 (MQ=255)          gggTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATAACACCGATAGAt                                                            <  1:4219930/71‑1 (MQ=255)               ttCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGt                                                       >  1:146097/1‑71 (MQ=255)                  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)                        cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCg                                                                  >  1:254371/1‑51 (MQ=255)                        cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGa                                                                 >  1:1655566/1‑52 (MQ=255)                        cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGa                                                                 >  1:3229867/1‑52 (MQ=255)                        cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                                   >  1:2858751/1‑66 (MQ=255)                        cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGAt                                              <  1:1468218/71‑1 (MQ=255)                         ggAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAt                                                                >  1:2928662/1‑52 (MQ=16)                         ggAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAt                                                                >  1:3848782/1‑52 (MQ=16)                         ggAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                                   >  1:4980736/1‑65 (MQ=255)                            aGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                                           <  1:2476638/70‑1 (MQ=255)                             ggagTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGt                                         <  1:2964678/68‑1 (MQ=255)                                gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                              >  1:2007762/1‑45 (MQ=16)                                gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                              >  1:1750265/1‑45 (MQ=16)                                gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                              >  1:5579429/1‑45 (MQ=16)                                gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                              >  1:417144/1‑45 (MQ=16)                                gTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCgg                                                        >  1:4654558/1‑53 (MQ=255)                                 tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtc                                                               <  1:3875433/45‑2 (MQ=12)                                 tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTggg                                          >  1:4588029/1‑66 (MQ=255)                                 tttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                      <  1:151642/70‑1 (MQ=255)                                  ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                      >  1:2138134/1‑69 (MQ=255)                                  ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                      >  1:1515688/1‑69 (MQ=255)                                   tttGGTAACCAGGGTACGACCGTCGAGGCCCATacca                                                                     >  1:1993855/1‑37 (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:2647042/1‑70 (MQ=255)                                            cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                           >  1:3660695/1‑70 (MQ=255)                                            cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                                   >  1:2796722/1‑60 (MQ=255)                                               ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                       <  1:1293650/71‑1 (MQ=255)                                                 tACGACCGTCGAGGCCCATACCACCGATAGATTCgg                                                        >  1:4678409/1‑36 (MQ=255)                                                 tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                     <  1:2079290/71‑1 (MQ=255)                                                 tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                     <  1:3889893/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:4124624/1‑71 (MQ=255)                                                      ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                >  1:572776/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:4986670/55‑1 (MQ=255)                                                         tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                             <  1:3307660/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:3358088/46‑1 (MQ=255)                                                                 cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                              <  1:4291754/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:3820930/71‑1 (MQ=255)                                                                   tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg   <  1:1712763/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)                                                                       |                                                                    TGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGA  >  Exported/609164‑609301

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