Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,040	T→C	13.3%	V738V (GTA→GTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,040	0	T	C	13.3%	185.4 / 17.7	75	V738V (GTA→GTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (37/28);  new base C (4/6);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCA  >  Exported/607970‑608107                                                                        |                                                                    aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                                     <  1:3026259/71‑1 (MQ=255)  ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                                     >  1:4396549/1‑70 (MQ=255)  ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                                    >  1:2300988/1‑71 (MQ=255)  ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGATGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                                     >  1:3494523/1‑70 (MQ=255)     cATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACt                                                                  <  1:4815956/70‑1 (MQ=255)       tAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGGCAGCGAGTTGAAACGTACTGc                                                                <  1:3749529/70‑1 (MQ=255)        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCg                                                               >  1:3536577/1‑70 (MQ=255)           ttGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                                                           <  1:2617571/71‑1 (MQ=255)             gAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGAGAGTTGAAACGTACTGCGTAGcc                                                          <  1:83408/70‑1 (MQ=255)               gtaaaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGa                                                       <  1:350587/68‑1 (MQ=255)                tGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                                     <  1:1294776/56‑1 (MQ=255)                  aaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                                                     <  1:404675/70‑3 (MQ=255)                  aaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                                                     <  1:2969174/70‑3 (MQ=255)                  aaGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                     >  1:5315911/1‑70 (MQ=255)                  aaGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                     >  1:227011/1‑70 (MQ=255)                      tACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac                                                  >  1:3279435/1‑69 (MQ=255)                      tACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                 >  1:398507/1‑70 (MQ=255)                          aGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACTGATgg                                            <  1:87811/71‑1 (MQ=255)                          aGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg                                            <  1:3814349/71‑1 (MQ=255)                            taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                <  1:5158368/65‑1 (MQ=255)                            taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                <  1:3326983/65‑1 (MQ=255)                            taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                                                 <  1:366473/64‑1 (MQ=255)                             aataaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGCATTGTCa                                         <  1:569736/71‑1 (MQ=255)                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                              >  1:2715948/1‑65 (MQ=255)                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                              >  1:3593189/1‑65 (MQ=255)                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                              >  1:5538430/1‑65 (MQ=255)                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                              >  1:531529/1‑65 (MQ=255)                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                        <  1:472346/71‑1 (MQ=255)                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                        >  1:1064665/1‑71 (MQ=255)                              ataaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAg                                        >  1:4040379/1‑71 (MQ=255)                               taaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                         <  1:3450180/69‑1 (MQ=255)                                aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                                    >  1:3235411/1‑41 (MQ=255)                                aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                                    >  1:739180/1‑41 (MQ=255)                                aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                                    >  1:2849456/1‑41 (MQ=255)                                aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAg                                                            <  1:1961296/49‑1 (MQ=255)                                 aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGt                                                              >  1:1115134/1‑46 (MQ=255)                                 aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                >  1:2975706/1‑60 (MQ=255)                                 aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                         >  1:2340301/1‑67 (MQ=255)                                   gTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                   >  1:5527442/1‑71 (MQ=255)                                   gTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                   >  1:3563115/1‑71 (MQ=255)                                    tCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                           <  1:5225139/62‑1 (MQ=255)                                     ccTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTAc                                                                   >  1:2986568/1‑37 (MQ=255)                                      ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                         >  1:753132/1‑62 (MQ=255)                                      ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                 >  1:1178530/1‑70 (MQ=255)                                      ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                 >  1:3375138/1‑70 (MQ=255)                                      ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                 >  1:2966344/1‑70 (MQ=255)                                      ctgctgca‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                 >  1:562899/1‑70 (MQ=255)                                       tgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                      <  1:1232043/64‑1 (MQ=255)                                         ctgctgGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                                                >  1:4867945/1‑37 (MQ=38)                                            ct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                                          <  1:1052531/39‑1 (MQ=255)                                            ct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                          <  1:1119205/71‑1 (MQ=255)                                                ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                        <  1:1488244/70‑1 (MQ=255)                                                 tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                      >  1:2144533/1‑71 (MQ=255)                                                 tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                      >  1:5385737/1‑71 (MQ=255)                                                 tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                      >  1:3291109/1‑71 (MQ=255)                                                  cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                                       >  1:4501310/1‑37 (MQ=255)                                                  cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                                       >  1:5045808/1‑37 (MQ=255)                                                  cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                     >  1:2851405/1‑39 (MQ=255)                                                  cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                                     >  1:2000940/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:1067665/71‑1 (MQ=255)                                                        tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc               <  1:4341197/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:1169838/53‑1 (MQ=255)                                                                ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                         <  1:1855915/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)                                                                        |                                                                    ATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCA  >  Exported/607970‑608107

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