Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,073	T→A	13.6%	D394V (GAC→GTC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,073	0	T	A	13.6%	202.0 / 22.9	81	D394V (GAC→GTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base A (7/4);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

ATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTT  >  Exported/609004‑609143                                                                      |                                                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:1082422/71‑1 (MQ=255) aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:5588921/71‑1 (MQ=255) aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       >  1:2302400/1‑71 (MQ=255) aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:1551045/71‑1 (MQ=255)  tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCg                                                                      <  1:5383566/71‑1 (MQ=255)    gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt                                                                    <  1:2751295/71‑1 (MQ=255)    gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa                                                                     <  1:940828/70‑1 (MQ=255)    gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGAt                                                                    >  1:2137280/1‑71 (MQ=255)       tAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGa                                                                 <  1:3301044/71‑1 (MQ=255)          tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                                              <  1:1115463/68‑1 (MQ=255)          tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                                 >  1:1970741/4‑68 (MQ=255)          tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                                 >  1:2152665/4‑68 (MQ=255)          tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                                 >  1:3663707/4‑68 (MQ=255)            aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:3179307/1‑71 (MQ=255)            aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:445539/1‑71 (MQ=255)            aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:705335/1‑71 (MQ=255)            aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:3379225/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:2198931/70‑1 (MQ=255)                           cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                              <  1:5337889/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)                                   tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:611434/37‑1 (MQ=255)                                   tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:16782/37‑1 (MQ=255)                                   tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                           <  1:2152486/65‑1 (MQ=255)                                   tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                     <  1:4737692/71‑1 (MQ=255)                                   tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                     <  1:4503201/71‑1 (MQ=255)                                   tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                     <  1:2073074/71‑1 (MQ=255)                                    cAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGtt                                      >  1:4818782/1‑69 (MQ=255)                                     aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:3710638/35‑1 (MQ=255)                                     aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGa                                                                     >  1:5544640/1‑37 (MQ=255)                                     aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                                         >  1:2111126/1‑65 (MQ=255)                                     aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                                         >  1:1385118/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:701365/71‑1 (MQ=255)                                              ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt                          <  1:5655573/71‑1 (MQ=255)                                                ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                  >  1:2721395/1‑45 (MQ=255)                                                ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                  >  1:3134255/1‑45 (MQ=255)                                                ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                         <  1:5499503/70‑1 (MQ=255)                                                ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                         >  1:1640732/1‑70 (MQ=255)                                                ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                         <  1:986105/70‑1 (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:2829398/68‑1 (MQ=255)                                                   aTACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                        <  1:4208587/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)                                                                      |                                                                       ATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTT  >  Exported/609004‑609143

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