Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,055	T→G	17.9%	V733V (GTA→GTC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,055	0	T	G	17.9%	134.5 / 29.6	67	V733V (GTA→GTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base G (9/3);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATC  >  Exported/607987‑608124                                                                     |                                                                      aaGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                                      >  1:5315911/1‑70 (MQ=255) aaGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                                      >  1:227011/1‑70 (MQ=255)     tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac                                                                   >  1:3279435/1‑69 (MQ=255)     tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                  >  1:398507/1‑70 (MQ=255)         aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACTGATgg                                                             <  1:87811/71‑1 (MQ=255)         aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg                                                             <  1:3814349/71‑1 (MQ=255)           taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                                 <  1:5158368/65‑1 (MQ=255)           taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                                 <  1:3326983/65‑1 (MQ=255)           taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                                                                  <  1:366473/64‑1 (MQ=255)            aataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGCATTGTCa                                                          <  1:569736/71‑1 (MQ=255)             ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                                               >  1:2715948/1‑65 (MQ=255)             ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                                               >  1:3593189/1‑65 (MQ=255)             ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                                               >  1:531529/1‑65 (MQ=255)             ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                                               >  1:5538430/1‑65 (MQ=255)             ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                         >  1:1064665/1‑71 (MQ=255)             ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                         <  1:472346/71‑1 (MQ=255)             ataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAg                                                         >  1:4040379/1‑71 (MQ=255)              taaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                                          <  1:3450180/69‑1 (MQ=255)                aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                                 >  1:2975706/1‑60 (MQ=255)                aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                                          >  1:2340301/1‑67 (MQ=255)                  gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                    >  1:5527442/1‑71 (MQ=255)                  gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                    >  1:3563115/1‑71 (MQ=255)                   tCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                            <  1:5225139/62‑1 (MQ=255)                     ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                                          >  1:753132/1‑62 (MQ=255)                     ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                  >  1:2966344/1‑70 (MQ=255)                     ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                  >  1:1178530/1‑70 (MQ=255)                     ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                  >  1:3375138/1‑70 (MQ=255)                     ctgctgcagTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                  >  1:562899/1‑70 (MQ=255)                      tgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                       <  1:1232043/64‑1 (MQ=255)                           ctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                                           <  1:1119205/71‑1 (MQ=255)                              ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                         <  1:1488244/70‑1 (MQ=255)                               tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                       >  1:3291109/1‑71 (MQ=255)                               tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                       >  1:2144533/1‑71 (MQ=255)                               tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                       >  1:5385737/1‑71 (MQ=255)                                cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                                      >  1:2000940/1‑39 (MQ=255)                                cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                                      >  1:2851405/1‑39 (MQ=255)                                cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                         >  1:4205704/1‑52 (MQ=255)                                    aGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                                           <  1:5332446/62‑1 (MQ=255)                                     gTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCgg                                       >  1:3904211/1‑65 (MQ=9)                                     gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc                                 <  1:2381112/71‑1 (MQ=255)                                     gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc                                 >  1:4713787/1‑71 (MQ=255)                                      tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc                                <  1:4341197/71‑1 (MQ=255)                                      tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc                                <  1:1067665/71‑1 (MQ=255)                                        aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                                    <  1:862411/65‑1 (MQ=255)                                        aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                                    <  1:5478999/65‑1 (MQ=255)                                        aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                              >  1:4465790/1‑71 (MQ=255)                                         gCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                            <  1:2426693/40‑1 (MQ=255)                                          cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATcg                            >  1:2599153/1‑71 (MQ=255)                                            cgTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTcc                                         <  1:5500462/55‑1 (MQ=12)                                              ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                                          <  1:1855915/53‑1 (MQ=255)                                              ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                                          <  1:1169838/53‑1 (MQ=255)                                                gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCa                               >  1:3090896/1‑62 (MQ=255)                                                gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                       <  1:1841906/70‑1 (MQ=255)                                                  aaCGTACTGCGTAGCCAGATACAGGGATGGTCAGCTGCGGATATTTTTcc                                         <  1:4784597/50‑1 (MQ=255)                                                   aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                   <  1:2574704/71‑1 (MQ=255)                                                   aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                     >  1:4723853/1‑69 (MQ=255)                                                      tACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttt                                             <  1:1550897/42‑1 (MQ=255)                                                        cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt              >  1:2496566/3‑71 (MQ=12)                                                        cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAt               >  1:2145640/3‑70 (MQ=12)                                                          gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                         <  1:4827450/42‑1 (MQ=255)                                                          gCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTc            >  1:3894100/1‑71 (MQ=16)                                                           cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                              >  1:4987806/1‑52 (MQ=255)                                                           cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt             <  1:2534990/69‑1 (MQ=255)                                                           cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa           >  1:686846/1‑71 (MQ=255)                                                              aGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga        >  1:611247/1‑70 (MQ=12)                                                               gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCgat       >  1:204410/1‑69 (MQ=12)                                                                  agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtc     <  1:3476824/70‑1 (MQ=9)                                                                  aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca    >  1:4707475/1‑71 (MQ=255)                                                                     tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc  <  1:4459557/70‑1 (MQ=255)                                                                     tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc  <  1:3986411/70‑1 (MQ=255)                                                                     |                                                                      AAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATC  >  Exported/607987‑608124

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