Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,304	A→C	29.0%	T650T (ACT→ACG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,304	0	A	C	29.0%	74.9 / 64.0	76	T650T (ACT→ACG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base C (16/6);  total (46/30)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGC  >  Exported/608234‑608374                                                                       |                                                                       cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                                        <  1:184186/71‑1 (MQ=255)      ggAACGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                                                    >  1:3755348/1‑70 (MQ=255)      ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                   <  1:1390292/71‑1 (MQ=255)      ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                   <  1:5653632/71‑1 (MQ=255)        aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                    >  1:4607131/1‑68 (MQ=255)         aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                                  >  1:513774/1‑67 (MQ=255)         aCGATTGAGAAGGTGTACGAGATCCCATATTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                                  >  1:4234740/1‑67 (MQ=255)           gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                                        <  1:4060079/61‑1 (MQ=255)           gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                                        <  1:5255274/61‑1 (MQ=255)            aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                             <  1:1327759/71‑1 (MQ=255)             ttgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                                             >  1:5678493/3‑70 (MQ=255)             tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                             <  1:4135293/70‑1 (MQ=255)             tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                             >  1:1762916/1‑70 (MQ=255)             tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                             >  1:3348759/1‑70 (MQ=255)              agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:2969125/1‑71 (MQ=255)              agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:4378990/1‑71 (MQ=255)              agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:5573347/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)                      tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                    <  1:2288181/54‑1 (MQ=255)                      tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                    <  1:5068583/54‑1 (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)                          cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                                                    <  1:813430/49‑1 (MQ=255)                           gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                              >  1:4114886/1‑68 (MQ=255)                           gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                              <  1:834402/71‑4 (MQ=255)                           gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                              >  1:1693457/1‑68 (MQ=255)                           gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:4826911/1‑58 (MQ=255)                           gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                           >  1:3596066/1‑58 (MQ=255)                              aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                                                      <  1:3959556/44‑1 (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)                                      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:3138578/1‑70 (MQ=255)                                                    aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                     <  1:945225/71‑2 (MQ=255)                                                    aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                     >  1:2844125/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:1762056/71‑1 (MQ=255)                                                         aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg                <  1:3073683/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:2368787/1‑58 (MQ=255)                                                                    ggAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAccc                  >  1:1922208/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)                                                                       |                                                                       CGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGC  >  Exported/608234‑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: ‑