Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,391	T→C	20.3%	P621P (CCA→CCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,391	0	T	C	20.3%	132.3 / 28.7	69	P621P (CCA→CCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base C (5/9);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGG  >  Exported/608321‑608457                                                                       |                                                                   cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                                                    >  1:5464753/1‑71 (MQ=255) cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                                                    >  1:5213759/1‑71 (MQ=255)   ccTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCttg                                                                  <  1:2552236/71‑3 (MQ=255)      ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                                               >  1:2646935/1‑71 (MQ=255)      ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                                               >  1:4138635/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)                        aTCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTgg                                                                  >  1:2202266/1‑50 (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:2862864/64‑1 (MQ=255)                                             aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                               <  1:5010659/64‑1 (MQ=255)                                             aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGAt                        >  1:3351370/1‑71 (MQ=255)                                               cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATgg                      >  1:3361145/1‑71 (MQ=255)                                                 gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                     >  1:1669343/1‑70 (MQ=255)                                                 gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                     >  1:2058523/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:2419006/1‑71 (MQ=255)                                                         gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa            >  1:1599268/1‑71 (MQ=255)                                                         gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa            >  1:1605881/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: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)                                                                       |                                                                   CACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGG  >  Exported/608321‑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: ‑