Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	281,469	T→C	33.3%	A182A (GCA→GCG)	ampG ←	muropeptide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	281,469	0	T	C	33.3%	52.5 / 38.3	72	A182A (GCA→GCG)	ampG	muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base T (27/21);  new base C (22/2);  total (49/23)
Fisher's exact test for biased strand distribution p-value = 2.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TTGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTAG  >  minE/281401‑281531                                                                     |                                                               ttGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGt                                                                >  1:2012831/1‑69 (MQ=255) ttGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGt                                                                >  1:2066407/1‑69 (MQ=255) ttGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGt                                                                >  1:2178764/1‑69 (MQ=255) ttGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGt                                                                >  1:1535870/1‑69 (MQ=255) ttGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGt                                                                >  1:544178/1‑69 (MQ=255) ttGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGt                                                                >  1:2435458/1‑69 (MQ=255) ttGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGt                                                                >  1:2439342/1‑69 (MQ=255) ttGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGt                                                                >  1:877067/1‑69 (MQ=255) ttGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAAAAGt                                                                >  1:3125965/1‑69 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:2643198/1‑68 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:689290/1‑68 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:998911/1‑68 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:656745/1‑68 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:2756104/1‑68 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:2931662/1‑68 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:3160279/1‑68 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:3194630/1‑68 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:3391996/1‑68 (MQ=255)        caGGAATGGTGTCGGTGGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCgccgcc                                                          >  1:1425029/1‑68 (MQ=255)          ggAATGGTGTCGGTTGGTTCTGGTGCGAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                       >  1:2327439/1‑69 (MQ=255)          ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACGGTGCCGCCAtt                                                       >  1:3076646/1‑69 (MQ=255)          ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                       >  1:3311158/1‑69 (MQ=255)          ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                       >  1:2702835/1‑69 (MQ=255)          ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                       >  1:2881519/1‑69 (MQ=255)          ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                       >  1:2600785/1‑69 (MQ=255)          ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                       >  1:2966599/1‑69 (MQ=255)          ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                       >  1:3193327/1‑69 (MQ=255)          ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                       >  1:478768/1‑69 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:2551033/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:2769641/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:1356333/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:1224604/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:1958790/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:2024702/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:881674/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:3041941/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:2052063/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:3106592/1‑67 (MQ=255)             aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                      >  1:2187735/1‑67 (MQ=255)                gtgtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAg                                                 <  1:1792466/69‑1 (MQ=255)                gtgtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAg                                                 <  1:3254765/69‑1 (MQ=255)                gtgtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAg                                                 <  1:2286368/69‑1 (MQ=255)                  gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                               <  1:3050349/69‑1 (MQ=255)                  gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                               <  1:115679/69‑1 (MQ=255)                  gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                               <  1:412492/69‑1 (MQ=255)                  gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                               <  1:186186/69‑1 (MQ=255)                  gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                               <  1:1875300/69‑1 (MQ=255)                  gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                               <  1:1687028/69‑1 (MQ=255)                  gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                               <  1:882946/69‑1 (MQ=255)                  gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                               <  1:1506766/69‑1 (MQ=255)                     ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                    <  1:181923/61‑1 (MQ=255)                     ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                    <  1:582603/61‑1 (MQ=255)                     ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                    <  1:3409504/61‑1 (MQ=255)                     ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                    <  1:1098489/61‑1 (MQ=255)                     ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                    <  1:1266137/61‑1 (MQ=255)                     ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACAt                                            >  1:2261835/1‑69 (MQ=255)                             cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                    >  1:1515012/1‑69 (MQ=255)                             cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                    >  1:3105012/1‑69 (MQ=255)                             cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                    >  1:3137295/1‑69 (MQ=255)                             cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                    >  1:249397/1‑69 (MQ=255)                             cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                    >  1:1084346/1‑69 (MQ=255)                             cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                    >  1:814367/1‑69 (MQ=255)                             cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                    >  1:1170801/1‑69 (MQ=255)                             cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                    >  1:335578/1‑69 (MQ=255)                             cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                    >  1:1415965/1‑69 (MQ=255)                                gTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGccag                                  >  1:1654413/1‑68 (MQ=255)                                gTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGccag                                  >  1:753285/1‑68 (MQ=255)                                gTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGccag                                  >  1:2824454/1‑68 (MQ=255)                                   caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                              >  1:1045652/1‑69 (MQ=255)                                   caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                              >  1:361394/1‑69 (MQ=255)                                   caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                              >  1:2491317/1‑69 (MQ=255)                                   caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                              >  1:2108664/1‑69 (MQ=255)                                           gTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTa                      <  1:3219127/69‑1 (MQ=255)                                           gTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTa                      <  1:1488934/69‑1 (MQ=255)                                           gTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTa                      <  1:1013808/69‑1 (MQ=255)                                                    ataCAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTgccagc             <  1:2824793/69‑1 (MQ=255)                                                    ataCAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTgccagc             >  1:2615649/1‑69 (MQ=255)                                                      aCAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTgccagc             >  1:2838540/1‑67 (MQ=255)                                                       cAGGGGATCAACAGTGTCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCcac          <  1:2651231/69‑1 (MQ=255)                                                               cAACAGCGCCGCCATCAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTaa  <  1:2616209/69‑2 (MQ=255)                                                                aaCAGCGCCGCCATCAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTaa  <  1:2592896/68‑2 (MQ=255)                                                                     |                                                               TTGGGCACAGGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTAG  >  minE/281401‑281531

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