Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,310	A→T	28.6%	S648S (TCT→TCA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,310	0	A	T	28.6%	93.9 / 46.5	63	S648S (TCT→TCA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (27/18);  new base T (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 = 8.92e-01

ATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCAC  >  Exported/608245‑608380                                                                  |                                                                       aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                                   <  1:1327759/71‑1 (MQ=255)  ttgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                                                   >  1:5678493/3‑70 (MQ=255)  tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                                   >  1:1762916/1‑70 (MQ=255)  tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                                   <  1:4135293/70‑1 (MQ=255)  tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                                   >  1:3348759/1‑70 (MQ=255)   agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                                 >  1:5573347/1‑71 (MQ=255)   agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                                 >  1:2969125/1‑71 (MQ=255)   agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                                 >  1:4378990/1‑71 (MQ=255)      gaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                              <  1:920718/71‑1 (MQ=255)       aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                                             <  1:3033685/71‑1 (MQ=255)       aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                                             <  1:3026919/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:834402/71‑4 (MQ=255)                gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                                    >  1:4114886/1‑68 (MQ=255)                gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                                    >  1:1693457/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:489970/1‑71 (MQ=255)                       cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                                             >  1:196632/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)                            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:321985/1‑70 (MQ=255)                                         aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                           >  1:3138578/1‑70 (MQ=255)                                         aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                                  >  1:2937546/1‑48 (MQ=255)                                         aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGc                            >  1:1711072/1‑70 (MQ=255)                                           ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                                       <  1:3407023/57‑1 (MQ=255)                                           ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                                       <  1:2386775/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)                                                                  agagGCTACTGCACCTTTCTGGTCACGACCGTGCAt                                     >  1:2995951/1‑36 (MQ=255)                                                                  agagGCTACTGCACCTTTCTGGTCACGACCGTGCAt                                     >  1:3128837/1‑36 (MQ=255)                                                                  agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac  <  1:2490337/71‑1 (MQ=255)                                                                  agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac  <  1:1399283/71‑1 (MQ=255)                                                                  |                                                                       ATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCAC  >  Exported/608245‑608380

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