Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,010	C→T	10.1%	Q748Q (CAG→CAA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,010	0	C	T	10.1%	212.0 / 11.2	79	Q748Q (CAG→CAA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (43/27);  new base T (3/5);  total (47/32)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGA  >  Exported/607940‑608076                                                                       |                                                                    cTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgc                                                                     >  1:5226863/1‑71 (MQ=255)   ttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctg                                                                   >  1:2236666/1‑71 (MQ=255)    tttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgc                                                                  >  1:4449022/1‑71 (MQ=255)      tACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgct‑g                                                               >  1:3002028/1‑71 (MQ=255)         tGTACTATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctg                                                                   <  1:970631/65‑1 (MQ=255)           tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑Gtt                                                             >  1:3642488/1‑68 (MQ=255)           tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCtt                                                          >  1:2897703/1‑71 (MQ=255)              gATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAg                                                       <  1:4346096/71‑1 (MQ=255)               aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑Gt                                                              >  1:1604128/1‑63 (MQ=255)               aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑Gt                                                              >  1:1914988/1‑63 (MQ=255)               aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑Gt                                                              >  1:1146529/1‑63 (MQ=255)                tttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTc                                                     <  1:355191/71‑1 (MQ=255)                 ttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgct                                                                 >  1:2868653/1‑59 (MQ=255)                 ttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgct                                                                 >  1:3999315/1‑59 (MQ=255)                    aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGc                                                  <  1:2762980/70‑1 (MQ=255)                    aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGc                                                  <  1:2413541/70‑1 (MQ=255)                        aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                             <  1:645509/71‑1 (MQ=255)                        aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                             >  1:2591745/1‑71 (MQ=255)                        aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                             >  1:229397/1‑71 (MQ=255)                        aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                             >  1:4352840/1‑71 (MQ=255)                        aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                             >  1:973397/1‑71 (MQ=255)                        aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                             <  1:2151974/71‑1 (MQ=255)                         aaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCtt                                                          <  1:894144/57‑1 (MQ=255)                           tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGa                                                <  1:1309261/65‑1 (MQ=255)                            cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                         >  1:2170650/1‑71 (MQ=255)                            cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                         >  1:4477486/1‑71 (MQ=255)                            cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                         >  1:982748/1‑71 (MQ=255)                               aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                         <  1:4419698/68‑1 (MQ=255)                               aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                       >  1:2298041/1‑70 (MQ=255)                               aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                       <  1:1962975/70‑1 (MQ=255)                               aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                      <  1:3026259/71‑1 (MQ=255)                                ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                      >  1:4396549/1‑70 (MQ=255)                                ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                     >  1:2300988/1‑71 (MQ=255)                                ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGATGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                      >  1:3494523/1‑70 (MQ=255)                                   cATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACt                                   <  1:4815956/70‑1 (MQ=255)                                    aTAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTa                                                        <  1:3725039/48‑1 (MQ=255)                                    aTAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTa                                                        <  1:756690/48‑1 (MQ=255)                                    aTAGATTGAGTGAAGGTACGAGTAATAAAGTCCTGCTGCT‑GTTCTTTa                                                        <  1:4201361/48‑1 (MQ=255)                                     tAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGGCAGCGAGTTGAAACGTACTGc                                 <  1:3749529/70‑1 (MQ=255)                                      aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                              >  1:71531/1‑56 (MQ=255)                                      aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                              >  1:3197352/1‑56 (MQ=255)                                      aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                              >  1:176854/1‑56 (MQ=255)                                      aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                              >  1:210154/1‑56 (MQ=255)                                      aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                              >  1:2244981/1‑56 (MQ=255)                                      aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                              >  1:2338224/1‑56 (MQ=255)                                      aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                              >  1:2561674/1‑56 (MQ=255)                                      aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCg                                >  1:3536577/1‑70 (MQ=255)                                         ttGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                            <  1:2617571/71‑1 (MQ=255)                                          tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                       >  1:1801633/1‑59 (MQ=255)                                          tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                       >  1:1953076/1‑59 (MQ=255)                                           gAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGAGAGTTGAAACGTACTGCGTAGcc                           <  1:83408/70‑1 (MQ=255)                                             gtaaaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGa                        <  1:350587/68‑1 (MQ=255)                                              tGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                      <  1:1294776/56‑1 (MQ=255)                                                aaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                      <  1:404675/70‑3 (MQ=255)                                                aaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                      <  1:2969174/70‑3 (MQ=255)                                                aaGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                      >  1:5315911/1‑70 (MQ=255)                                                aaGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                      >  1:227011/1‑70 (MQ=255)                                                 aGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                       >  1:805418/1‑52 (MQ=255)                                                 aGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                       >  1:5010636/1‑52 (MQ=255)                                                    tACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac                   >  1:3279435/1‑69 (MQ=255)                                                    tACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                  >  1:398507/1‑70 (MQ=255)                                                        aGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACTGATgg             <  1:87811/71‑1 (MQ=255)                                                        aGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg             <  1:3814349/71‑1 (MQ=255)                                                          taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                 <  1:3326983/65‑1 (MQ=255)                                                          taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                 <  1:5158368/65‑1 (MQ=255)                                                          taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                  <  1:366473/64‑1 (MQ=255)                                                           aataaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGCATTGTCa          <  1:569736/71‑1 (MQ=255)                                                            ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt               >  1:2715948/1‑65 (MQ=255)                                                            ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt               >  1:3593189/1‑65 (MQ=255)                                                            ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt               >  1:531529/1‑65 (MQ=255)                                                            ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt               >  1:5538430/1‑65 (MQ=255)                                                            ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg         <  1:472346/71‑1 (MQ=255)                                                            ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg         >  1:1064665/1‑71 (MQ=255)                                                            ataaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAg         >  1:4040379/1‑71 (MQ=255)                                                             taaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa          <  1:3450180/69‑1 (MQ=255)                                                              aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                     >  1:2849456/1‑41 (MQ=255)                                                              aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                     >  1:739180/1‑41 (MQ=255)                                                              aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                     >  1:3235411/1‑41 (MQ=255)                                                              aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAg                             <  1:1961296/49‑1 (MQ=255)                                                               aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGt                               >  1:1115134/1‑46 (MQ=255)                                                               aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                 >  1:2975706/1‑60 (MQ=255)                                                               aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa          >  1:2340301/1‑67 (MQ=255)                                                                 gTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg    >  1:5527442/1‑71 (MQ=255)                                                                 gTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg    >  1:3563115/1‑71 (MQ=255)                                                                  tCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt            <  1:5225139/62‑1 (MQ=255)                                                                   ccTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTAc                                    >  1:2986568/1‑37 (MQ=255)                                                                    ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa          >  1:753132/1‑62 (MQ=255)                                                                    ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa  >  1:3375138/1‑70 (MQ=255)                                                                    ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa  >  1:2966344/1‑70 (MQ=255)                                                                    ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa  >  1:1178530/1‑70 (MQ=255)                                                                    ctgctgca‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa  >  1:562899/1‑70 (MQ=255)                                                                     tgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt       <  1:1232043/64‑1 (MQ=255)                                                                       ctgctgGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                 >  1:4867945/1‑37 (MQ=38)                                                                       |                                                                    CTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGA  >  Exported/607940‑608076

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