Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,833	G→A	8.4%	L246L (CTG→CTA)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,833	0	G	A	8.4%	246.9 / 9.5	83	L246L (CTG→CTA)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/39);  new base A (4/3);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCT  >  minE/2570765‑2570901                                                                     |                                                                     gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCg                                                                    <  1:3418881/71‑1 (MQ=255) gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCg                                                                    <  1:43089/71‑1 (MQ=255)   gcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                                                  >  1:1557229/1‑71 (MQ=255)    ccGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                                                  <  1:2240903/70‑1 (MQ=255)      gAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                                                               >  1:3315257/1‑71 (MQ=255)       aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCg                                                                    <  1:2927951/65‑1 (MQ=255)       aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGt                                                              <  1:3936123/71‑1 (MQ=255)       aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGt                                                              <  1:358038/71‑1 (MQ=255)        gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGt                                                                   >  1:2357414/1‑65 (MQ=255)         cAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGa                                                            <  1:1659244/71‑1 (MQ=255)         cAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGa                                                            <  1:220481/71‑1 (MQ=255)           aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATg                                                          <  1:2915519/71‑1 (MQ=255)           aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATg                                                          >  1:895186/1‑71 (MQ=255)            tGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGAt                                                           <  1:1699239/69‑1 (MQ=255)             ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATgg                                                         >  1:418068/1‑70 (MQ=255)             ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt                                                        >  1:3459218/1‑71 (MQ=255)              gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                                                       >  1:3838290/1‑71 (MQ=255)               cTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCt                                                      <  1:2534870/71‑1 (MQ=255)                  cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCa                                                   >  1:3808903/1‑71 (MQ=255)                   gTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                    <  1:632607/69‑1 (MQ=255)                     tACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATc                                                 >  1:2892101/1‑70 (MQ=255)                       cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATcccc                                              <  1:2316874/71‑1 (MQ=255)                       cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                               <  1:1413914/70‑1 (MQ=255)                       cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                               <  1:315582/70‑1 (MQ=255)                       cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                               <  1:3493113/70‑1 (MQ=255)                        tCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGTTGCGTCCGGTGAt                                                           <  1:2296213/57‑1 (MQ=255)                        tCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGAt                                                           <  1:2331357/57‑1 (MQ=255)                        tCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt                                                        <  1:2933217/60‑1 (MQ=255)                         cAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGa                                            >  1:3094513/1‑71 (MQ=255)                          aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGat                                           >  1:1277992/1‑70 (MQ=255)                          aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGat                                           >  1:1626206/1‑70 (MQ=255)                          aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg                                           >  1:3257161/1‑71 (MQ=255)                            gTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGtgt                                                      <  1:289122/58‑3 (MQ=255)                            gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                    >  1:2573784/1‑60 (MQ=255)                            gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa                                         <  1:2080786/71‑1 (MQ=255)                             tACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaa                                        <  1:3020262/71‑1 (MQ=255)                              aCTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGc                                                                     >  1:2165715/1‑41 (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:722732/71‑1 (MQ=255)                                                    gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAtt                 <  1:3456558/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:3238906/1‑70 (MQ=255)                                                            cTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa          >  1:1227510/1‑70 (MQ=255)                                                             ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc                                    <  1:3686817/43‑1 (MQ=255)                                                             ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa        >  1:369339/1‑71 (MQ=255)                                                             ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa        >  1:2131217/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:2138665/1‑71 (MQ=255)                                                                ccGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTc     <  1:1738223/71‑1 (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:1415037/1‑71 (MQ=255)                                                                   cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCt  >  1:3393766/1‑71 (MQ=255)                                                                     gCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                  <  1:1683833/53‑1 (MQ=255)                                                                     |                                                                     GCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCT  >  minE/2570765‑2570901

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