Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	530,638	T→C	53.8%	E105E (GAA→GAG)	ybhP ←	predicted DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	530,638	0	T	C	53.8%	22.9 / 126.9	93	E105E (GAA→GAG)	ybhP	predicted DNase
Reads supporting (aligned to +/- strand):  ref base T (13/30);  new base C (26/24);  total (39/54)
Fisher's exact test for biased strand distribution p-value = 3.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGCGGCAGTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTGCGACCGT  >  minE/530572‑530705                                                                   |                                                                    aaTGCGGCAGTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCaa                                                                 <  1:2355729/71‑1 (MQ=255)     cGGCAGTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATgg                                                              >  1:2774231/1‑70 (MQ=255)     cGGCAGTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATgg                                                              >  1:1088677/1‑70 (MQ=255)      ggCAGTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGa                                                            >  1:2275830/1‑71 (MQ=255)         aGTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCa                                                                  >  1:2569141/1‑62 (MQ=255)          gTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAAc                                                        <  1:481442/71‑1 (MQ=255)          gTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAAc                                                        <  1:1025640/71‑1 (MQ=255)          gTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAAc                                                        <  1:574992/71‑1 (MQ=255)           tAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAAcg                                                       >  1:1044833/1‑71 (MQ=255)             gagCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAAcgc                                                      <  1:1994399/70‑1 (MQ=255)             gagCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAAcgc                                                      <  1:2079276/70‑1 (MQ=255)              agCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGa                                                    <  1:830026/71‑1 (MQ=255)                     cgcgCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCaa                                                                 <  1:2879039/51‑1 (MQ=255)                     cgcgCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCaa                                                                 <  1:1977561/51‑1 (MQ=255)                     cgcgCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCaa                                                                 <  1:1661230/51‑1 (MQ=255)                     cgcgCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAAc                                                        >  1:1805258/1‑60 (MQ=255)                     cgcgCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACg                                             >  1:1983614/1‑71 (MQ=255)                     cgcgCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACg                                             >  1:1455620/1‑71 (MQ=255)                       cgcTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGc                                           >  1:1671377/1‑71 (MQ=255)                        gcTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCg                                          <  1:1183968/71‑1 (MQ=255)                         cTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGt                                                <  1:603369/64‑1 (MQ=255)                         cTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCg                                          >  1:3581415/1‑70 (MQ=255)                         cTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGt                                         >  1:2995473/1‑71 (MQ=255)                           tttCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTg                                       >  1:3049564/1‑71 (MQ=255)                            ttCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATgg                                                              <  1:2753258/47‑1 (MQ=255)                            ttCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGtt                                        <  1:1378169/69‑1 (MQ=255)                            ttCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTc                                                                   <  1:2830000/42‑1 (MQ=255)                             tCCGCACCATCGACCGAAACATCGCGATTCTCATAATGCTc                                                                   <  1:3423472/41‑1 (MQ=255)                               cGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGAt                                                           >  1:524213/1‑47 (MQ=255)                               cGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGAt                                                           >  1:3371739/1‑47 (MQ=255)                               cGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGAt                                                           >  1:1105245/1‑47 (MQ=255)                               cGCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGt                                   <  1:2501436/71‑1 (MQ=255)                                gCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGt                                   >  1:1126574/1‑70 (MQ=255)                                gCACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTg                                  <  1:2186769/71‑1 (MQ=255)                                 cACCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGa                                 >  1:1067408/1‑71 (MQ=255)                                  aCCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGa                                 >  1:2837514/1‑70 (MQ=255)                                  aCCATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGa                                 >  1:1437648/1‑70 (MQ=255)                                   ccATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATg                               <  1:3225719/71‑1 (MQ=255)                                   ccATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATg                               <  1:3726709/71‑1 (MQ=255)                                   ccATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATg                               <  1:1673974/71‑1 (MQ=255)                                   ccATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGACAACGCGACAGTACCGCGTTGCCGTGATg                               <  1:3320668/71‑1 (MQ=255)                                   ccATCGACCGAAACATCGCGATTCTCATAATGCTc                                                                   <  1:575998/35‑1 (MQ=255)                                   ccATCGACCGAAACATCGCGATTCTCATAATGCTc                                                                   <  1:882771/35‑1 (MQ=255)                                    cATCGACCGAAACATCGCGATTCTCATAATGCTCATTGGGATAACGCGACAGTACGGc                                           <  1:1595817/58‑1 (MQ=255)                                    cATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGc                                           <  1:2250389/58‑1 (MQ=255)                                    cATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGc                                           <  1:8554/58‑1 (MQ=255)                                    cATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGc                                           <  1:3196697/58‑1 (MQ=255)                                    cATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGc                                           <  1:1697870/58‑1 (MQ=255)                                    cATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGc                              >  1:1435071/1‑71 (MQ=255)                                    cATCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGc                              >  1:194685/1‑71 (MQ=255)                                     aTCGACCGAAACATCGCGATTCTCATAATGCTCaa                                                                 >  1:348051/1‑35 (MQ=255)                                     aTCGACCGAAACATCGCGATTCTCATAATGCTCaa                                                                 >  1:153228/1‑35 (MQ=255)                                      tCGACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGccc                            >  1:98812/1‑71 (MQ=255)                                        gACCGAAACATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCt                           <  1:3561211/70‑1 (MQ=255)                                           cGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTcc                        <  1:3531882/70‑1 (MQ=255)                                            gAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCgg                      <  1:1483522/71‑1 (MQ=255)                                            gAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCgg                      <  1:2570513/71‑1 (MQ=255)                                             aaaCATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGt                                                <  1:27566/44‑1 (MQ=255)                                              aaCATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGcc                             >  1:3245699/1‑62 (MQ=255)                                              aaCATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCggg                     >  1:1325406/1‑70 (MQ=255)                                              aaCATCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCggg                     >  1:377415/1‑70 (MQ=255)                                                 aTCGCGATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCt                           <  1:2706447/61‑1 (MQ=255)                                                   cgcgATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGtat                  <  1:1257038/68‑1 (MQ=255)                                                   cgcgATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGtat                  <  1:2092740/68‑1 (MQ=255)                                                    gcgATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGc                              >  1:2558634/1‑55 (MQ=255)                                                     cgATTCTCATAATGTTCAATGGGATAACGCGACAGt                                                <  1:2823636/36‑1 (MQ=255)                                                      gATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCa            >  1:548030/1‑71 (MQ=255)                                                      gATTCTCATAATGCTCAATGGGATAACGCGACAGTACgg                                            >  1:735926/1‑39 (MQ=255)                                                      gATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGa                                 >  1:2750009/1‑50 (MQ=255)                                                      gATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGa                                 >  1:720968/1‑50 (MQ=255)                                                      gATTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCa            <  1:1802148/71‑1 (MQ=255)                                                       aTTCTCATAATGTTCAATGGGATAACGCGACAGTAc                                              >  1:3567938/1‑36 (MQ=255)                                                       aTTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCAt           <  1:2245004/71‑1 (MQ=255)                                                       aTTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCAt           <  1:3586062/71‑1 (MQ=255)                                                       aTTCTCATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGta                   <  1:1362303/63‑1 (MQ=255)                                                        ttCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGcc                                     <  1:134125/44‑1 (MQ=255)                                                           tcATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTg         <  1:2775130/69‑1 (MQ=255)                                                           tcATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTg         <  1:3435538/69‑1 (MQ=255)                                                            cATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCgg                      <  1:2373632/55‑1 (MQ=255)                                                            cATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCgg                      <  1:3525467/55‑1 (MQ=255)                                                            cATAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTGCGa      >  1:421813/1‑71 (MQ=255)                                                             aTAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTGCAAc     >  1:653878/1‑71 (MQ=255)                                                             aTAATGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGt                                   <  1:1577963/41‑1 (MQ=255)                                                               aaTGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCggg                     <  1:1404832/53‑1 (MQ=255)                                                               aaTGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCggg                     <  1:248963/53‑1 (MQ=255)                                                               aaTGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTGCGACCg   >  1:1625688/1‑71 (MQ=255)                                                               aaTGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTGCGACCg   >  1:525943/1‑71 (MQ=255)                                                               aaTGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTGCGACCg   >  1:2707439/1‑71 (MQ=255)                                                               aaTGCTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTg         >  1:33493/1‑65 (MQ=255)                                                                aTGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTGCGACCg   <  1:328167/70‑1 (MQ=255)                                                                aTGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTGCGACCGt  <  1:2536513/71‑1 (MQ=255)                                                                  gTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGa                                 <  1:1642613/38‑1 (MQ=255)                                                                  gTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGa                                 <  1:258265/38‑1 (MQ=255)                                                                   |                                                                    AATGCGGCAGTAGAGCACGCCGCGCTTTTCCGCACCATCGACCGAAACATCGCGATTCTCATAATGTTCAATGGGATAACGCGACAGTACGGCGTTGCCGTGATGCCCTTCCGGGTATACGGCATTGCGACCGT  >  minE/530572‑530705

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