Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,376	A→G	21.3%	A626A (GCT→GCC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,376	0	A	G	21.3%	112.3 / 27.1	61	A626A (GCT→GCC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (27/21);  new base G (5/8);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

AGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAA  >  Exported/608310‑608443                                                                   |                                                                    agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                                 <  1:1399283/71‑1 (MQ=255) agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                                 <  1:2490337/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:3947601/1‑71 (MQ=255)                      gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                                            >  1:5106630/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:3273572/71‑1 (MQ=255)                                             cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                     >  1:5220220/1‑71 (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:2340756/49‑1 (MQ=255)                                                              gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                          <  1:3983295/49‑1 (MQ=255)                                                              gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                    <  1:210484/55‑1 (MQ=255)                                                              gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                    <  1:4115035/55‑1 (MQ=255)                                                               cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg    >  1:2490241/1‑70 (MQ=255)                                                                gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa  >  1:781163/1‑71 (MQ=255)                                                                   |                                                                    AGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAA  >  Exported/608310‑608443

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