Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,073	C→A	24.8%	P727P (CCG→CCT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,073	0	C	A	24.8%	156.0 / 27.9	64	P727P (CCG→CCT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (27/21);  new base A (7/9);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCT  >  minE/608004‑608143                                                                      |                                                                       gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                                       >  1:5527442/1‑71 (MQ=255) gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                                       >  1:3563115/1‑71 (MQ=255)    ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                                     >  1:2966344/1‑70 (MQ=255)    ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                                     >  1:3375138/1‑70 (MQ=255)    ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                                     >  1:1178530/1‑70 (MQ=255)    ctgctgcagTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                                     >  1:562899/1‑70 (MQ=255)          ctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                                                              <  1:1119205/71‑1 (MQ=255)             ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                                            <  1:1488244/70‑1 (MQ=255)              tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                                          >  1:5385737/1‑71 (MQ=255)              tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                                          >  1:2144533/1‑71 (MQ=255)              tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                                          >  1:3291109/1‑71 (MQ=255)                   aGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                                                              <  1:5332446/62‑1 (MQ=255)                    gTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCgg                                                          >  1:3904211/1‑65 (MQ=9)                    gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc                                                    >  1:4713787/1‑71 (MQ=255)                    gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc                                                    <  1:2381112/71‑1 (MQ=255)                     tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc                                                   <  1:4341197/71‑1 (MQ=255)                     tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc                                                   <  1:1067665/71‑1 (MQ=255)                       aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                                                       <  1:5478999/65‑1 (MQ=255)                       aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                                                       <  1:862411/65‑1 (MQ=255)                       aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                                                 >  1:4465790/1‑71 (MQ=255)                         cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATcg                                               >  1:2599153/1‑71 (MQ=255)                           cgTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTcc                                                            <  1:5500462/55‑1 (MQ=12)                             ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                                                             <  1:1169838/53‑1 (MQ=255)                             ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                                                             <  1:1855915/53‑1 (MQ=255)                               gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCa                                                  >  1:3090896/1‑62 (MQ=255)                               gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                          <  1:1841906/70‑1 (MQ=255)                                 aaCGTACTGCGTAGCCAGATACAGGGATGGTCAGCTGCGGATATTTTTcc                                                            <  1:4784597/50‑1 (MQ=255)                                  aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                                      <  1:2574704/71‑1 (MQ=255)                                  aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                        >  1:4723853/1‑69 (MQ=255)                                     tACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttt                                                                <  1:1550897/42‑1 (MQ=255)                                       cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                                 >  1:2496566/3‑71 (MQ=12)                                       cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAt                                  >  1:2145640/3‑70 (MQ=12)                                         gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                                            <  1:4827450/42‑1 (MQ=255)                                         gCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTc                               >  1:3894100/1‑71 (MQ=16)                                          cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                                                 >  1:4987806/1‑52 (MQ=255)                                          cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                <  1:2534990/69‑1 (MQ=255)                                          cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                              >  1:686846/1‑71 (MQ=255)                                             aGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga                           >  1:611247/1‑70 (MQ=12)                                              gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCgat                          >  1:204410/1‑69 (MQ=12)                                                 agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtc                        <  1:3476824/70‑1 (MQ=9)                                                 aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                       >  1:4707475/1‑71 (MQ=255)                                                  gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                      <  1:1890052/68‑1 (MQ=9)                                                  gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                      <  1:1208217/68‑1 (MQ=9)                                                  gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                      <  1:3047823/68‑1 (MQ=9)                                                    tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                     <  1:4459557/70‑1 (MQ=255)                                                    tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                     <  1:3986411/70‑1 (MQ=255)                                                    gacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcatta                    >  1:2583400/2‑69 (MQ=9)                                                     acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtc                        >  1:2752349/1‑66 (MQ=255)                                                     acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                   >  1:125761/1‑71 (MQ=255)                                                     acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                   >  1:378263/1‑71 (MQ=255)                                                      cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                  <  1:3994210/71‑1 (MQ=12)                                                      cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                  <  1:2878476/71‑1 (MQ=12)                                                      cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                  <  1:5515929/71‑1 (MQ=12)                                                      cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                  <  1:2610403/71‑1 (MQ=12)                                                      cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACg                  >  1:1900190/1‑71 (MQ=255)                                                        cGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                <  1:724340/71‑1 (MQ=255)                                                         ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                >  1:2256425/1‑70 (MQ=255)                                                          gATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa              >  1:3215964/1‑71 (MQ=255)                                                              gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa              <  1:388756/67‑1 (MQ=255)                                                              gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGt           <  1:623050/70‑1 (MQ=255)                                                                cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                     >  1:3268764/1‑58 (MQ=255)                                                                cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt          >  1:672858/1‑69 (MQ=255)                                                                 aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa              >  1:5346504/1‑64 (MQ=255)                                                                 aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg       >  1:4238926/1‑71 (MQ=255)                                                                      cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg    <  1:3561424/69‑1 (MQ=255)                                                                      cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt  >  1:4815554/1‑71 (MQ=255)                                                                      cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt  <  1:2347595/71‑1 (MQ=255)                                                                      |                                                                       GTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCT  >  minE/608004‑608143

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑