Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,078	G→C	10.1%	S392S (TCC→TCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,078	0	G	C	10.1%	218.1 / 12.9	81	S392S (TCC→TCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base C (5/3);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

TAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTT  >  Exported/609010‑609148                                                                     |                                                                       tAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGa                                                                      <  1:3301044/71‑1 (MQ=255)    tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                                                   <  1:1115463/68‑1 (MQ=255)    tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                                      >  1:3663707/4‑68 (MQ=255)    tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                                      >  1:2152665/4‑68 (MQ=255)    tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                                      >  1:1970741/4‑68 (MQ=255)      aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                 >  1:445539/1‑71 (MQ=255)      aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                 >  1:3379225/1‑71 (MQ=255)      aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                 >  1:3179307/1‑71 (MQ=255)      aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                 >  1:705335/1‑71 (MQ=255)       tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                 >  1:3116851/1‑70 (MQ=255)       tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                 >  1:3110437/1‑70 (MQ=255)        cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                 >  1:2781174/1‑69 (MQ=255)         gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                                                               <  1:845802/70‑1 (MQ=255)         gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAg                                                              >  1:5548385/1‑71 (MQ=255)           tgttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                                            >  1:456444/1‑71 (MQ=255)                    ccGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAt                                                      <  1:1150458/68‑1 (MQ=255)                     cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                                   <  1:5337889/70‑1 (MQ=255)                     cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                                   <  1:2198931/70‑1 (MQ=255)                       gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                <  1:3077275/71‑1 (MQ=255)                       gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                <  1:5231204/71‑1 (MQ=255)                       gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGAAACTTTAGCGGCGAATTTCtt                                                 >  1:249047/1‑70 (MQ=255)                        aCGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                                               <  1:3909367/71‑1 (MQ=255)                         cGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                >  1:1036772/1‑69 (MQ=255)                         cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                                                   >  1:1559307/1‑50 (MQ=12)                         cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTg                                                >  1:5175063/1‑69 (MQ=255)                         cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGa                                               <  1:1731496/70‑1 (MQ=255)                           cATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                                            >  1:1500114/1‑55 (MQ=255)                           cATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGAAgg                                            <  1:1224327/71‑2 (MQ=255)                            aTCAGATCGTCGTTTTCATACTGCAAAGAGGATGTGACGATAGACACCTGCGCGGCATATTTTTTGaa                                              >  1:3717520/1‑68 (MQ=255)                             tCAGTTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                          <  1:3847743/71‑1 (MQ=255)                             tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                <  1:2152486/65‑1 (MQ=255)                             tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                          <  1:2073074/71‑1 (MQ=255)                             tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                          <  1:4503201/71‑1 (MQ=255)                             tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                          <  1:4737692/71‑1 (MQ=255)                              cAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGtt                                           >  1:4818782/1‑69 (MQ=255)                               aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                                              >  1:1385118/1‑65 (MQ=255)                               aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                                              >  1:2111126/1‑65 (MQ=255)                               aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                                               >  1:2873580/1‑64 (MQ=255)                                 gTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                                                               <  1:5636241/46‑1 (MQ=255)                                  tCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                                     <  1:4179864/71‑1 (MQ=255)                                   cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                 >  1:3452629/1‑42 (MQ=255)                                   cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                                     >  1:5694392/1‑70 (MQ=255)                                    aTCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGc                                    <  1:5530304/70‑1 (MQ=255)                                     tCGTTCTCATTCTGCAGAGAAGAGGTGTCGATGGACACTTTAGGGGCGAATttgttg                                                <  1:3695301/57‑1 (MQ=255)                                     tCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAg                                  <  1:1416829/71‑1 (MQ=255)                                     tCGTTCTCATACTGCAGAGAAGAGCTGTCGATGGa                                                                      <  1:4657500/35‑1 (MQ=255)                                       gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGGGGCAGtt                                >  1:145938/1‑71 (MQ=255)                                       gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtt                                >  1:3218500/1‑71 (MQ=255)                                        ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt                               <  1:5655573/71‑1 (MQ=255)                                        ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt                               <  1:701365/71‑1 (MQ=255)                                          ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                       >  1:3134255/1‑45 (MQ=255)                                          ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                       >  1:2721395/1‑45 (MQ=255)                                          ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                              <  1:5499503/70‑1 (MQ=255)                                          ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                              <  1:986105/70‑1 (MQ=255)                                          ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                              >  1:1640732/1‑70 (MQ=255)                                           tcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                <  1:2005458/51‑1 (MQ=255)                                           tcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                          <  1:2449680/57‑1 (MQ=255)                                            cATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                             >  1:3973005/1‑69 (MQ=255)                                            cATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAATTCTTGAAGTTCAGCGGCAGttttt                             >  1:3973945/1‑69 (MQ=255)                                             aTACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                             <  1:4208587/68‑1 (MQ=255)                                             aTACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                             <  1:2829398/68‑1 (MQ=255)                                                   cAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGt                                            >  1:3659735/1‑47 (MQ=255)                                                   cAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACca                      >  1:4997603/1‑69 (MQ=255)                                                    agagAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAg                   <  1:2205120/71‑1 (MQ=255)                                                     gagaAGAGGTGTCGATGGGCACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGa                  <  1:2248125/71‑1 (MQ=255)                                                     gagaAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCa                                   >  1:439162/1‑54 (MQ=255)                                                     gagaAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCa                                   >  1:3008616/1‑54 (MQ=255)                                                      agaagaGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAg                   <  1:1388042/69‑1 (MQ=255)                                                         agagGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATgg              <  1:143769/71‑1 (MQ=255)                                                             gtgtCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGt             <  1:944340/68‑1 (MQ=255)                                                              tgtCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                                     <  1:5172726/43‑1 (MQ=255)                                                               gtCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGt        >  1:5575644/1‑71 (MQ=255)                                                                tCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGtt       >  1:3521916/1‑71 (MQ=255)                                                                tCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGt        <  1:2777022/70‑1 (MQ=255)                                                                 cGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGaa                 >  1:5057120/1‑60 (MQ=255)                                                                  gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAg                                        >  1:4742184/1‑36 (MQ=255)                                                                  gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATg         <  1:4903770/67‑1 (MQ=255)                                                                  gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATg         <  1:3906427/67‑1 (MQ=255)                                                                  gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATg         <  1:4558037/67‑1 (MQ=255)                                                                  gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCg     >  1:1249222/1‑71 (MQ=255)                                                                    tGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACcac                     <  1:4296934/53‑1 (MQ=255)                                                                     ggACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGtt  >  1:908052/1‑71 (MQ=255)                                                                     |                                                                       TAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTT  >  Exported/609010‑609148

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