Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,341	T→G	16.7%	M305L (ATG→CTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,341	0	T	G	16.7%	122.1 / 26.5	66	M305L (ATG→CTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (33/22);  new base G (3/8);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 9.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

ACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAA  >  Exported/609272‑609407                                                                      |                                                                   aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt                                                                   >  1:1719359/1‑71 (MQ=255) aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt                                                                   >  1:4089161/1‑71 (MQ=255)  cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                   >  1:2999037/1‑70 (MQ=255)       cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCGACCATTTCCtac                                                             >  1:3684893/1‑69 (MQ=255)       cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                                             >  1:3693472/1‑71 (MQ=255)       cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                                             >  1:1981534/1‑71 (MQ=255)          caaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgc                                                           >  1:1560493/4‑70 (MQ=255)           gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGctt                                                         >  1:4107251/1‑71 (MQ=255)           aaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcc                                                          <  1:939224/68‑2 (MQ=255)              ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGcttctt                                                      >  1:5283728/1‑71 (MQ=255)               tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                   <  1:4064226/57‑1 (MQ=255)               tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                   <  1:5318468/57‑1 (MQ=255)                 cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                   <  1:236672/71‑1 (MQ=255)                 cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                   <  1:4578514/71‑1 (MQ=255)                   gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                              >  1:1729511/3‑58 (MQ=255)                   gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                              >  1:4909761/3‑58 (MQ=255)                   gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAttt                                                                  <  1:4530147/54‑1 (MQ=255)                      tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                         <  1:4567928/60‑3 (MQ=255)                      tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                         <  1:3122764/60‑3 (MQ=255)                      tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                         <  1:339874/60‑3 (MQ=255)                      tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                         <  1:1205687/60‑3 (MQ=255)                      tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                         <  1:1301603/60‑3 (MQ=255)                      tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                         <  1:2205378/60‑3 (MQ=255)                      tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                         <  1:4251922/60‑3 (MQ=255)                        cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                 >  1:299515/1‑50 (MQ=255)                        cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                 >  1:1787213/1‑50 (MQ=255)                           aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                 >  1:1682521/1‑47 (MQ=255)                           aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                 >  1:4362599/1‑47 (MQ=255)                           aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                 >  1:2102882/1‑47 (MQ=255)                            ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                   <  1:811417/44‑1 (MQ=255)                            ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                   <  1:242077/44‑1 (MQ=255)                            ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                        >  1:1523225/1‑71 (MQ=255)                            ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                        >  1:299446/1‑71 (MQ=255)                            ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                        <  1:3169799/71‑1 (MQ=255)                             gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCt                                          >  1:4128862/1‑68 (MQ=255)                             gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                                       >  1:5601422/1‑71 (MQ=255)                               aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                        >  1:5283690/1‑68 (MQ=255)                               aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                                       >  1:2070553/1‑69 (MQ=255)                                  gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATc                                           <  1:4481714/62‑1 (MQ=255)                                  gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                                  >  1:990197/1‑71 (MQ=255)                                  gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                                  <  1:1970934/71‑1 (MQ=255)                                     ccATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                 >  1:3614193/1‑37 (MQ=255)                                       aTACGCAGTTTCTTGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                   <  1:30087/49‑1 (MQ=255)                                          cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc                          <  1:297722/71‑1 (MQ=255)                                          cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc                          <  1:4365751/71‑1 (MQ=255)                                              gTTTCATGACCAGGTGGTCTACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGt                      <  1:5111201/71‑1 (MQ=255)                                               tttCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCttt                                <  1:960013/60‑1 (MQ=255)                                                 tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                                         >  1:4065374/1‑49 (MQ=255)                                                 tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                                         >  1:422164/1‑49 (MQ=255)                                                 tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                    <  1:2058404/70‑1 (MQ=255)                                                 tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                    <  1:1026504/70‑1 (MQ=255)                                                 tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                    <  1:1665705/70‑1 (MQ=255)                                                 tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                    <  1:4254341/70‑1 (MQ=255)                                                 tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                    >  1:3027650/1‑70 (MQ=255)                                                 tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                    >  1:728656/1‑70 (MQ=255)                                                 tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                    <  1:333195/70‑1 (MQ=255)                                                    tGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATg                >  1:178792/1‑71 (MQ=255)                                                        cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                              >  1:2169634/1‑53 (MQ=255)                                                        cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg                             >  1:1121064/1‑54 (MQ=255)                                                        cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg                             >  1:2332258/1‑54 (MQ=255)                                                          ggtggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                              >  1:680557/1‑51 (MQ=255)                                                           gtggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACAt         >  1:5601343/1‑71 (MQ=255)                                                             ggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCa                                     >  1:2663644/1‑41 (MQ=255)                                                              gtCAACCATTTCCTGCTCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATcc       <  1:1769834/70‑1 (MQ=255)                                                              gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATc        >  1:4725418/1‑69 (MQ=255)                                                              gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATc        >  1:2417452/1‑69 (MQ=255)                                                              gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCa      <  1:4826755/71‑1 (MQ=255)                                                                cAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTa              <  1:5674027/61‑1 (MQ=255)                                                                 aatcATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAgg    >  1:3897824/4‑70 (MQ=255)                                                                 aaCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAgg    >  1:4826953/1‑70 (MQ=255)                                                                 aaCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAgg    <  1:3425056/70‑1 (MQ=255)                                                                 aaCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGa   <  1:1648544/71‑1 (MQ=255)                                                                 aaCCATTTACTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGa   <  1:2019610/71‑1 (MQ=255)                                                                  aCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGGCAGCTTTCAGGTCACGTTc                    <  1:275840/53‑1 (MQ=255)                                                                      tttCTTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                              <  1:5024448/39‑1 (MQ=255)                                                                      tttCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCttt                                <  1:3752275/37‑1 (MQ=255)                                                                      tttCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCttt                                <  1:4387532/37‑1 (MQ=255)                                                                      tttCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTc                               <  1:2860847/38‑1 (MQ=255)                                                                      tttCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATg                <  1:5584527/53‑1 (MQ=255)                                                                      tttCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGaa  <  1:2833334/67‑1 (MQ=255)                                                                      tttCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGa   <  1:5449981/66‑1 (MQ=255)                                                                      |                                                                   ACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAA  >  Exported/609272‑609407

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