Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,854	C→A	7.9%	I253I (ATC→ATA)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,854	0	C	A	7.9%	244.6 / 4.3	75	I253I (ATC→ATA)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base A (4/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

TACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGAT  >  minE/2570785‑2570923                                                                      |                                                                      tACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATc                                                                       >  1:2892101/1‑70 (MQ=255)   cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATcccc                                                                    <  1:2316874/71‑1 (MQ=255)   cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                                     <  1:3493113/70‑1 (MQ=255)   cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                                     <  1:315582/70‑1 (MQ=255)   cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                                     <  1:1413914/70‑1 (MQ=255)     cAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGa                                                                  >  1:3094513/1‑71 (MQ=255)      aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGat                                                                 >  1:1626206/1‑70 (MQ=255)      aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGat                                                                 >  1:1277992/1‑70 (MQ=255)      aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg                                                                 >  1:3257161/1‑71 (MQ=255)        gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa                                                               <  1:2080786/71‑1 (MQ=255)         tACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaa                                                              <  1:3020262/71‑1 (MQ=255)           cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                                                            <  1:2419635/71‑1 (MQ=255)           cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                             <  1:70292/70‑1 (MQ=255)              aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACg                                                         <  1:660131/71‑1 (MQ=255)               tCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCGAAAACgg                                                        <  1:3952017/71‑1 (MQ=255)                   cccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                                     <  1:1866897/54‑1 (MQ=255)                   cccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCaaa                                                    <  1:3843591/71‑1 (MQ=255)                   cccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCaa                                                     <  1:2235528/70‑1 (MQ=255)                    ccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAt                                                   >  1:2956086/1‑71 (MQ=255)                      aGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTATTCATCCCCGAGCAAAAACGGCAAATTg                                                 <  1:1975554/71‑1 (MQ=255)                        gTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTc                                               >  1:2704787/1‑71 (MQ=255)                          aCATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGc                                             >  1:2167269/1‑71 (MQ=255)                              aCGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCaa                                         >  1:1104172/1‑71 (MQ=255)                               cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa                                                               >  1:1540175/1‑48 (MQ=255)                               cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                        >  1:1085065/1‑71 (MQ=255)                               cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                        >  1:1774530/1‑71 (MQ=255)                                gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCaa                                                     <  1:2368150/57‑1 (MQ=255)                                gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAtt                                       <  1:3594239/71‑1 (MQ=255)                                gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAtt                                       <  1:3456558/71‑1 (MQ=255)                                gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAtt                                       <  1:722732/71‑1 (MQ=255)                                 cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                      <  1:2847195/71‑1 (MQ=255)                                  aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                      <  1:2036316/70‑1 (MQ=255)                                  aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTc                                     >  1:3330861/1‑71 (MQ=255)                                     tGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaaa                                                           >  1:1228381/1‑46 (MQ=255)                                       ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaaa                                                           >  1:1194273/1‑44 (MQ=255)                                       ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                                                            <  1:3066398/43‑1 (MQ=255)                                       ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCa                                 >  1:2607976/1‑70 (MQ=255)                                        cTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaa                               <  1:1956051/71‑1 (MQ=255)                                        cTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                >  1:1227510/1‑70 (MQ=255)                                        cTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                >  1:3238906/1‑70 (MQ=255)                                         ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc                                                          <  1:3686817/43‑1 (MQ=255)                                         ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                              >  1:2131217/1‑71 (MQ=255)                                         ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                              >  1:369339/1‑71 (MQ=255)                                          tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAAc                             >  1:3042317/1‑71 (MQ=255)                                          tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAAc                             >  1:1613731/1‑71 (MQ=255)                                           gCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACt                            >  1:1021290/1‑71 (MQ=255)                                           gCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACt                            >  1:958001/1‑71 (MQ=255)                                            ccGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAtt                                                  >  1:1523492/1‑48 (MQ=255)                                            ccGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTc                           <  1:1738223/71‑1 (MQ=255)                                            ccGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTc                           >  1:2138665/1‑71 (MQ=255)                                            ccGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACTc                           <  1:2482413/71‑1 (MQ=255)                                            acGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAtt                                                  >  1:983289/2‑48 (MQ=255)                                             cGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCg                          >  1:441072/1‑71 (MQ=255)                                             cGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACTcc                          <  1:2263165/71‑2 (MQ=255)                                               cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAt                                                   >  1:1356220/1‑44 (MQ=255)                                               cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCt                        >  1:3393766/1‑71 (MQ=255)                                               cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCt                        >  1:1415037/1‑71 (MQ=255)                                                 gCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                        <  1:1683833/53‑1 (MQ=255)                                                  cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTg                       <  1:1516537/69‑1 (MQ=255)                                                  cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGt                      >  1:1658581/1‑70 (MQ=255)                                                  cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGt                      >  1:1653882/1‑70 (MQ=255)                                                    tCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                   <  1:255181/71‑1 (MQ=255)                                                    tCCGGTGATGGTCTTCATCCCCGAGAAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                   >  1:2941404/1‑71 (MQ=255)                                                     ccGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGc                                             >  1:3330836/1‑44 (MQ=255)                                                     ccGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                   <  1:3767516/70‑1 (MQ=255)                                                      cGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt                 <  1:3473696/71‑1 (MQ=255)                                                      cGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt                 <  1:2072660/71‑1 (MQ=255)                                                      aggTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGt                      <  1:2417400/65‑1 (MQ=255)                                                        gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                >  1:2052476/1‑54 (MQ=255)                                                        gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt               >  1:1661046/1‑71 (MQ=255)                                                        gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt               >  1:972126/1‑71 (MQ=255)                                                          gATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc                <  1:1366071/68‑1 (MQ=255)                                                          gATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc                <  1:3912398/68‑1 (MQ=255)                                                              gTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAcc         >  1:2632932/1‑71 (MQ=255)                                                               tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt        <  1:1480792/71‑1 (MQ=255)                                                               tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt        <  1:388802/71‑1 (MQ=255)                                                               tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt        >  1:1281247/1‑71 (MQ=255)                                                                  tCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTgg      >  1:2152602/1‑70 (MQ=255)                                                                    aTCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCga   >  1:3344326/1‑71 (MQ=255)                                                                     tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat  <  1:2523396/71‑1 (MQ=255)                                                                     tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat  <  1:3304425/71‑1 (MQ=255)                                                                     accccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt                 >  1:2165273/2‑56 (MQ=255)                                                                      |                                                                      TACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGAT  >  minE/2570785‑2570923

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