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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,394	0	G	C	20.3%	137.3 / 30.2	69	T620T (ACC→ACG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base C (5/9);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

TTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGG  >  Exported/608326‑608457                                                                     |                                                                ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                                               >  1:4138635/1‑71 (MQ=255) ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                                               >  1:2646935/1‑71 (MQ=255)  tgcggTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCgt                                                               <  1:1831786/67‑1 (MQ=255)  tCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                                               <  1:2347125/70‑1 (MQ=255)      gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgtgt                                                             <  1:502969/68‑1 (MQ=255)      gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                                                          >  1:3947601/1‑71 (MQ=255)      gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                                                          >  1:5106630/1‑71 (MQ=255)      gTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAc                                                          <  1:1035538/71‑1 (MQ=255)       tCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCgta                                                              >  1:4677054/1‑65 (MQ=255)         aCGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTggg                                                                 >  1:928333/1‑61 (MQ=255)         aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccgg                                                       >  1:4037640/1‑69 (MQ=255)         aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccg                                                        <  1:2079948/70‑2 (MQ=255)           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: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: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:1605881/1‑71 (MQ=255)                                                    gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa            >  1:1599268/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: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:858281/71‑1 (MQ=255)                                                              cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg  <  1:3416414/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)                                                                     |                                                                TTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGG  >  Exported/608326‑608457

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