Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,436	G→A	14.4%	T606T (ACC→ACT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,436	0	G	A	14.4%	200.1 / 21.6	83	T606T (ACC→ACT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (38/33);  new base A (6/6);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

CGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACG  >  Exported/608367‑608504                                                                      |                                                                     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)      cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                                                                 >  1:2490241/1‑70 (MQ=255)       gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                                               >  1:781163/1‑71 (MQ=255)           gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                                           >  1:1605881/1‑71 (MQ=255)           gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                                           >  1:2419006/1‑71 (MQ=255)           gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                                           >  1:1599268/1‑71 (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)               gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACt                                                         >  1:779657/1‑69 (MQ=255)                 cgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGt                                                     >  1:1411661/1‑71 (MQ=255)                  gacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTc                                                    <  1:2463970/71‑1 (MQ=255)                   acCGGCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCg                                                   <  1:998088/71‑1 (MQ=255)                    ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg                                                  <  1:1327464/71‑1 (MQ=255)                    ccGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                                                  >  1:3422674/1‑71 (MQ=255)                     cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                                                 <  1:3416414/71‑1 (MQ=255)                     cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                                                 <  1:858281/71‑1 (MQ=255)                     cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                                                 <  1:3959433/71‑1 (MQ=255)                      gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                             >  1:5010318/1‑58 (MQ=255)                      gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                             >  1:1857874/1‑58 (MQ=255)                         tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                                               <  1:3540171/53‑1 (MQ=255)                         tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                                               <  1:3858501/53‑1 (MQ=255)                                ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGta                                                    >  1:2900815/1‑56 (MQ=255)                                ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                                               <  1:1265460/46‑1 (MQ=255)                                ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg                                      >  1:4274323/1‑71 (MQ=255)                                ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGCAt                                               <  1:4533783/62‑1 (MQ=255)                                 tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                                                  <  1:4226166/58‑1 (MQ=255)                                 tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGt                                          <  1:3955736/66‑1 (MQ=255)                                  aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                     >  1:115970/1‑70 (MQ=255)                                      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                     >  1:3634530/1‑66 (MQ=255)                                      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                     >  1:1062633/1‑66 (MQ=255)                                      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                >  1:2879597/1‑71 (MQ=255)                                      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                >  1:4710746/1‑71 (MQ=255)                                      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                >  1:2541054/1‑71 (MQ=255)                                      gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGgt                                >  1:2102100/1‑71 (MQ=255)                                       ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                             <  1:5667625/41‑1 (MQ=255)                                       ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                             <  1:38625/41‑1 (MQ=255)                                       ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                                             <  1:2378405/41‑1 (MQ=255)                                       ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCAc                                       <  1:4317281/63‑1 (MQ=255)                                        tttCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACt                                                         >  1:4245046/1‑44 (MQ=12)                                         ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgtgt                              <  1:4831545/70‑1 (MQ=255)                                         ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg                             >  1:1910693/1‑71 (MQ=255)                                         ttCTTACCATACACAACGTTAGAAGAGATGGTCAGAACAGACTGAGTCGGGATAg                                             >  1:3934166/1‑55 (MQ=255)                                          tCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                                           <  1:2128709/40‑1 (MQ=255)                                              aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCa                           <  1:3568046/68‑1 (MQ=255)                                              aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt                         <  1:4449503/70‑1 (MQ=255)                                              aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt                         >  1:1920407/1‑70 (MQ=255)                                               ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATATCGTCACGGTAGGTGTGCAGttt                       <  1:82175/71‑1 (MQ=255)                                               ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAg                          <  1:5049511/68‑1 (MQ=255)                                                 aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGt                                    <  1:3255532/56‑1 (MQ=255)                                                 aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCt                     >  1:1650878/1‑71 (MQ=255)                                                   acacAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGa                                                >  1:2274774/1‑42 (MQ=255)                                                   acacAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGa                                                >  1:711807/1‑42 (MQ=255)                                                   acacAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGa                                                >  1:5409775/1‑42 (MQ=255)                                                    caccacGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAgg                                 >  1:471017/1‑56 (MQ=255)                                                    cacaACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCa                                        >  1:5395537/1‑49 (MQ=255)                                                    cacaACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGaa                  >  1:2005261/1‑71 (MQ=255)                                                       aaCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt               >  1:1961911/1‑71 (MQ=255)                                                       aaCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt               >  1:932886/1‑71 (MQ=255)                                                       aaCGTTAGAAGTGATGGTCAGAACAGACTCAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt               <  1:1590691/71‑1 (MQ=255)                                                        aCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt               >  1:5092046/1‑70 (MQ=255)                                                         cGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGGTTCTGAATTttc             <  1:4513366/71‑1 (MQ=255)                                                           ttAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt               <  1:478686/67‑1 (MQ=255)                                                           ttAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt               <  1:1665726/67‑1 (MQ=255)                                                           ttAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTttctt           <  1:3191103/71‑1 (MQ=255)                                                             aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa         >  1:2664794/1‑71 (MQ=255)                                                             aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa         >  1:1923857/1‑71 (MQ=255)                                                             aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa         >  1:4144598/1‑71 (MQ=255)                                                             aGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCa         >  1:3691980/1‑71 (MQ=255)                                                              gAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCAt        <  1:3737269/71‑1 (MQ=255)                                                              gAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTGAATTTTCTTCAt        <  1:903448/71‑1 (MQ=255)                                                                  tGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTttc             <  1:1842762/62‑1 (MQ=255)                                                                  tGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGa      <  1:5430284/69‑1 (MQ=255)                                                                  tGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGa      <  1:2555051/69‑1 (MQ=255)                                                                   gATGGTCAGAACAGACTGAGTCGGGATAGAGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc   >  1:1607965/1‑71 (MQ=255)                                                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAg                                  >  1:1619658/1‑39 (MQ=255)                                                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc   <  1:5169601/70‑1 (MQ=255)                                                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc   <  1:2557183/70‑1 (MQ=255)                                                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg  >  1:4177682/1‑71 (MQ=255)                                                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg  <  1:411506/71‑1 (MQ=255)                                                                      |                                                                     CGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACG  >  Exported/608367‑608504

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