Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,193,222	T→A	47.2%	V147V (GTA→GTT)	fadD ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,193,222	0	T	A	47.2%	1.9 / 82.2	70	V147V (GTA→GTT)	fadD	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base A (21/12);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

TTTAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACA  >  minE/1193156‑1193288                                                                   |                                                                   tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTGAGAATTAc                                                                  >  1:1930767/1‑69 (MQ=255) tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:381919/1‑69 (MQ=255) tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:857753/1‑69 (MQ=255) tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:1294466/1‑69 (MQ=255) tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:581883/1‑69 (MQ=255) tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:2067728/1‑69 (MQ=255)       aacGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTg                                                            <  1:1013237/69‑1 (MQ=255)       aacGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTg                                                            <  1:2703162/69‑1 (MQ=255)       aacGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTg                                                            <  1:678085/69‑1 (MQ=255)                   tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:61011/69‑1 (MQ=255)                   tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:2068553/69‑1 (MQ=255)                   tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:199135/69‑1 (MQ=255)                   tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:1528251/69‑1 (MQ=255)                   tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGcttt                                                <  1:1399077/69‑1 (MQ=255)                    aCCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGttt                                                 <  1:1060670/67‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:2479076/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:3387999/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:3200109/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:2904346/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:506101/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:1123599/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:2526621/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:826259/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:1755477/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:1873045/68‑1 (MQ=255)                                   gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGGTTTATCAACCACTTTTTc                                 <  1:63189/68‑1 (MQ=255)                                      gATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTTTCAACCACTTTTTCCAgt                             >  1:131118/1‑69 (MQ=255)                                      gATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgt                             >  1:1465220/1‑69 (MQ=255)                                      gATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgt                             >  1:2117789/1‑69 (MQ=255)                                      gATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgt                             >  1:1846808/1‑69 (MQ=255)                                         aGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgtgt                           >  1:1376675/1‑68 (MQ=255)                                         aGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgtgt                           >  1:388671/1‑68 (MQ=255)                                                   ccATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:305981/37‑1 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:3034641/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2990387/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:1620330/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:3102201/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:805130/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:3280145/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2896607/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:285904/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2778637/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2716363/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2712581/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:266609/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2424192/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:508315/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2413383/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2285532/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:614062/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:1325561/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:85307/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAAACACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:409850/1‑68 (MQ=255)                                                      tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAAACACTTTTTCCAGAGTGCGGGCAAAGTTa              >  1:2500694/1‑68 (MQ=255)                                                       aCGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:992221/48‑1 (MQ=255)                                                        cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:3164954/47‑1 (MQ=255)                                                        cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:1005772/47‑1 (MQ=255)                                                        cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:988995/47‑1 (MQ=255)                                                        cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:683618/47‑1 (MQ=255)                                                        cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:602527/47‑1 (MQ=255)                                                        cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:2322562/47‑1 (MQ=255)                                                        cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:450390/47‑1 (MQ=255)                                                        cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:3306108/47‑1 (MQ=255)                                                        cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGac           >  1:2180596/1‑69 (MQ=255)                                                             cAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:662781/1‑69 (MQ=255)                                                             cAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:2976133/1‑69 (MQ=255)                                                              aGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCTAAGTTAGACACGAt      >  1:935892/1‑68 (MQ=255)                                                              aGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:2423991/1‑68 (MQ=255)                                                              aGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:2130187/1‑68 (MQ=255)                                                              aGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:1833096/1‑68 (MQ=255)                                                                  ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:2711638/68‑1 (MQ=255)                                                                  ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1839383/68‑1 (MQ=255)                                                                  ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1736746/68‑1 (MQ=255)                                                                  ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1629760/68‑1 (MQ=255)                                                                  ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1358054/68‑1 (MQ=255)                                                                  ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1185488/68‑1 (MQ=255)                                                                  ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:3135452/68‑1 (MQ=255)                                                                   |                                                                   TTTAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACA  >  minE/1193156‑1193288

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