Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,194	G→A	9.2%	A687V (GCT→GTT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,194	0	G	A	9.2%	264.0 / 11.6	87	A687V (GCT→GTT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (29/50);  new base A (5/3);  total (34/53)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGA  >  Exported/608128‑608263                                                                   |                                                                      ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcc                                                                     <  1:2545555/69‑1 (MQ=255)   aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                                                                 >  1:1161895/1‑71 (MQ=255)    aCGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGtt                                                                <  1:2650862/71‑1 (MQ=255)    aCGTTCAGGTGCTGACCACCTTCGATGTATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                                >  1:221652/1‑71 (MQ=255)    aCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                                >  1:2128310/1‑71 (MQ=255)    aCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                                <  1:3235564/71‑1 (MQ=255)     cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTg                                                               <  1:2547285/71‑1 (MQ=255)     cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTg                                                               >  1:3416905/1‑71 (MQ=255)      gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTgg                                                              <  1:2992485/71‑1 (MQ=255)      gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTgg                                                              <  1:1385402/71‑1 (MQ=255)          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                           <  1:2370936/70‑1 (MQ=255)          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                           <  1:1380264/70‑1 (MQ=255)          aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAACCTGGCTGGTCtg                                                          <  1:828901/71‑2 (MQ=39)             tGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAACCTGGCTg                                                               <  1:1656682/63‑1 (MQ=255)              gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCTATCAGACCAGCCAGGTTGGTCTTACg                                                       >  1:3258761/1‑70 (MQ=255)               cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGa                                                      <  1:434264/70‑1 (MQ=255)               cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAACCTGGCTGGt                                                             >  1:895472/1‑63 (MQ=255)               cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAACCTGGCTGGt                                                             >  1:3086116/1‑63 (MQ=255)               cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAACCTGGCTGGt                                                             >  1:1586726/1‑63 (MQ=255)                tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAAc                                                    >  1:1377271/1‑71 (MQ=255)                  accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                                                                 <  1:2799403/56‑1 (MQ=255)                  accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACtt                                                  <  1:3205139/71‑1 (MQ=255)                  accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                   <  1:1969002/70‑1 (MQ=255)                   ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                 <  1:599943/71‑1 (MQ=255)                   ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                 >  1:831423/1‑71 (MQ=255)                   ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACAATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                 <  1:1910382/71‑1 (MQ=255)                   ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACAATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                 <  1:1914503/71‑1 (MQ=255)                    caccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcg                                                <  1:2452448/71‑1 (MQ=255)                       cTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                             >  1:1584842/1‑71 (MQ=255)                        ttCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                                                                 <  1:2916824/50‑1 (MQ=255)                        ttCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                                                                 <  1:1420810/50‑1 (MQ=255)                        ttCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCtt                                                          <  1:769596/57‑1 (MQ=255)                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                            <  1:158074/69‑1 (MQ=255)                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                            <  1:3281764/69‑1 (MQ=255)                          cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                                          >  1:1257412/1‑71 (MQ=255)                           gatggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                         <  1:2778301/71‑1 (MQ=255)                            atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                 <  1:547807/62‑1 (MQ=255)                             tggatgCTTCGTGGTGTAAGTAACCATCCATCAGACCAGCCAGCTTGGTCTTACGa                                                      >  1:1817731/1‑56 (MQ=39)                              ggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                      >  1:1056648/1‑71 (MQ=255)                              ggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                      >  1:978465/1‑71 (MQ=255)                              ggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                      >  1:1084660/1‑71 (MQ=255)                               gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTAGTCTTTAcc                                     <  1:1333530/71‑1 (MQ=255)                                atgCTTCGTGGTTGAAGTAACCATCCATCAGACCAGCCAGGTTGGt                                                             <  1:955561/46‑1 (MQ=255)                                atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGt                                                             <  1:1958557/46‑1 (MQ=255)                                atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGt                                                             <  1:2996878/46‑1 (MQ=255)                                atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                           <  1:2407937/48‑1 (MQ=255)                                atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                           <  1:2632856/48‑1 (MQ=255)                                atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                                    <  1:3641242/71‑1 (MQ=255)                                atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                                    <  1:1049229/71‑1 (MQ=255)                                  gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                         >  1:1416806/1‑48 (MQ=255)                                  gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                         >  1:1191422/1‑48 (MQ=255)                                  gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                         >  1:621257/1‑48 (MQ=255)                                  gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                         >  1:151684/1‑48 (MQ=255)                                  gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                         >  1:3438648/1‑48 (MQ=255)                                   cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGt                                 <  1:2006578/71‑1 (MQ=255)                                    ttCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                                >  1:3317690/1‑69 (MQ=255)                                    ttCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                                >  1:555346/1‑69 (MQ=255)                                        tggtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCa                                   >  1:2123364/1‑64 (MQ=255)                                        tggtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCa                                   >  1:1810574/1‑64 (MQ=255)                                        ctgtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                             <  1:1079716/68‑1 (MQ=255)                                           tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                           <  1:242005/37‑1 (MQ=255)                                           tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                             <  1:1129711/67‑1 (MQ=255)                                           tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                             <  1:2592002/67‑1 (MQ=255)                                           tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                             <  1:1732766/67‑1 (MQ=255)                                           tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                         <  1:1470584/71‑1 (MQ=255)                                            ggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCtt                                                          <  1:651479/37‑1 (MQ=255)                                             gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                           >  1:2901128/1‑35 (MQ=255)                                             gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                           >  1:271989/1‑35 (MQ=255)                                             gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAAc                                                    >  1:3705305/1‑42 (MQ=255)                                             gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                 <  1:2564473/45‑1 (MQ=255)                                             gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                         <  1:81576/69‑1 (MQ=255)                                               aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACGCAGTGCGTTCGGAAc                      >  1:3145403/1‑70 (MQ=255)                                               aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                         <  1:3243649/67‑1 (MQ=255)                                               aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                         <  1:2212706/67‑1 (MQ=255)                                                   aCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                           <  1:2868811/45‑1 (MQ=255)                                                     catccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                      <  1:2684017/64‑1 (MQ=255)                                                      atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag               >  1:2197010/1‑70 (MQ=255)                                                      atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag               >  1:2396999/1‑70 (MQ=255)                                                      atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga              <  1:983901/71‑1 (MQ=255)                                                      atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga              <  1:1198443/71‑1 (MQ=255)                                                      atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga              <  1:2909350/71‑1 (MQ=255)                                                        ccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAg            >  1:2451130/1‑71 (MQ=255)                                                         catcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAgg           >  1:3194138/1‑71 (MQ=255)                                                            cAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg         >  1:2079794/1‑70 (MQ=255)                                                             aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt        <  1:1355717/70‑1 (MQ=255)                                                             aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt        <  1:541844/70‑1 (MQ=255)                                                             aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt        <  1:2219934/70‑1 (MQ=255)                                                              gACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAg      <  1:189502/71‑1 (MQ=255)                                                                ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                          >  1:2560959/1‑49 (MQ=255)                                                                ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                          >  1:179306/1‑49 (MQ=255)                                                                ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                          >  1:2847790/1‑49 (MQ=255)                                                                ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAg      >  1:470641/1‑69 (MQ=255)                                                                 cagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGGACGATAGAGAAGGTGTAg      >  1:341976/1‑68 (MQ=255)                                                                  agccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAg      >  1:3269415/1‑67 (MQ=255)                                                                  agccagGTTGGTCTTACGAACTTCGTCGTCTATACCCAGTGCGTTCGGAACGATAGAGAAGGTGTa       >  1:2142812/1‑66 (MQ=255)                                                                   gccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGgaga  >  1:2956301/1‑70 (MQ=255)                                                                   |                                                                      TTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGA  >  Exported/608128‑608263

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