Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,327	C→G	15.9%	L309L (CTG→CTC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,327	0	C	G	15.9%	158.3 / 22.0	69	L309L (CTG→CTC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base G (3/8);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 9.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

TCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATG  >  Exported/609257‑609393                                                                       |                                                                   tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                                                    <  1:3888221/71‑1 (MQ=255)  cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                                                    <  1:2468755/70‑1 (MQ=255)  cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                                                    <  1:5527717/70‑1 (MQ=255)       cGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtg                                                               >  1:4023666/1‑70 (MQ=255)        gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                                                             >  1:3271539/1‑71 (MQ=255)        gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                                                             >  1:5303577/1‑71 (MQ=255)          cAGAGAACAGTTCATCGTATTCCTGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                           <  1:5031927/71‑1 (MQ=255)          cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                           <  1:278814/71‑1 (MQ=255)           agagAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAg                                                                  >  1:1972603/1‑63 (MQ=255)            gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                                          >  1:878281/1‑69 (MQ=255)            gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                                          >  1:2931471/1‑69 (MQ=255)              gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa                                                       <  1:1885619/71‑1 (MQ=255)                aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt                                                     >  1:1719359/1‑71 (MQ=255)                aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt                                                     >  1:4089161/1‑71 (MQ=255)                 cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                     >  1:2999037/1‑70 (MQ=255)                      cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCGACCATTTCCtac                                               >  1:3684893/1‑69 (MQ=255)                      cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                               >  1:3693472/1‑71 (MQ=255)                      cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                               >  1:1981534/1‑71 (MQ=255)                         caaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgc                                             >  1:1560493/4‑70 (MQ=255)                          gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGctt                                           >  1:4107251/1‑71 (MQ=255)                          aaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcc                                            <  1:939224/68‑2 (MQ=255)                             ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTcta                                                          >  1:958715/1‑51 (MQ=255)                             ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGcttctt                                        >  1:5283728/1‑71 (MQ=255)                              tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                     <  1:4064226/57‑1 (MQ=255)                              tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                     <  1:5318468/57‑1 (MQ=255)                               ccGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCaa                                                          >  1:1607554/1‑51 (MQ=255)                                cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                     <  1:236672/71‑1 (MQ=255)                                cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                     <  1:4578514/71‑1 (MQ=255)                                  gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                >  1:1729511/3‑58 (MQ=255)                                  gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                >  1:4909761/3‑58 (MQ=255)                                  gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                           <  1:3155037/47‑1 (MQ=255)                                  gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAttt                                                    <  1:4530147/54‑1 (MQ=255)                                     tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                           <  1:3122764/60‑3 (MQ=255)                                     tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                           <  1:4251922/60‑3 (MQ=255)                                     tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                           <  1:4567928/60‑3 (MQ=255)                                     tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                           <  1:339874/60‑3 (MQ=255)                                     tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                           <  1:1205687/60‑3 (MQ=255)                                     tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                           <  1:1301603/60‑3 (MQ=255)                                     tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                           <  1:2205378/60‑3 (MQ=255)                                       cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                   >  1:299515/1‑50 (MQ=255)                                       cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                   >  1:1787213/1‑50 (MQ=255)                                          aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                   >  1:4362599/1‑47 (MQ=255)                                          aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                   >  1:1682521/1‑47 (MQ=255)                                          aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                   >  1:2102882/1‑47 (MQ=255)                                           ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                     <  1:811417/44‑1 (MQ=255)                                           ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                     <  1:242077/44‑1 (MQ=255)                                           ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                          >  1:1523225/1‑71 (MQ=255)                                           ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                          <  1:3169799/71‑1 (MQ=255)                                           ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                          >  1:299446/1‑71 (MQ=255)                                            gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCt                            >  1:4128862/1‑68 (MQ=255)                                            gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                         >  1:5601422/1‑71 (MQ=255)                                              aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                          >  1:5283690/1‑68 (MQ=255)                                              aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                         >  1:2070553/1‑69 (MQ=255)                                                 gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATc                             <  1:4481714/62‑1 (MQ=255)                                                 gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                    >  1:990197/1‑71 (MQ=255)                                                 gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                    <  1:1970934/71‑1 (MQ=255)                                                    ccATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                   >  1:3614193/1‑37 (MQ=255)                                                      aTACGCAGTTTCTTGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                     <  1:30087/49‑1 (MQ=255)                                                         cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc            <  1:4365751/71‑1 (MQ=255)                                                         cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc            <  1:297722/71‑1 (MQ=255)                                                             gTTTCATGACCAGGTGGTCTACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGt        <  1:5111201/71‑1 (MQ=255)                                                              tttCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCttt                  <  1:960013/60‑1 (MQ=255)                                                                tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                           >  1:422164/1‑49 (MQ=255)                                                                tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                           >  1:4065374/1‑49 (MQ=255)                                                                tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc      <  1:4254341/70‑1 (MQ=255)                                                                tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc      <  1:1026504/70‑1 (MQ=255)                                                                tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc      <  1:1665705/70‑1 (MQ=255)                                                                tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc      >  1:728656/1‑70 (MQ=255)                                                                tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc      <  1:2058404/70‑1 (MQ=255)                                                                tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc      >  1:3027650/1‑70 (MQ=255)                                                                tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc      <  1:333195/70‑1 (MQ=255)                                                                   tGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATg  >  1:178792/1‑71 (MQ=255)                                                                       cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                >  1:2169634/1‑53 (MQ=255)                                                                       cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg               >  1:2332258/1‑54 (MQ=255)                                                                       cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg               >  1:1121064/1‑54 (MQ=255)                                                                       |                                                                   TCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATG  >  Exported/609257‑609393

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