Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,309	G→A	28.6%	L649L (CTG→TTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,309	0	G	A	28.6%	89.9 / 40.2	63	L649L (CTG→TTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (27/18);  new base A (13/5);  total (40/23)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGC  >  Exported/608239‑608374                                                                       |                                                                  ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                   <  1:1390292/71‑1 (MQ=255) ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                   <  1:5653632/71‑1 (MQ=255)       aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                             <  1:1327759/71‑1 (MQ=255)        ttgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                                             >  1:5678493/3‑70 (MQ=255)        tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                             >  1:3348759/1‑70 (MQ=255)        tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                             <  1:4135293/70‑1 (MQ=255)        tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                             >  1:1762916/1‑70 (MQ=255)         agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:4378990/1‑71 (MQ=255)         agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:5573347/1‑71 (MQ=255)         agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:2969125/1‑71 (MQ=255)            gaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                        <  1:920718/71‑1 (MQ=255)             aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                                       <  1:3026919/71‑1 (MQ=255)             aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                                       <  1:3033685/71‑1 (MQ=255)                gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                             <  1:15115/62‑1 (MQ=255)                gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCAcc                                                    >  1:397857/1‑71 (MQ=255)                  gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                                                  >  1:267150/1‑71 (MQ=255)                   tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                  >  1:324205/1‑54 (MQ=255)                     ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                                >  1:5311233/1‑54 (MQ=255)                     ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                               >  1:415944/1‑71 (MQ=255)                      gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                              >  1:1693457/1‑68 (MQ=255)                      gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                              <  1:834402/71‑4 (MQ=255)                      gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                              >  1:4114886/1‑68 (MQ=255)                      gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:4826911/1‑58 (MQ=255)                      gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:3596066/1‑58 (MQ=255)                           aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                         >  1:5691076/1‑71 (MQ=255)                            ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                                         >  1:3762863/1‑70 (MQ=255)                            ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                                         >  1:3321072/1‑70 (MQ=255)                            ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                             <  1:1633976/66‑1 (MQ=255)                            ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                         >  1:526394/1‑70 (MQ=255)                             cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                                       >  1:196632/1‑71 (MQ=255)                             cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                                       >  1:489970/1‑71 (MQ=255)                             cATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                  <  1:3901048/44‑1 (MQ=255)                              aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                                      <  1:1258203/71‑1 (MQ=255)                               tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGATCACGAc                                      >  1:432341/1‑70 (MQ=255)                                cTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAcca                                    >  1:5093733/1‑70 (MQ=255)                                 tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                  >  1:4230043/1‑40 (MQ=255)                                 tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                  <  1:1066885/40‑1 (MQ=255)                                 tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                   >  1:3917764/1‑39 (MQ=255)                                  ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                                   >  1:3836571/1‑70 (MQ=255)                                  ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                  >  1:1192396/1‑71 (MQ=255)                                  ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCCTg                                  >  1:4201025/1‑71 (MQ=255)                                   tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                     <  1:635675/67‑1 (MQ=255)                                   tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                                   <  1:2513253/69‑1 (MQ=255)                                   tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                  >  1:3385579/1‑70 (MQ=255)                                      cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                              <  1:3595044/71‑1 (MQ=255)                                       gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                             >  1:5471197/1‑71 (MQ=255)                                        tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                             <  1:2509113/70‑1 (MQ=255)                                          aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                                      >  1:3140718/1‑59 (MQ=255)                                          aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                                      >  1:4945308/1‑59 (MQ=255)                                          aGCAAACGGCAGTTTAGAAACGGAAGTCAGAGAGGCTACTCCACCTTTCTGGTAACGACCGTGCATCGGGt                          >  1:4693906/1‑71 (MQ=255)                                              aaCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGc                      >  1:2058665/1‑71 (MQ=255)                                               aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                     >  1:2844125/1‑70 (MQ=255)                                               aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                     <  1:945225/71‑2 (MQ=255)                                               aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                     >  1:3138578/1‑70 (MQ=255)                                               aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                     >  1:321985/1‑70 (MQ=255)                                               aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                            >  1:2937546/1‑48 (MQ=255)                                               aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGc                      >  1:1711072/1‑70 (MQ=255)                                                 ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                                 <  1:2386775/57‑1 (MQ=255)                                                 ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                                 <  1:3407023/57‑1 (MQ=255)                                                  gCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGccc                  <  1:1432672/71‑1 (MQ=255)                                                    aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg                <  1:3073683/71‑1 (MQ=255)                                                    aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg                <  1:1762056/71‑1 (MQ=255)                                                          gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa          <  1:2120626/71‑1 (MQ=255)                                                           cAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                         >  1:1451238/1‑39 (MQ=255)                                                              cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc      <  1:5113031/71‑1 (MQ=255)                                                               ggAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAccc                  >  1:1922208/1‑58 (MQ=255)                                                               ggAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAccc                  >  1:2368787/1‑58 (MQ=255)                                                                gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgc    <  1:3688099/71‑1 (MQ=255)                                                                  aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc  >  1:3105067/1‑71 (MQ=255)                                                                  aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc  >  1:980134/1‑71 (MQ=255)                                                                    tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                          <  1:5463821/45‑1 (MQ=255)                                                                       |                                                                  GGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGC  >  Exported/608239‑608374

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑