Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,295	A→C	26.1%	A653A (GCT→GCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,295	0	A	C	26.1%	102.5 / 64.1	88	A653A (GCT→GCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (39/26);  new base C (16/7);  total (55/33)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

ACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGG  >  Exported/608226‑608360                                                                      |                                                                  aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                                   <  1:5430160/70‑1 (MQ=255) aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                                   <  1:4161336/70‑1 (MQ=255)      gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                             >  1:5103641/1‑71 (MQ=255)       tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                            <  1:666583/71‑1 (MQ=255)       tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                             >  1:3162435/1‑70 (MQ=255)       tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                             >  1:2595823/1‑70 (MQ=255)         cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                            >  1:726501/1‑69 (MQ=255)         cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                          <  1:184186/71‑1 (MQ=255)              ggAACGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                                      >  1:3755348/1‑70 (MQ=255)              ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                     <  1:5653632/71‑1 (MQ=255)              ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                     <  1:1390292/71‑1 (MQ=255)                aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                      >  1:4607131/1‑68 (MQ=255)                 aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                    >  1:513774/1‑67 (MQ=255)                 aCGATTGAGAAGGTGTACGAGATCCCATATTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                    >  1:4234740/1‑67 (MQ=255)                   gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                          <  1:4060079/61‑1 (MQ=255)                   gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                          <  1:5255274/61‑1 (MQ=255)                    aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                               <  1:1327759/71‑1 (MQ=255)                     ttgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                               >  1:5678493/3‑70 (MQ=255)                     tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                               >  1:3348759/1‑70 (MQ=255)                     tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                               >  1:1762916/1‑70 (MQ=255)                     tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                               <  1:4135293/70‑1 (MQ=255)                      agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                             >  1:4378990/1‑71 (MQ=255)                      agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                             >  1:5573347/1‑71 (MQ=255)                      agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                             >  1:2969125/1‑71 (MQ=255)                       gagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGc                                                                 <  1:4579883/50‑1 (MQ=255)                         gaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                          <  1:920718/71‑1 (MQ=255)                          aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                         <  1:3026919/71‑1 (MQ=255)                          aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                         <  1:3033685/71‑1 (MQ=255)                            ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                              >  1:1318659/1‑48 (MQ=255)                            ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                              >  1:3370390/1‑48 (MQ=255)                            ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                              >  1:3013444/1‑48 (MQ=255)                            ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                              >  1:1242477/1‑48 (MQ=255)                             gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                               <  1:15115/62‑1 (MQ=255)                             gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCAcc                                      >  1:397857/1‑71 (MQ=255)                              tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                      <  1:5068583/54‑1 (MQ=255)                              tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                      <  1:2288181/54‑1 (MQ=255)                               gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                                    >  1:267150/1‑71 (MQ=255)                                tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                    >  1:324205/1‑54 (MQ=255)                                 aGGAGATACCATCTTAAGCGAAAGCAAACGGCAGTTTAGc                                                                 <  1:177471/40‑1 (MQ=38)                                  ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                                   >  1:1595621/1‑37 (MQ=255)                                  ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                  >  1:5311233/1‑54 (MQ=255)                                  ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                 >  1:415944/1‑71 (MQ=255)                                  cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                                      <  1:813430/49‑1 (MQ=255)                                   gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                >  1:4114886/1‑68 (MQ=255)                                   gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                >  1:1693457/1‑68 (MQ=255)                                   gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                <  1:834402/71‑4 (MQ=255)                                   gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                             >  1:4826911/1‑58 (MQ=255)                                   gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                             >  1:3596066/1‑58 (MQ=255)                                      aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                                        <  1:3959556/44‑1 (MQ=255)                                        aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                               >  1:498555/1‑35 (MQ=255)                                        aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                               >  1:2515299/1‑35 (MQ=255)                                        aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                               >  1:2090346/1‑35 (MQ=255)                                        aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGa                                                           >  1:585966/1‑39 (MQ=255)                                        aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                           >  1:5691076/1‑71 (MQ=255)                                         ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                           >  1:3762863/1‑70 (MQ=255)                                         ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                           >  1:3321072/1‑70 (MQ=255)                                         ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                               <  1:1633976/66‑1 (MQ=255)                                         ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                           >  1:526394/1‑70 (MQ=255)                                          cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                         >  1:489970/1‑71 (MQ=255)                                          cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                         >  1:196632/1‑71 (MQ=255)                                          cATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                    <  1:3901048/44‑1 (MQ=255)                                           aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                        <  1:1258203/71‑1 (MQ=255)                                            tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGATCACGAc                        >  1:432341/1‑70 (MQ=255)                                             cTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAcca                      >  1:5093733/1‑70 (MQ=255)                                              tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                    >  1:4230043/1‑40 (MQ=255)                                              tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                    <  1:1066885/40‑1 (MQ=255)                                              tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                     >  1:3917764/1‑39 (MQ=255)                                               ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                     >  1:3836571/1‑70 (MQ=255)                                               ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                    >  1:1192396/1‑71 (MQ=255)                                               ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCCTg                    >  1:4201025/1‑71 (MQ=255)                                                tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                       <  1:635675/67‑1 (MQ=255)                                                tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                     <  1:2513253/69‑1 (MQ=255)                                                tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                    >  1:3385579/1‑70 (MQ=255)                                                   cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                <  1:3595044/71‑1 (MQ=255)                                                    gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg               >  1:5471197/1‑71 (MQ=255)                                                     tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg               <  1:2509113/70‑1 (MQ=255)                                                       aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                        >  1:4945308/1‑59 (MQ=255)                                                       aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                        >  1:3140718/1‑59 (MQ=255)                                                       aGCAAACGGCAGTTTAGAAACGGAAGTCAGAGAGGCTACTCCACCTTTCTGGTAACGACCGTGCATCGGGt            >  1:4693906/1‑71 (MQ=255)                                                           aaCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGc        >  1:2058665/1‑71 (MQ=255)                                                            aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg       <  1:945225/71‑2 (MQ=255)                                                            aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg       >  1:2844125/1‑70 (MQ=255)                                                            aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg       >  1:3138578/1‑70 (MQ=255)                                                            aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg       >  1:321985/1‑70 (MQ=255)                                                            aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                              >  1:2937546/1‑48 (MQ=255)                                                            aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGc        >  1:1711072/1‑70 (MQ=255)                                                              ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                   <  1:2386775/57‑1 (MQ=255)                                                              ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                   <  1:3407023/57‑1 (MQ=255)                                                               gCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGccc    <  1:1432672/71‑1 (MQ=255)                                                                 aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg  <  1:3073683/71‑1 (MQ=255)                                                                 aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg  <  1:1762056/71‑1 (MQ=255)                                                                      |                                                                  ACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGG  >  Exported/608226‑608360

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