Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,715	G→A	7.2%	L124L (CTG→TTG)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,715	0	G	A	7.2%	296.6 / 8.2	97	L124L (CTG→TTG)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (54/36);  new base A (5/2);  total (59/38)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATAGAGGAAT  >  minE/454646‑454782                                                                      |                                                                    ccgccgCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                                                    >  1:2965736/1‑71 (MQ=255) ccgccgCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                                                    >  1:2391635/1‑71 (MQ=255)  cgccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                                                    >  1:3608211/1‑70 (MQ=255)   gccgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                                                    >  1:3541892/1‑69 (MQ=255)    cagcATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAg                                                                 >  1:2326625/3‑71 (MQ=255)     cgcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGc                                                                >  1:2477066/1‑71 (MQ=255)      gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                                               >  1:1069030/1‑71 (MQ=255)      gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                                               <  1:3281296/71‑1 (MQ=255)      gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                                               >  1:405055/1‑71 (MQ=255)      gcATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                                               <  1:829819/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)          aCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACg                                                           >  1:3177783/1‑71 (MQ=255)          aCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTAACg                                                           >  1:165886/1‑71 (MQ=255)           cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                                          >  1:3877817/1‑71 (MQ=255)           cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                                          <  1:501747/71‑1 (MQ=255)           cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                                          <  1:1233663/71‑1 (MQ=255)              cAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                                          >  1:2540867/1‑68 (MQ=255)              cAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                                          >  1:1202482/1‑68 (MQ=255)               aGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                                                          >  1:2064997/1‑67 (MQ=255)                gCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                                               <  1:3103711/61‑1 (MQ=255)                 cGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtt                                                    >  1:3321801/1‑71 (MQ=255)                 cGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtt                                                    >  1:402654/1‑71 (MQ=255)                  gagaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAtt                                                    >  1:1724687/1‑70 (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)                        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:3615566/1‑61 (MQ=255)                          gCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAt                                                     >  1:2097698/1‑61 (MQ=255)                           ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                                               <  1:484344/50‑1 (MQ=255)                           ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                                             >  1:50847/1‑52 (MQ=255)                           ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                                             >  1:783409/1‑52 (MQ=255)                           ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                                             >  1:968076/1‑52 (MQ=255)                            cATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCa                                         <  1:1801350/71‑1 (MQ=255)                             aTCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCAt                                        <  1:813998/71‑1 (MQ=255)                              tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCa                                         >  1:246003/1‑69 (MQ=255)                              tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATa                                       <  1:3780608/71‑1 (MQ=255)                              tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATa                                       <  1:2751043/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)                                         aTCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGa                            <  1:2758569/71‑1 (MQ=255)                                          tCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGa                            >  1:709944/1‑70 (MQ=255)                                          tCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGAt                           <  1:476305/71‑1 (MQ=255)                                             gTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGtt                                                         <  1:2172007/38‑1 (MQ=255)                                             gTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAttt                                                   <  1:2049359/44‑1 (MQ=255)                                             gTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAttt                                                   <  1:645129/44‑1 (MQ=255)                                                 gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa                                                >  1:2483421/1‑43 (MQ=255)                                                 gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGcc                    <  1:816773/71‑1 (MQ=255)                                                 gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGcc                    <  1:3879563/71‑1 (MQ=255)                                                 gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGcc                    <  1:386479/71‑1 (MQ=255)                                                 gCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGc                     <  1:1096232/70‑1 (MQ=255)                                                 ccTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGc                     <  1:3247900/69‑1 (MQ=255)                                                  cTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGc                     >  1:281619/1‑69 (MQ=255)                                                    cATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAgg                 >  1:2851010/1‑71 (MQ=255)                                                    cATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAgg                 >  1:3742428/1‑71 (MQ=255)                                                    cATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAgg                 >  1:311638/1‑71 (MQ=255)                                                    cATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAgg                 >  1:1456961/1‑71 (MQ=255)                                                    cATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAgg                 >  1:316507/1‑71 (MQ=255)                                                     aTACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTATACAGGCCATAGCCGCAGACGATAATCGCCAgg                 <  1:3081533/70‑1 (MQ=255)                                                       aCAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAg                  <  1:2941641/67‑1 (MQ=255)                                                        cAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa                                                >  1:599382/1‑36 (MQ=255)                                                        cAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa                                                >  1:1924399/1‑36 (MQ=255)                                                        cAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGa               >  1:2952482/1‑69 (MQ=255)                                                        cAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGa               >  1:3560681/1‑69 (MQ=255)                                                        cAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGa               >  1:1946835/1‑69 (MQ=255)                                                        cAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGa               >  1:355510/1‑69 (MQ=255)                                                        cAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGAc              >  1:2563548/1‑70 (MQ=255)                                                        cAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACa             <  1:3939566/71‑1 (MQ=255)                                                        cAGCTCGCCGAGCAGACAGCTAACGTTAGATTTAAACa                                              >  1:3810630/1‑38 (MQ=255)                                                          gCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGa           <  1:77958/71‑1 (MQ=255)                                                          gCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGa           <  1:1802884/71‑1 (MQ=255)                                                           cTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAg                                             >  1:1652159/1‑36 (MQ=255)                                                           cTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGa                               >  1:2389554/1‑50 (MQ=255)                                                              gCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATa         >  1:1249316/1‑69 (MQ=255)                                                              gCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGAAAGATa         >  1:902123/1‑69 (MQ=255)                                                               ccGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATagag      >  1:741722/1‑71 (MQ=255)                                                               ccGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGAAGACGATAATCGCCAGGGACAGATagag      >  1:244683/1‑71 (MQ=255)                                                                cGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATAGAgg     <  1:929960/71‑1 (MQ=255)                                                                  aGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGAc              >  1:3633405/1‑60 (MQ=255)                                                                   gCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATAGAGGAAt  <  1:3100512/71‑1 (MQ=255)                                                                      |                                                                    CCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATAGAGGAAT  >  minE/454646‑454782

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