Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,093,811	A→G	7.0%	C347C (TGT→TGC)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,093,811	0	A	G	7.0%	159.2 / 2.0	57	C347C (TGT→TGC)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base A (33/20);  new base G (2/2);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATA  >  Exported/2093744‑2093876                                                                    |                                                                   ccgccgCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACa                                                                  <  1:161852/70‑1 (MQ=255)  cgccgcGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCacg                                                                >  1:4570746/1‑69 (MQ=255)         ttGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                                          >  1:4963297/1‑70 (MQ=255)          tGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCAcgccacgcca                                                        <  1:543722/71‑1 (MQ=255)           gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAAc                                                                   >  1:1229003/1‑59 (MQ=255)           gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                                        <  1:4668797/70‑1 (MQ=255)           gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                                        <  1:4730031/70‑1 (MQ=255)               aCCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCACGCCACGCCaactt                                                    >  1:3027797/1‑67 (MQ=255)                 cTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑ccc                                                >  1:4067880/1‑71 (MQ=255)                    cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAtc                                                     >  1:5441428/1‑64 (MQ=255)                    cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCt                                             <  1:4555033/71‑1 (MQ=255)                    cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCt                                             <  1:756842/71‑1 (MQ=255)                      ccGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGa                                           <  1:663082/71‑1 (MQ=255)                      ccGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGa                                           <  1:344475/71‑1 (MQ=255)                        ggACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa                                                                    >  1:3686379/1‑45 (MQ=255)                        ggACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACa                                         <  1:1509830/71‑1 (MQ=255)                         gACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAg                                        >  1:258567/1‑71 (MQ=255)                             gccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCg                                     >  1:2350007/1‑70 (MQ=255)                              ccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGa                                    >  1:2655589/1‑70 (MQ=255)                              ccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGa                                    <  1:384699/70‑1 (MQ=255)                                  ccAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTcc                               <  1:995537/71‑1 (MQ=255)                                   cAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                                          >  1:552915/1‑44 (MQ=255)                                   cAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                                          >  1:3114831/1‑44 (MQ=255)                                       agtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                                        <  1:4597969/42‑1 (MQ=255)                                        gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCAc                                                            >  1:3106350/1‑37 (MQ=255)                                        gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCAc                                                            >  1:998413/1‑37 (MQ=255)                                        gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCAc                                                            >  1:1670102/1‑37 (MQ=255)                                        gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAt                                                      >  1:756109/1‑43 (MQ=255)                                        gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAt                                                      >  1:5150474/1‑43 (MQ=255)                                        gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCgcgc                           >  1:3478369/1‑69 (MQ=255)                                        gtTCAGTTAACACTGCCGCCGCCATTGAACAGGCCAc                                                            >  1:5085282/1‑37 (MQ=255)                                         ttCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGa                                    <  1:2961925/59‑1 (MQ=255)                                         ttCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTc                                <  1:1246004/63‑1 (MQ=255)                                          tCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACa                                         <  1:4327679/53‑1 (MQ=255)                                           cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCaga                       <  1:155535/70‑1 (MQ=255)                                           cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCa                         >  1:4659274/1‑68 (MQ=255)                                           cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCa                         >  1:1812151/1‑68 (MQ=255)                                           cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCa                         >  1:5442082/1‑68 (MQ=255)                                           cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCa                         >  1:4780405/1‑68 (MQ=255)                                                aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCa                         >  1:2918784/1‑63 (MQ=255)                                                aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAgaga                     >  1:4206147/1‑67 (MQ=255)                                                aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAgaga                     >  1:3914397/1‑67 (MQ=255)                                                aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATgg                  <  1:2754508/70‑1 (MQ=255)                                                aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATg                   >  1:2734480/1‑69 (MQ=255)                                                   ccTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGc              <  1:4921248/71‑1 (MQ=255)                                                       tcgccgccATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGtt           >  1:4804364/2‑70 (MQ=255)                                                       ccgccgccATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGtt           <  1:1938807/70‑1 (MQ=255)                                                        cgccgccATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCa         <  1:2045170/71‑1 (MQ=255)                                                         gccgccATTGAGCACGCCACGCCAACTT‑CACCCTGGCATCCCACTTCAGCACCaga                       <  1:5632257/56‑1 (MQ=255)                                                         gccgccATTGAACAGGCCACGCCAATCTCCCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCa         >  1:5512852/1‑71 (MQ=255)                                                               aTTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGAt                    >  1:241526/1‑53 (MQ=255)                                                                ttGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATg                   <  1:921595/53‑1 (MQ=255)                                                                ttGAACAGGCCACGCCAATCT‑CACCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTtata  >  1:1819504/1‑70 (MQ=255)                                                                 tGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAgaga                     >  1:689356/1‑50 (MQ=255)                                                                  gAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCaga                       >  1:2958131/1‑47 (MQ=255)                                                                  gAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATGGa                 >  1:2173376/1‑53 (MQ=255)                                                                  gAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCAttt      >  1:1605944/1‑64 (MQ=255)                                                                  gAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTtata  >  1:5608419/1‑68 (MQ=255)                                                                    aCAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAgaga                     >  1:3301418/1‑47 (MQ=255)                                                                    |                                                                   CCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCT‑CCCCCTGACAGCCGACTTCCGCGCCAGAGATGGAGGCGTTCATTTTATA  >  Exported/2093744‑2093876

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