Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,372	C→T	25.0%	A628T (GCG→ACG) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,372	0	C	T	25.0%	88.8 / 24.6	52	A628T (GCG→ACG) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (22/17);  new base T (5/8);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAG  >  Exported/608302‑608441                                                                       |                                                                      gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgc                                                                       <  1:3688099/71‑1 (MQ=255)   aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                                                                     >  1:980134/1‑71 (MQ=255)   aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                                                                     >  1:3105067/1‑71 (MQ=255)         agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                               <  1:2490337/71‑1 (MQ=255)         agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                               <  1:1399283/71‑1 (MQ=255)            ggCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                             >  1:2161674/1‑70 (MQ=255)             gCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacgacg                                                           <  1:3342709/71‑1 (MQ=255)                 cggCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCg                                                       >  1:616796/3‑71 (MQ=255)                  ggCACCTTTGAGGTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCg                                                       <  1:276826/69‑1 (MQ=255)                   gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                             >  1:3639095/1‑63 (MQ=255)                   gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                             >  1:3379603/1‑63 (MQ=255)                   gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACgacgac                                                         >  1:933964/1‑67 (MQ=255)                    cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACgacgac                                                         >  1:5403555/1‑66 (MQ=255)                    cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                                    >  1:5213759/1‑71 (MQ=255)                    cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                                    >  1:5464753/1‑71 (MQ=255)                      ccTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCttg                                                  <  1:2552236/71‑3 (MQ=255)                         ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg                                                                  >  1:2966184/1‑52 (MQ=255)                         ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                               >  1:4138635/1‑71 (MQ=255)                         ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                               >  1:2646935/1‑71 (MQ=255)                          tgcggTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCgt                                               <  1:1831786/67‑1 (MQ=255)                          tCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                               <  1:2347125/70‑1 (MQ=255)                           cTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                             >  1:2409864/1‑55 (MQ=255)                              gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgtgt                                             <  1:502969/68‑1 (MQ=255)                              gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                                          >  1:5106630/1‑71 (MQ=255)                              gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                                          >  1:3947601/1‑71 (MQ=255)                              gTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAc                                          <  1:1035538/71‑1 (MQ=255)                               tCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCgta                                              >  1:4677054/1‑65 (MQ=255)                                 aCGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTggg                                                 >  1:928333/1‑61 (MQ=255)                                 aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccgg                                       >  1:4037640/1‑69 (MQ=255)                                 aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccg                                        <  1:2079948/70‑2 (MQ=255)                                   gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                                     <  1:1515024/71‑1 (MQ=255)                                   gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                                     <  1:3841163/71‑1 (MQ=255)                                    aCCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGttt                                    <  1:4754005/71‑1 (MQ=255)                                       gTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCt                                 <  1:2869956/71‑1 (MQ=255)                                        tGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                                     >  1:5672208/1‑66 (MQ=255)                                         gCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTa                               >  1:4342200/1‑71 (MQ=255)                                          cATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCtggc                              <  1:1104124/71‑4 (MQ=255)                                           aTCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTgg                                                  >  1:2202266/1‑50 (MQ=255)                                              gggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATa                          <  1:1332168/71‑1 (MQ=255)                                                 ttAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCg                                                       <  1:1314312/39‑1 (MQ=255)                                                   aGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACaa                     >  1:4628740/1‑71 (MQ=255)                                                     cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                   >  1:5220220/1‑71 (MQ=255)                                                     cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                   <  1:3273572/71‑1 (MQ=255)                                                      aCCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                  <  1:1479695/71‑1 (MQ=255)                                                       cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                  >  1:2916497/1‑70 (MQ=255)                                                         cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg               <  1:1196942/71‑1 (MQ=255)                                                         cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg               <  1:4386550/71‑1 (MQ=255)                                                           gTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGttt                                    <  1:5369907/48‑1 (MQ=255)                                                             ccGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAcc                             >  1:321984/1‑53 (MQ=255)                                                              cGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGt           >  1:3633108/1‑70 (MQ=255)                                                              cGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg          <  1:899285/71‑1 (MQ=255)                                                                aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg               <  1:5010659/64‑1 (MQ=255)                                                                aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg               <  1:2862864/64‑1 (MQ=255)                                                                aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGAt        >  1:3351370/1‑71 (MQ=255)                                                                  cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATgg      >  1:3361145/1‑71 (MQ=255)                                                                    gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt     >  1:2058523/1‑70 (MQ=255)                                                                    gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt     >  1:1669343/1‑70 (MQ=255)                                                                    gcgcgcCAGCACGACGACCGTCTGGGGTGTTAACCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt     >  1:2981590/1‑70 (MQ=255)                                                                      gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                        <  1:3983295/49‑1 (MQ=255)                                                                      gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                        <  1:2340756/49‑1 (MQ=255)                                                                      gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                  <  1:4115035/55‑1 (MQ=255)                                                                      gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                  <  1:210484/55‑1 (MQ=255)                                                                       cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg  >  1:2490241/1‑70 (MQ=255)                                                                       |                                                                      GAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAG  >  Exported/608302‑608441

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