Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,170	A→G	12.7%	Y362H (TAC→CAC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,170	0	A	G	12.7%	156.0 / 12.5	63	Y362H (TAC→CAC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (6/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

GTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGAT  >  Exported/609105‑609239                                                                  |                                                                      gTTCAGCGGCAGTTTTTCAGACCACAGAATGGGCATGTTCGGTTCCGGAGACGGACCCATGGGGTACAggg                                                                  <  1:630309/71‑1 (MQ=255)  ttCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGgt                                                                 <  1:1283183/71‑1 (MQ=255)    cAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGgtg                                                                <  1:3811779/70‑1 (MQ=255)          cAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAgg                                                          >  1:2101799/1‑70 (MQ=255)             tttttCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACa                                                                     <  1:5682322/56‑1 (MQ=255)              ttttCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTAc                                                                      >  1:5387739/1‑54 (MQ=255)               tttCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACa                                                                     >  1:3332314/1‑54 (MQ=255)                ttCCGACCAAAGAATGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAg                                                                    <  1:1147016/54‑1 (MQ=255)                ttCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACgg                                                    <  1:5219778/70‑1 (MQ=255)                 tCAGACCACAGAATGGTCATGTTCGGCTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGa                                                   <  1:3006246/70‑1 (MQ=255)                      ccACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGt                                             >  1:5574060/1‑71 (MQ=255)                        acaGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGttt                                           <  1:5498938/71‑1 (MQ=255)                         caGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGtttt                                          <  1:3007695/71‑1 (MQ=255)                           gAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAg                                                                    <  1:5652737/43‑1 (MQ=255)                           gAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGttttt                                         >  1:324917/1‑70 (MQ=255)                            aaTGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGtttt                                          <  1:4267762/68‑1 (MQ=255)                            aaTGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGtttt                                          <  1:333433/68‑1 (MQ=255)                              tGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACa                                                                     >  1:1275916/1‑39 (MQ=255)                                 tCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGgtg                                                                >  1:1903412/1‑41 (MQ=255)                                 tCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACggg                                                   >  1:4442518/1‑53 (MQ=255)                                 tCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACgg                                                    >  1:3008823/1‑53 (MQ=255)                                   aTGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTa                                     >  1:3061706/1‑66 (MQ=255)                                   aTGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAg                                >  1:1055736/1‑71 (MQ=255)                                    tGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTc                                                             <  1:4275667/41‑1 (MQ=255)                                    tGTTCGGTTCCGGAGACGGACCCATGATGTACAGGGTGTTc                                                             <  1:4775169/41‑1 (MQ=255)                                     gTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAggg                              >  1:1284140/1‑71 (MQ=255)                                      ttCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAg                                                 <  1:3889965/51‑1 (MQ=255)                                      ttCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGtt                                            >  1:3038035/1‑56 (MQ=255)                                      ttCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTaa                                    >  1:260483/1‑64 (MQ=255)                                      ttCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAggg                              >  1:675934/1‑70 (MQ=255)                                       tCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTa                            >  1:359129/1‑71 (MQ=255)                                       tCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTa                            >  1:251769/1‑71 (MQ=255)                                        cGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTAc                           >  1:5628799/1‑71 (MQ=255)                                         ggTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGa                                                   <  1:1371647/46‑1 (MQ=255)                                         ggTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGa                                                   <  1:663643/46‑1 (MQ=255)                                         ggTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGa                                                   <  1:3070475/46‑1 (MQ=255)                                             ccGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCg                      >  1:3917869/1‑71 (MQ=255)                                              cGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCg                      <  1:4329675/70‑1 (MQ=255)                                               ggAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTc                    <  1:4430501/71‑1 (MQ=255)                                               ggAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTc                    <  1:5267283/71‑1 (MQ=255)                                                  gaCGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAg                 <  1:1688737/71‑1 (MQ=255)                                                     ggccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc              >  1:88760/4‑71 (MQ=255)                                                     ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAgg                               >  1:3280802/1‑54 (MQ=255)                                                     ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc              >  1:4277363/1‑71 (MQ=255)                                                     ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc              <  1:831753/71‑1 (MQ=255)                                                     ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc              <  1:4683451/71‑1 (MQ=255)                                                      gccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc              >  1:5585024/3‑70 (MQ=255)                                                      gccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc              >  1:1969608/3‑70 (MQ=255)                                                        cccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCAt           >  1:2876034/1‑71 (MQ=255)                                                         ccATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGa                  >  1:1727957/1‑63 (MQ=255)                                                         ccATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGa                  >  1:3538263/1‑63 (MQ=255)                                                          cATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac         <  1:1498445/71‑1 (MQ=255)                                                          cATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCAtcc         <  1:28585/71‑3 (MQ=255)                                                           aTGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac         >  1:3725049/1‑70 (MQ=255)                                                            tGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATacca       >  1:1355781/1‑71 (MQ=255)                                                              gtgtGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATcccgcc     >  1:3131886/1‑71 (MQ=255)                                                              gtgtGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATcccgc      >  1:4181715/1‑68 (MQ=255)                                                              gtgtACAGGGTGTTCAGGAAACGGAAGCTGTTTTTg                                        >  1:2325332/1‑36 (MQ=255)                                                              gtgtACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATaccac      >  1:1664215/1‑70 (MQ=255)                                                                 tACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAAcc                                  <  1:5088391/39‑1 (MQ=255)                                                                 tACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAAcc                                  <  1:3399252/39‑1 (MQ=255)                                                                 tACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAAcc                                  <  1:3078938/39‑1 (MQ=255)                                                                 tACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGAt  <  1:2213120/71‑1 (MQ=255)                                                                  |                                                                      GTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGAT  >  Exported/609105‑609239

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑