Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,298	A→C	25.9%	V652V (GTT→GTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,298	0	A	C	25.9%	103.6 / 61.5	85	V652V (GTT→GTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (37/26);  new base C (16/6);  total (53/32)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAA  >  Exported/608231‑608366                                                                    |                                                                     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:1390292/71‑1 (MQ=255)         ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                           <  1:5653632/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:5255274/61‑1 (MQ=255)              gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                                <  1:4060079/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)                    gaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                <  1:920718/71‑1 (MQ=255)                     aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                               <  1:3033685/71‑1 (MQ=255)                     aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                               <  1:3026919/71‑1 (MQ=255)                       ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                                    >  1:3370390/1‑48 (MQ=255)                       ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                                    >  1:1242477/1‑48 (MQ=255)                       ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                                    >  1:1318659/1‑48 (MQ=255)                       ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                                    >  1:3013444/1‑48 (MQ=255)                        gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                     <  1:15115/62‑1 (MQ=255)                        gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCAcc                                            >  1:397857/1‑71 (MQ=255)                         tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                            <  1:2288181/54‑1 (MQ=255)                         tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                            <  1:5068583/54‑1 (MQ=255)                          gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                                          >  1:267150/1‑71 (MQ=255)                           tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                          >  1:324205/1‑54 (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:2090346/1‑35 (MQ=255)                                   aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                                     >  1:498555/1‑35 (MQ=255)                                   aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                                     >  1:2515299/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:196632/1‑71 (MQ=255)                                     cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                               >  1:489970/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)                                                                  gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa  <  1:2120626/71‑1 (MQ=255)                                                                   cAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                 >  1:1451238/1‑39 (MQ=255)                                                                    |                                                                     GTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAA  >  Exported/608231‑608366

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