Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,338	C→T	14.8%	V306I (GTT→ATT) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,338	0	C	T	14.8%	174.6 / 19.3	74	V306I (GTT→ATT) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (35/28);  new base T (3/8);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGA  >  minE/609270‑609406                                                                     |                                                                     gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa                                                                    <  1:1885619/71‑1 (MQ=255)   aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt                                                                  >  1:1719359/1‑71 (MQ=255)   aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt                                                                  >  1:4089161/1‑71 (MQ=255)    cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  >  1:2999037/1‑70 (MQ=255)         cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCGACCATTTCCtac                                                            >  1:3684893/1‑69 (MQ=255)         cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                                            >  1:1981534/1‑71 (MQ=255)         cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                                            >  1:3693472/1‑71 (MQ=255)            caaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgc                                                          >  1:1560493/4‑70 (MQ=255)             gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGctt                                                        >  1:4107251/1‑71 (MQ=255)             aaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcc                                                         <  1:939224/68‑2 (MQ=255)                ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGcttctt                                                     >  1:5283728/1‑71 (MQ=255)                 tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  <  1:4064226/57‑1 (MQ=255)                 tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  <  1:5318468/57‑1 (MQ=255)                   cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                  <  1:4578514/71‑1 (MQ=255)                   cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                  <  1:236672/71‑1 (MQ=255)                     gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                             >  1:4909761/3‑58 (MQ=255)                     gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                             >  1:1729511/3‑58 (MQ=255)                     gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAttt                                                                 <  1:4530147/54‑1 (MQ=255)                        tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:4251922/60‑3 (MQ=255)                        tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:1301603/60‑3 (MQ=255)                        tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:339874/60‑3 (MQ=255)                        tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:4567928/60‑3 (MQ=255)                        tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:1205687/60‑3 (MQ=255)                        tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:2205378/60‑3 (MQ=255)                        tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:3122764/60‑3 (MQ=255)                          cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:299515/1‑50 (MQ=255)                          cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:1787213/1‑50 (MQ=255)                             aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:2102882/1‑47 (MQ=255)                             aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:1682521/1‑47 (MQ=255)                             aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:4362599/1‑47 (MQ=255)                              ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  <  1:242077/44‑1 (MQ=255)                              ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  <  1:811417/44‑1 (MQ=255)                              ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                       >  1:299446/1‑71 (MQ=255)                              ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                       >  1:1523225/1‑71 (MQ=255)                              ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                       <  1:3169799/71‑1 (MQ=255)                               gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCt                                         >  1:4128862/1‑68 (MQ=255)                               gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                                      >  1:5601422/1‑71 (MQ=255)                                 aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                       >  1:5283690/1‑68 (MQ=255)                                 aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                                      >  1:2070553/1‑69 (MQ=255)                                    gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATc                                          <  1:4481714/62‑1 (MQ=255)                                    gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                                 >  1:990197/1‑71 (MQ=255)                                    gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                                 <  1:1970934/71‑1 (MQ=255)                                       ccATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:3614193/1‑37 (MQ=255)                                         aTACGCAGTTTCTTGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                  <  1:30087/49‑1 (MQ=255)                                            cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc                         <  1:4365751/71‑1 (MQ=255)                                            cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc                         <  1:297722/71‑1 (MQ=255)                                                gTTTCATGACCAGGTGGTCTACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGt                     <  1:5111201/71‑1 (MQ=255)                                                 tttCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCttt                               <  1:960013/60‑1 (MQ=255)                                                   tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                                        >  1:4065374/1‑49 (MQ=255)                                                   tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                                        >  1:422164/1‑49 (MQ=255)                                                   tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   >  1:728656/1‑70 (MQ=255)                                                   tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:333195/70‑1 (MQ=255)                                                   tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:4254341/70‑1 (MQ=255)                                                   tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:1026504/70‑1 (MQ=255)                                                   tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   >  1:3027650/1‑70 (MQ=255)                                                   tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:1665705/70‑1 (MQ=255)                                                   tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:2058404/70‑1 (MQ=255)                                                      tGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATg               >  1:178792/1‑71 (MQ=255)                                                          cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                             >  1:2169634/1‑53 (MQ=255)                                                          cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg                            >  1:1121064/1‑54 (MQ=255)                                                          cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg                            >  1:2332258/1‑54 (MQ=255)                                                            ggtggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                             >  1:680557/1‑51 (MQ=255)                                                             gtggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACAt        >  1:5601343/1‑71 (MQ=255)                                                               ggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCa                                    >  1:2663644/1‑41 (MQ=255)                                                                gtCAACCATTTCCTGCTCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATcc      <  1:1769834/70‑1 (MQ=255)                                                                gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATc       >  1:2417452/1‑69 (MQ=255)                                                                gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATc       >  1:4725418/1‑69 (MQ=255)                                                                gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCa     <  1:4826755/71‑1 (MQ=255)                                                                  cAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTa             <  1:5674027/61‑1 (MQ=255)                                                                   aaCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAgg   >  1:4826953/1‑70 (MQ=255)                                                                   aaCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAgg   <  1:3425056/70‑1 (MQ=255)                                                                   aaCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGa  <  1:1648544/71‑1 (MQ=255)                                                                   aaCCATTTACTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGa  <  1:2019610/71‑1 (MQ=255)                                                                    aCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGGCAGCTTTCAGGTCACGTTc                   <  1:275840/53‑1 (MQ=255)                                                                     |                                                                     GAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGA  >  minE/609270‑609406

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