Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,403	C→G	20.8%	T617T (ACG→ACC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,403	0	C	G	20.8%	149.4 / 32.1	72	T617T (ACG→ACC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base G (6/9);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGG  >  Exported/608336‑608469                                                                    |                                                                   gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                                                                 <  1:1515024/71‑1 (MQ=255) gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                                                                 <  1:3841163/71‑1 (MQ=255)  aCCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGttt                                                                <  1:4754005/71‑1 (MQ=255)     gTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCt                                                             <  1:2869956/71‑1 (MQ=255)      tGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                                                                 >  1:5672208/1‑66 (MQ=255)       gCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTa                                                           >  1:4342200/1‑71 (MQ=255)        cATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCtggc                                                          <  1:1104124/71‑4 (MQ=255)            gggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATa                                                      <  1:1332168/71‑1 (MQ=255)                 aGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACaa                                                 >  1:4628740/1‑71 (MQ=255)                   cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                               <  1:3273572/71‑1 (MQ=255)                   cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                               >  1:5220220/1‑71 (MQ=255)                    aCCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                              <  1:1479695/71‑1 (MQ=255)                     cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                              >  1:2916497/1‑70 (MQ=255)                       cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                           <  1:4386550/71‑1 (MQ=255)                       cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                           <  1:1196942/71‑1 (MQ=255)                         gTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGttt                                                                <  1:5369907/48‑1 (MQ=255)                           ccGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAcc                                                         >  1:321984/1‑53 (MQ=255)                            cGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                       >  1:3633108/1‑70 (MQ=255)                            cGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                      <  1:899285/71‑1 (MQ=255)                              aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                           <  1:5010659/64‑1 (MQ=255)                              aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                           <  1:2862864/64‑1 (MQ=255)                              aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGAt                                    >  1:3351370/1‑71 (MQ=255)                                cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATgg                                  >  1:3361145/1‑71 (MQ=255)                                  gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                 >  1:2058523/1‑70 (MQ=255)                                  gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                 >  1:1669343/1‑70 (MQ=255)                                  gcgcgcCAGCACGACGACCGTCTGGGGTGTTAACCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                 >  1:2981590/1‑70 (MQ=255)                                    gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                                                    <  1:2340756/49‑1 (MQ=255)                                    gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                                                    <  1:3983295/49‑1 (MQ=255)                                    gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                              <  1:210484/55‑1 (MQ=255)                                    gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                              <  1:4115035/55‑1 (MQ=255)                                     cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                              >  1:2490241/1‑70 (MQ=255)                                      gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                            >  1:781163/1‑71 (MQ=255)                                          gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                                   <  1:1390336/44‑1 (MQ=255)                                          gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                      <  1:392881/57‑1 (MQ=255)                                          gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                        >  1:2419006/1‑71 (MQ=255)                                          gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                        >  1:1605881/1‑71 (MQ=255)                                          gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                        >  1:1599268/1‑71 (MQ=255)                                           cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                                   <  1:3516012/43‑1 (MQ=255)                                           cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCa                               >  1:2961006/1‑63 (MQ=255)                                           cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGAc                       <  1:5056468/71‑1 (MQ=255)                                             cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                     <  1:1131484/71‑1 (MQ=255)                                             cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                     <  1:985849/71‑1 (MQ=255)                                              gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTATaa                                         >  1:964193/1‑50 (MQ=255)                                              gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACt                      >  1:779657/1‑69 (MQ=255)                                               acgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                           <  1:5523749/47‑1 (MQ=255)                                                cgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGt                  >  1:1411661/1‑71 (MQ=255)                                                 gacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTc                 <  1:2463970/71‑1 (MQ=255)                                                  acCGGCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCg                <  1:998088/71‑1 (MQ=255)                                                   ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTa                                            >  1:2208333/1‑42 (MQ=255)                                                   ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg               <  1:1327464/71‑1 (MQ=255)                                                   ccGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg               >  1:3422674/1‑71 (MQ=255)                                                    cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg              <  1:3959433/71‑1 (MQ=255)                                                    cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg              <  1:3416414/71‑1 (MQ=255)                                                    cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg              <  1:858281/71‑1 (MQ=255)                                                     gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                          >  1:5010318/1‑58 (MQ=255)                                                     gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                          >  1:1857874/1‑58 (MQ=255)                                                        tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                            <  1:3540171/53‑1 (MQ=255)                                                        tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                            <  1:3858501/53‑1 (MQ=255)                                                         ggggTGTTACCCGTTTTCTTACCATACACAACGTTAg                                           <  1:1368028/37‑1 (MQ=255)                                                               ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                       >  1:5545369/1‑35 (MQ=12)                                                               ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGa                                     <  1:18036/37‑1 (MQ=12)                                                               ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGta                 >  1:2900815/1‑56 (MQ=255)                                                               ttACCCGTTTTCTTACCATACACAACGTTAGAAGt                                       <  1:3801679/35‑1 (MQ=255)                                                               ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                            <  1:1265460/46‑1 (MQ=255)                                                               ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg   >  1:4274323/1‑71 (MQ=255)                                                               ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGCAt            <  1:4533783/62‑1 (MQ=255)                                                                tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg               <  1:4226166/58‑1 (MQ=255)                                                                tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGt       <  1:3955736/66‑1 (MQ=255)                                                                 aCCGGTTTTCTGGCCGTACACCACGTTAGAAGTGAta                                   >  1:4473590/1‑36 (MQ=12)                                                                 aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                     >  1:3691654/1‑35 (MQ=255)                                                                 aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                     >  1:1242288/1‑35 (MQ=255)                                                                 aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg  >  1:115970/1‑70 (MQ=255)                                                                    |                                                                   GACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGG  >  Exported/608336‑608469

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