Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,683	C→A	7.4%	P134P (CCG→CCT) ‡	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,683	0	C	A	7.4%	312.2 / 5.6	94	P134P (CCG→CCT) ‡	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (48/39);  new base A (5/2);  total (53/41)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

AGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGA  >  minE/454613‑454753                                                                       |                                                                       aGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGAcgc                                                                        <  1:2857253/71‑1 (MQ=255)    cGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGa                                                                     <  1:2740831/71‑1 (MQ=255)     gATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGAt                                                                    >  1:828148/1‑71 (MQ=255)       tAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGAt                                                                    <  1:1862105/69‑1 (MQ=255)       tAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCg                                                                  >  1:3902234/1‑71 (MQ=255)       tAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCg                                                                  >  1:1665292/1‑71 (MQ=255)         ggTGCGATAATAGACCCCACTTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGt                                                                <  1:549847/71‑1 (MQ=255)         ggTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCgg                                                                 <  1:1987309/70‑1 (MQ=255)         ggTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCgg                                                                 <  1:900590/70‑1 (MQ=255)          gTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCgg                                                                      >  1:1140057/1‑64 (MQ=255)          gTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGt                                                                <  1:3295401/70‑1 (MQ=255)          gTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGt                                                                <  1:524166/70‑1 (MQ=255)             cGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGc                                                            >  1:1172243/1‑71 (MQ=255)             cGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGc                                                            >  1:1569829/1‑71 (MQ=255)               ataataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGAcgcg                                                                       <  1:1952022/58‑1 (MQ=255)                taataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                         >  1:2161032/1‑71 (MQ=255)                 aataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCgg                                                                 >  1:1899291/1‑62 (MQ=255)                     gACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGtt                                                               >  1:2036705/1‑60 (MQ=255)                     gACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                         <  1:184754/66‑1 (MQ=255)                        cccACGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATa                                                       <  1:859548/65‑1 (MQ=255)                        cccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTc                                                 >  1:3847756/1‑71 (MQ=255)                          cACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGc                                               >  1:608391/1‑71 (MQ=255)                           aCGTAACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTg                                                              >  1:831538/1‑55 (MQ=255)                            cGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCg                                             <  1:1312282/71‑1 (MQ=255)                             gTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGa                                            <  1:2009645/71‑1 (MQ=255)                              ttACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                           <  1:276893/71‑1 (MQ=255)                                aCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                           >  1:1799032/1‑69 (MQ=255)                                aCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                           >  1:373191/1‑69 (MQ=255)                                 cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                           >  1:2578201/1‑68 (MQ=255)                                 cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa                                         <  1:2559078/70‑1 (MQ=255)                                  ccgccgCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                       >  1:2965736/1‑71 (MQ=255)                                  ccgccgCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                       >  1:2391635/1‑71 (MQ=255)                                   cgccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                       >  1:3608211/1‑70 (MQ=255)                                    gccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGcc                                              <  1:2241698/62‑1 (MQ=255)                                    gccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                       >  1:3541892/1‑69 (MQ=255)                                    cccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCTGTTGGCTCATACAGCTCGcc                                              <  1:3619297/61‑1 (MQ=255)                                     cagcATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAg                                    >  1:2326625/3‑71 (MQ=255)                                      cgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGc                                   >  1:2477066/1‑71 (MQ=255)                                       gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                  <  1:829819/71‑1 (MQ=255)                                       gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                  >  1:405055/1‑71 (MQ=255)                                       gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                  >  1:1069030/1‑71 (MQ=255)                                       gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                  <  1:3281296/71‑1 (MQ=255)                                       gcATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCt                                  >  1:3766968/1‑71 (MQ=255)                                        cATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCtt                                 <  1:3802070/71‑1 (MQ=255)                                         aTACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCtga                                <  1:1096568/71‑3 (MQ=255)                                           aCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGc                                                   <  1:2158513/50‑1 (MQ=255)                                           aCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACg                              >  1:3177783/1‑71 (MQ=255)                                           aCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTAACg                              >  1:165886/1‑71 (MQ=255)                                            cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                             >  1:3877817/1‑71 (MQ=255)                                            cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                             <  1:1233663/71‑1 (MQ=255)                                            cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                             <  1:501747/71‑1 (MQ=255)                                             aTCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATAc                                                      <  1:1081706/45‑1 (MQ=255)                                               cAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                             >  1:1202482/1‑68 (MQ=255)                                               cAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                             >  1:2540867/1‑68 (MQ=255)                                                aGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                             >  1:2064997/1‑67 (MQ=255)                                                 gCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                         <  1:3110265/38‑1 (MQ=255)                                                 gCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                  <  1:3103711/61‑1 (MQ=255)                                                  cGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCg                                             >  1:2622821/1‑49 (MQ=255)                                                  cGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtt                       >  1:402654/1‑71 (MQ=255)                                                  cGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtt                       >  1:3321801/1‑71 (MQ=255)                                                   gagaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                           <  1:3643367/50‑1 (MQ=255)                                                   gagaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAtt                       >  1:1724687/1‑70 (MQ=255)                                                    agaAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                         >  1:78781/1‑35 (MQ=255)                                                    agaAACCGCCATCACGACGCGGATCGGTTGGCTCATAc                                                      >  1:2698670/1‑38 (MQ=255)                                                    agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGc                                                   <  1:3723749/41‑1 (MQ=255)                                                    agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAt                        >  1:1792505/1‑68 (MQ=255)                                                    aaaaaCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtttt                     <  1:407424/69‑2 (MQ=255)                                                     gaAACCGCCATCACGACGCGGATCGGTTGGCTCATTCAGCTCGCCGAGCAGACAGCTTACGTTAGAt                        >  1:2123033/1‑67 (MQ=255)                                                     gaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaa                    <  1:559351/71‑1 (MQ=255)                                                        aCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGtt                            >  1:2765125/1‑60 (MQ=255)                                                        aCCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCg                                                >  1:1895018/1‑40 (MQ=255)                                                         ccGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAg                <  1:1102261/71‑1 (MQ=255)                                                         ccGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAg                <  1:683214/71‑1 (MQ=255)                                                          cGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                  <  1:2614407/52‑1 (MQ=255)                                                          cGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAgg               >  1:2730951/1‑71 (MQ=255)                                                          cGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAg                <  1:3693585/70‑1 (MQ=255)                                                          cGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCAGAGCAGACAGCTTACGTTAGATTTAAACAgg               >  1:2724593/1‑71 (MQ=255)                                                           gCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAt                        >  1:2097698/1‑61 (MQ=255)                                                           gCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAt                        >  1:3615566/1‑61 (MQ=255)                                                            ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                  <  1:484344/50‑1 (MQ=255)                                                            ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                >  1:783409/1‑52 (MQ=255)                                                            ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                >  1:968076/1‑52 (MQ=255)                                                            ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                >  1:50847/1‑52 (MQ=255)                                                             cATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa                                         <  1:3866321/42‑1 (MQ=255)                                                             cATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCa            <  1:1801350/71‑1 (MQ=255)                                                              aTCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGc                                          >  1:3345456/1‑40 (MQ=255)                                                              aTCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCAt           <  1:813998/71‑1 (MQ=255)                                                               tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCa            >  1:246003/1‑69 (MQ=255)                                                               tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATa          <  1:2751043/71‑1 (MQ=255)                                                               tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATa          <  1:3780608/71‑1 (MQ=255)                                                                cACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACa                 >  1:968035/1‑63 (MQ=255)                                                                 acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa                   >  1:1682499/1‑60 (MQ=255)                                                                 acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGc        >  1:3451288/1‑71 (MQ=255)                                                                   gacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCg      <  1:5780/71‑1 (MQ=255)                                                                    acgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAg                                    >  1:2425160/1‑40 (MQ=255)                                                                    acgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACATCTTACGTTAGATTTa                     <  1:2108123/55‑1 (MQ=255)                                                                     cgcgGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAg                                    >  1:3792344/1‑39 (MQ=255)                                                                     cgcgGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAg                                    >  1:870356/1‑39 (MQ=255)                                                                     cgcgGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCg      >  1:1228929/1‑69 (MQ=255)                                                                       cgGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGa  <  1:2779507/71‑1 (MQ=255)                                                                       |                                                                       AGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGA  >  minE/454613‑454753

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