Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,454	C→A	10.3%	P600P (CCG→CCT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,454	0	C	A	10.3%	241.0 / 9.0	77	P600P (CCG→CCT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base A (4/4);  total (36/41)
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.40e-01

GACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCA  >  Exported/608384‑608517                                                                       |                                                                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:858281/71‑1 (MQ=255)    cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                                                              <  1:3959433/71‑1 (MQ=255)    cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                                                              <  1:3416414/71‑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:1062633/1‑66 (MQ=255)                     gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                                  >  1:3634530/1‑66 (MQ=255)                     gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                             >  1:4710746/1‑71 (MQ=255)                     gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                             >  1:2541054/1‑71 (MQ=255)                     gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                             >  1:2879597/1‑71 (MQ=255)                     gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGgt                                             >  1:2102100/1‑71 (MQ=255)                      ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCAc                                                    <  1:4317281/63‑1 (MQ=255)                        ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgtgt                                           <  1:4831545/70‑1 (MQ=255)                        ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg                                          >  1:1910693/1‑71 (MQ=255)                        ttCTTACCATACACAACGTTAGAAGAGATGGTCAGAACAGACTGAGTCGGGATAg                                                          >  1:3934166/1‑55 (MQ=255)                             aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCa                                        <  1:3568046/68‑1 (MQ=255)                             aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt                                      >  1:1920407/1‑70 (MQ=255)                             aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt                                      <  1:4449503/70‑1 (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:711807/1‑42 (MQ=255)                                  acacAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGa                                                             >  1:5409775/1‑42 (MQ=255)                                  acacAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGa                                                             >  1:2274774/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:1923857/1‑71 (MQ=255)                                            aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa                      >  1:2664794/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:2557183/70‑1 (MQ=255)                                                   aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                <  1:5169601/70‑1 (MQ=255)                                                   aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg               >  1:4177682/1‑71 (MQ=255)                                                   aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg               <  1:411506/71‑1 (MQ=255)                                                      gTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGGGTGCAGTTTCTGAATTTTCTTCATGAAACGtt             <  1:470384/70‑1 (MQ=255)                                                       tCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCt           <  1:836843/71‑1 (MQ=255)                                                        cAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                            >  1:3275993/1‑53 (MQ=255)                                                        cAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                            >  1:512016/1‑53 (MQ=255)                                                         aGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTa          >  1:457559/1‑70 (MQ=255)                                                          gAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTAcc        <  1:4320318/71‑1 (MQ=255)                                                           aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                <  1:727378/62‑1 (MQ=255)                                                           aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                <  1:3001829/62‑1 (MQ=255)                                                           aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                <  1:4152402/62‑1 (MQ=255)                                                            aCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAg                                       >  1:5564359/1‑38 (MQ=255)                                                            aCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTg                                 >  1:2980204/1‑44 (MQ=255)                                                             cAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAgg     >  1:2523699/1‑71 (MQ=255)                                                                aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                <  1:1551903/57‑1 (MQ=255)                                                                aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCa  <  1:5237242/71‑1 (MQ=255)                                                                aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCa  <  1:712001/71‑1 (MQ=255)                                                                aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGATCTACCAGGTc   >  1:2196096/1‑70 (MQ=255)                                                                 cTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                            <  1:4296603/44‑1 (MQ=255)                                                                 cTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                            <  1:4267765/44‑1 (MQ=255)                                                                  tGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                            <  1:1389460/43‑1 (MQ=255)                                                                      tCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTc   <  1:2804514/64‑1 (MQ=255)                                                                       cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                            <  1:5526022/38‑1 (MQ=255)                                                                       cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGt              <  1:1309733/52‑1 (MQ=255)                                                                       cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGt              <  1:3884306/52‑1 (MQ=255)                                                                       |                                                                GACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCA  >  Exported/608384‑608517

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