Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,212,854	T→A	5.5%	I5N (ATT→AAT)	yobH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,212,854	0	T	A	5.5%	259.8 / 5.7	90	I5N (ATT→AAT)	yobH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/51);  new base A (3/2);  total (37/53)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TGCTTCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGTTTTGA  >  minE/1212784‑1212919                                                                       |                                                                  tgcttCAGCGCATATTGACGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                   >  1:3156554/1‑71 (MQ=255)  gcttCAGCGCATATTGGCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                   <  1:3277791/70‑1 (MQ=255)  gcttCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                   <  1:2389046/70‑1 (MQ=255)  gcttCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                   <  1:2525969/70‑1 (MQ=255)      cAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGca                                                              >  1:523720/1‑71 (MQ=255)      cAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGca                                                              >  1:1701447/1‑71 (MQ=255)       aGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGcac                                                             <  1:2620396/71‑1 (MQ=255)       aGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGca                                                              >  1:1169397/1‑70 (MQ=255)       aGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGca                                                              >  1:3447645/1‑70 (MQ=255)       aGCGCATATTGCCGCTAAGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGAATGATGAATCGcat                                                             <  1:1251357/71‑2 (MQ=255)        gcgcATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGAAATGCGATTCATCATTCgaa                                                              >  1:678781/1‑67 (MQ=255)             tatTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTAt                                                       <  1:1491297/71‑1 (MQ=255)              atTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTAt                                                       <  1:3568419/70‑1 (MQ=255)                tGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                    >  1:3044916/1‑71 (MQ=255)                 gCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                    <  1:478818/70‑1 (MQ=255)                 gCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                    <  1:429803/70‑1 (MQ=255)                   cGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGAt                                                 <  1:443524/71‑1 (MQ=255)                   cGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGAt                                                 <  1:896662/71‑1 (MQ=255)                    gCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                    <  1:1584533/67‑1 (MQ=255)                    gCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGAt                                                 <  1:1345729/70‑1 (MQ=255)                        cGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCAc                                            <  1:713195/71‑1 (MQ=255)                            tAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         >  1:3144332/1‑70 (MQ=255)                            tAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         >  1:2579033/1‑70 (MQ=255)                            tAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         >  1:2463037/1‑70 (MQ=255)                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                   <  1:18404/43‑1 (MQ=255)                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                        <  1:3463506/70‑1 (MQ=255)                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                        <  1:1521084/70‑1 (MQ=255)                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                        <  1:2026097/70‑1 (MQ=255)                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                        <  1:171884/70‑1 (MQ=255)                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                        <  1:3490180/70‑1 (MQ=255)                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                        <  1:126628/70‑1 (MQ=255)                             aaGCGAACATAAAAACAGAAGAGGTTGTAATGCGATGCATCATTCGCACAGTTATGCTGATAGCACTGGt                                        <  1:2440365/70‑1 (MQ=255)                              aGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTa                                       <  1:2730639/70‑1 (MQ=255)                              aGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTa                                       <  1:681767/70‑1 (MQ=255)                                  aaCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGAt                                  <  1:375387/71‑1 (MQ=255)                                   aCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTAt                                                       >  1:1311326/1‑49 (MQ=255)                                   aCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTAt                                                       >  1:934189/1‑49 (MQ=255)                                   aCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGa                                   >  1:888625/1‑69 (MQ=255)                                   aCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGa                                   >  1:2195428/1‑69 (MQ=255)                                    cATAAAAAGAGAAGATGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAg                                <  1:2618588/71‑1 (MQ=255)                                    cATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAg                                >  1:705242/1‑71 (MQ=255)                                      tAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGa                              <  1:571527/71‑1 (MQ=255)                                       aaaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                             >  1:2184047/1‑71 (MQ=255)                                        aaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                             >  1:528763/1‑70 (MQ=255)                                        aaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                             >  1:40880/1‑70 (MQ=255)                                        aaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                             >  1:111163/1‑70 (MQ=255)                                        aaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                             >  1:1843003/1‑70 (MQ=255)                                          aagagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGattat                          >  1:2751613/1‑71 (MQ=255)                                          aagagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGattat                          <  1:2496561/71‑1 (MQ=255)                                          aagagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGattat                          <  1:1842050/71‑1 (MQ=255)                                          aagagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTAATAGCACTGGTATGGATAGtattat                          <  1:52290/71‑1 (MQ=255)                                          aagagaagagGTTGTAATGCGAATGATGAATCGcat                                                             >  1:559331/1‑35 (MQ=9)                                          aagagaagagGTTGTAATGCGAATGATGAATCGcat                                                             >  1:2090826/1‑35 (MQ=9)                                          aagagaagagGTTGTAATGCGAATGATGAATCGca                                                              >  1:2551730/1‑35 (MQ=12)                                           agagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGattatt                         <  1:659146/71‑1 (MQ=255)                                            gagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGAttatta                        <  1:2558826/71‑1 (MQ=255)                                             agaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTAc                       >  1:2550229/1‑71 (MQ=255)                                              gaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                    <  1:978167/41‑1 (MQ=255)                                               aagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACt                      <  1:3134944/70‑1 (MQ=255)                                                agagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTTGTATGGATAGGattat                          >  1:1774031/1‑65 (MQ=255)                                                agagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTTGTATGGATAGGattat                          >  1:1383582/1‑65 (MQ=255)                                                agagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTTGTATGGATAGGattat                          >  1:791758/1‑65 (MQ=255)                                                agagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACt                      >  1:3594056/1‑69 (MQ=255)                                                 gagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGatt                            <  1:2105977/62‑1 (MQ=255)                                                 gagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGatt                            <  1:2101532/62‑1 (MQ=255)                                                 gagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGatt                            <  1:106768/62‑1 (MQ=255)                                                 gagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAg                   <  1:728624/71‑1 (MQ=255)                                                   ggTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCg                 <  1:2855628/71‑1 (MQ=255)                                                       gTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGata             >  1:907658/1‑71 (MQ=255)                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATTGATAGGATTATTACTCAGCGGatat            <  1:2131645/71‑1 (MQ=255)                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGa                                   <  1:981172/48‑1 (MQ=255)                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGa                                   <  1:1879091/48‑1 (MQ=255)                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGa                                   <  1:3235325/48‑1 (MQ=255)                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGatat            <  1:3319188/71‑1 (MQ=255)                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGatat            <  1:1886657/71‑1 (MQ=255)                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGatat            <  1:1933167/71‑1 (MQ=255)                                                         aaTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGatat            <  1:87869/70‑1 (MQ=255)                                                          gtGCGAATGATGAATCGCACAGTTATGCTGATAGCACTGGTATg                                     <  1:3482005/43‑1 (MQ=18)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         >  1:1356725/1‑40 (MQ=255)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         >  1:618864/1‑40 (MQ=255)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         >  1:809168/1‑40 (MQ=255)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         >  1:2079543/1‑40 (MQ=255)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         >  1:2024988/1‑40 (MQ=255)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTg                                          >  1:2711318/1‑39 (MQ=255)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATg           <  1:936376/70‑1 (MQ=255)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATg           <  1:1840887/70‑1 (MQ=255)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATg           <  1:418561/70‑1 (MQ=255)                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATg           <  1:2201186/70‑1 (MQ=255)                                                           tGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         <  1:935766/39‑1 (MQ=255)                                                            gCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                         <  1:3285753/38‑1 (MQ=255)                                                            gCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCg        >  1:526194/1‑71 (MQ=255)                                                             cGATTGATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGt       <  1:1054421/71‑1 (MQ=255)                                                               aTTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGttt     <  1:692002/71‑1 (MQ=255)                                                               aTTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGttt     <  1:2265151/71‑1 (MQ=255)                                                                ttCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGtttt    >  1:2748453/1‑71 (MQ=255)                                                                    tcatTCGCACAGTTATGCTGATAGCACTGTTATGGATAGGATTATTACTCAGCGGat              <  1:3572123/57‑1 (MQ=255)                                                                      atTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATgg          >  1:3637227/1‑59 (MQ=255)                                                                       ttCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGTTTTGa  <  1:3045441/66‑1 (MQ=255)                                                                       |                                                                  TGCTTCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGTTTTGA  >  minE/1212784‑1212919

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