Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,397	G→A	18.7%	N619N (AAC→AAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,397	0	G	A	18.7%	136.2 / 22.6	64	N619N (AAC→AAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base A (4/8);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGG  >  Exported/608331‑608457                                                                   |                                                             gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgtgt                                                             <  1:502969/68‑1 (MQ=255) gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                                                          >  1:5106630/1‑71 (MQ=255) gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                                                          >  1:3947601/1‑71 (MQ=255) gTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAc                                                          <  1:1035538/71‑1 (MQ=255)    aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccgg                                                       >  1:4037640/1‑69 (MQ=255)    aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccg                                                        <  1:2079948/70‑2 (MQ=255)      gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                                                     <  1:3841163/71‑1 (MQ=255)      gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                                                     <  1:1515024/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:5220220/1‑71 (MQ=255)                        cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                   <  1:3273572/71‑1 (MQ=255)                         aCCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                  <  1:1479695/71‑1 (MQ=255)                          cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                  >  1:2916497/1‑70 (MQ=255)                            cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                               <  1:1196942/71‑1 (MQ=255)                            cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                               <  1:4386550/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:3983295/49‑1 (MQ=255)                                         gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                                        <  1:2340756/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:1599268/1‑71 (MQ=255)                                               gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa            >  1:1605881/1‑71 (MQ=255)                                               gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa            >  1:2419006/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:985849/71‑1 (MQ=255)                                                  cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg         <  1:1131484/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:3858501/53‑1 (MQ=255)                                                             tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                <  1:3540171/53‑1 (MQ=255)                                                              ggggTGTTACCCGTTTTCTTACCATACACAACGTTAg                               <  1:1368028/37‑1 (MQ=255)                                                                   |                                                             GTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGG  >  Exported/608331‑608457

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