Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	876,263	G→A	29.9%	G493G (GGG→GGA)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	876,263	0	G	A	29.9%	79.6 / 35.0	80	G493G (GGG→GGA)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (17/39);  new base A (17/7);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 1.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

GGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGA  >  minE/876198‑876331                                                                  |                                                                     ggTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGttt                                                                   <  1:3303645/69‑1 (MQ=255) ggTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGttt                                                                   <  1:2268351/69‑1 (MQ=255) ggTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGttt                                                                   <  1:2552332/69‑1 (MQ=255) ggTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGttt                                                                   <  1:721667/69‑1 (MQ=255)           tCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTgg                                                         <  1:412675/69‑1 (MQ=255)           tCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTgg                                                         <  1:3073527/69‑1 (MQ=255)           tCGTTTATGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTgg                                                         <  1:1327335/69‑1 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACTCCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:1132208/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTggg                                                        >  1:460025/1‑68 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:1615814/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:1594954/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:2802219/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:446739/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:1724466/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:1778759/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:2311240/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:2050867/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:2154143/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:1010506/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:1159566/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:156093/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:3114323/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:944064/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGATGAACTGGGATTTAACCTTGTGGGa                                                       >  1:304942/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGa                                                       >  1:2882007/1‑69 (MQ=255)             gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGa                                                       >  1:824851/1‑69 (MQ=255)                      tATCTGGCAAAAGCGAAACTTACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGt                                              <  1:484118/69‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGt                                                 <  1:922244/66‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGt                                                 <  1:1377784/66‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGt                                                 <  1:1311753/66‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGt                                                 <  1:1974412/66‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGt                                                 <  1:467665/66‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGt                                                 <  1:2637498/66‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGt                                              <  1:3011809/69‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGt                                              <  1:3175654/69‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGt                                              <  1:1152599/69‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGt                                              <  1:121897/69‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGt                                              <  1:1903493/69‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGt                                              <  1:165820/69‑1 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGG‑TTAACCTTGTGGGa                                                       >  1:1680618/1‑59 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGG‑TTAACCTTGTGGGATACGGTTg                                               >  1:1322761/1‑67 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGG‑TTAACCTTGTGGGATACGGTTGt                                              >  1:2047967/1‑68 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGG‑TTAACCTTGTGGGATACGGTTGt                                              >  1:1458512/1‑68 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGG‑TTAACCTTGTGGGATACGGTTGt                                              >  1:167346/1‑68 (MQ=255)                      tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGG‑TTAACCTTGTGGGATACGGTTGt                                              >  1:2442836/1‑68 (MQ=255)                      tATCAGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGG‑TTAACCTTGTGGGATACGGTTGt                                              >  1:1725104/1‑68 (MQ=255)                       aTCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGt                                              <  1:1910633/68‑1 (MQ=255)                        tCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGt                                                 <  1:2399914/64‑1 (MQ=255)                              aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTg                                      <  1:2774411/69‑1 (MQ=255)                              aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTg                                      <  1:2246190/69‑1 (MQ=255)                              aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTg                                      <  1:2359411/69‑1 (MQ=255)                              aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTg                                      <  1:1008022/69‑1 (MQ=255)                              aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTg                                      <  1:2593231/69‑1 (MQ=255)                              aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTg                                      <  1:2583051/69‑1 (MQ=255)                              aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTg                                      <  1:2494610/69‑1 (MQ=255)                               aaaGCGAAACTGACACCGTATCTCGACGAACTGGGGGTTAACCTTGTGGGATACGGTTGTACCACCTGt                                     <  1:3418912/69‑1 (MQ=255)                                aaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGt                                     <  1:1592279/68‑1 (MQ=255)                                  gCGAAACTGACACCGTATCTCGACGAACTGGGGttt                                                                   <  1:2611971/36‑1 (MQ=255)                                        cTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAAc                            >  1:764086/1‑69 (MQ=255)                                        cTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAAc                            >  1:1428313/1‑69 (MQ=255)                                                 tATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTg                                      <  1:1672745/50‑1 (MQ=255)                                                     tcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGc               >  1:2385924/1‑69 (MQ=255)                                                     tcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGc               >  1:3154823/1‑69 (MQ=255)                                                     tcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGc               >  1:3055483/1‑69 (MQ=255)                                                     tcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGc               >  1:176844/1‑69 (MQ=255)                                                     tcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGc               >  1:216869/1‑69 (MQ=255)                                                     tcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGc               >  1:2153826/1‑69 (MQ=255)                                                         tgaACTGGGATTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGACCGTTGccc             <  1:508288/66‑1 (MQ=255)                                                         tgaACTGGGATTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGACCGTTGccc             <  1:1390605/66‑1 (MQ=255)                                                         tgaACTGGGATTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGACCGTTGccc             <  1:1523588/66‑1 (MQ=255)                                                         tgaACTGGGATTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGACCGTTGccc             <  1:770327/66‑1 (MQ=255)                                                         tgaACTGGGATTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGACCGTTGccc             <  1:3160903/66‑1 (MQ=255)                                                         tgaACTGGGATTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGACCGTTGccc             <  1:2459569/66‑1 (MQ=255)                                                         tgaACTGGGATTTAACCTTGTGGGATACGGTTGTACCACCTGTAATGGTAACTCTGGACCGTTGccc             <  1:96470/66‑1 (MQ=255)                                                             cTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCt       <  1:584074/69‑1 (MQ=255)                                                             cTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCt       <  1:2339486/69‑1 (MQ=255)                                                             cTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCt       <  1:89091/69‑1 (MQ=255)                                                             cTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCt       <  1:2589979/69‑1 (MQ=255)                                                             cTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCt       <  1:1252034/69‑1 (MQ=255)                                                              tGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCt       <  1:1620482/68‑1 (MQ=255)                                                                  gTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGa  <  1:492048/69‑1 (MQ=255)                                                                  gTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGa  <  1:454357/69‑1 (MQ=255)                                                                  |                                                                     GGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGA  >  minE/876198‑876331

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