Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	282,462	C→T	17.9%	A59T (GCA→ACA)	yajG ←	predicted lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	282,462	0	C	T	17.9%	190.7 / 32.8	78	A59T (GCA→ACA)	yajG	predicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (36/28);  new base T (8/6);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

CTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCT  >  minE/282393‑282530                                                                      |                                                                     cttGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGc                                                                      >  1:1069619/1‑70 (MQ=255)      agcagGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAgcgc                                                                >  1:2719960/1‑71 (MQ=255)       gcagGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGcc                                                               >  1:2186526/1‑71 (MQ=255)         agGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGcc                                                               >  1:1331275/1‑69 (MQ=255)          ggAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGa                                                            <  1:2939590/71‑1 (MQ=255)          ggAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGa                                                            >  1:1602975/1‑71 (MQ=255)            aaaCGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGa                                                            <  1:2002307/69‑1 (MQ=255)              aCGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCgg                                                        <  1:757635/71‑1 (MQ=255)                gCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCAGCGCCTGAt                                                           >  1:871934/1‑66 (MQ=255)                gCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCAGCGCCTGAt                                                           >  1:2270406/1‑66 (MQ=255)                gCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTa                                                      >  1:106940/1‑71 (MQ=255)                 cAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTAc                                                     >  1:3801407/1‑71 (MQ=255)                   gATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCAGCGCCTGATc                                                          <  1:2269860/64‑1 (MQ=255)                   gATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCAGCGCCTGATc                                                          <  1:287223/64‑1 (MQ=255)                    aTCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCAGCGCCTGATCgg                                                        <  1:421485/65‑1 (MQ=255)                    aTCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCt                                                  <  1:3963094/71‑1 (MQ=255)                    aTCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCt                                                  >  1:604603/1‑71 (MQ=255)                     tCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCAGCGCCTGATCGGTACGCTg                                                 <  1:2097120/71‑1 (MQ=255)                      cgcgGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCAGCGCCTGATCGGTACGCTg                                                 <  1:407177/70‑1 (MQ=255)                      cgcgGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCGGCGCCTGATCGGTACGCTGa                                                >  1:3851905/1‑71 (MQ=255)                      cgcgGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTAc                                                     >  1:961133/1‑66 (MQ=255)                      cgcgGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTAc                                                     >  1:149266/1‑66 (MQ=255)                      cgcgGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGa                                                >  1:3907630/1‑71 (MQ=255)                       gcgGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGa                                                <  1:4116317/70‑1 (MQ=255)                       gcgGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGAt                                               <  1:2979619/71‑1 (MQ=255)                         gggAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCt                                                              >  1:1913321/1‑54 (MQ=255)                                ggTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCAt                                      <  1:579801/71‑1 (MQ=255)                                ggTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCAt                                      <  1:2316596/71‑1 (MQ=255)                                 gTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCAtt                                     <  1:1993380/71‑1 (MQ=255)                                 gTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCAtt                                     <  1:1712359/71‑1 (MQ=255)                                 gTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCAtt                                     <  1:410970/71‑1 (MQ=255)                                 gTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCAtt                                     <  1:2559659/71‑1 (MQ=255)                                   cAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCa                                                                    >  1:1916990/1‑38 (MQ=255)                                   cAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCa                                                                    >  1:3602383/1‑38 (MQ=255)                                   cAGGGTAACGATTTGATTATCGCGGGTGACTTTTGTCa                                                                    >  1:457659/1‑38 (MQ=255)                                   cAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTaa                                   >  1:3551658/1‑71 (MQ=255)                                   cAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTaa                                   >  1:1066993/1‑71 (MQ=255)                                    aGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTaa                                   <  1:4074390/70‑1 (MQ=255)                                    aGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAAt                                  >  1:2975466/1‑71 (MQ=255)                                       gTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCt                               <  1:1032969/71‑1 (MQ=255)                                       gTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCt                               <  1:2407208/71‑1 (MQ=255)                                        tAACGATTTGATTATCGCGGGTGACTTTTGTCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATgat                               >  1:635583/1‑68 (MQ=255)                                        tAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCAc                                         <  1:276542/60‑1 (MQ=255)                                        tAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCAc                                         <  1:1386655/60‑1 (MQ=255)                                        tAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCtt                              <  1:3912785/71‑1 (MQ=255)                                        tAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCtt                              <  1:362214/71‑1 (MQ=255)                                         aaCGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTaa                                   >  1:2559202/1‑65 (MQ=255)                                          aCGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTaa                                   >  1:418862/1‑64 (MQ=255)                                          aCGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTa                             <  1:1403204/70‑1 (MQ=255)                                          aCGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTAc                            >  1:2298961/1‑71 (MQ=255)                                          aCGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTAc                            >  1:2409520/1‑71 (MQ=255)                                           cGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCgg                                                        <  1:2036183/42‑1 (MQ=255)                                            gATTTGATTATCGCTGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGAt                                               >  1:2295176/1‑50 (MQ=255)                                            gATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGAt                                               >  1:3954410/1‑50 (MQ=255)                                            gATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGAt                                               >  1:1196651/1‑50 (MQ=255)                                            gATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACgg                          <  1:2351883/71‑1 (MQ=255)                                               ttGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAAt                                  >  1:1168418/1‑60 (MQ=255)                                                   ttATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGt                         >  1:1152289/1‑65 (MQ=255)                                                   ttATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGt                         >  1:2858532/1‑65 (MQ=255)                                                   ttATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGt                         >  1:3432829/1‑65 (MQ=255)                                                     aTCGCGGGTGACTTTTGTCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGcc                   >  1:511560/1‑69 (MQ=255)                                                     aTCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGccc                  <  1:2499396/70‑1 (MQ=255)                                                     aTCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGccc                  <  1:2504248/70‑1 (MQ=255)                                                      tCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCa                 <  1:3320635/70‑1 (MQ=255)                                                            gTGACTTTTGTCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCt          <  1:3165808/71‑1 (MQ=255)                                                              gACTTTTGTCAGCGCCTGATAGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTg        >  1:671644/1‑71 (MQ=255)                                                                cTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGa      <  1:1025786/71‑1 (MQ=255)                                                                cTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGa      >  1:818986/1‑71 (MQ=255)                                                                cTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGa      >  1:1522316/1‑71 (MQ=255)                                                                 ttttGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGAt     >  1:952939/1‑71 (MQ=255)                                                                 ttttGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGAt     >  1:3680794/1‑71 (MQ=255)                                                                   ttGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGAt     <  1:1644032/69‑1 (MQ=255)                                                                   ttGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATcc   >  1:578264/1‑71 (MQ=255)                                                                   ttGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATcc   >  1:1522555/1‑71 (MQ=255)                                                                   ttGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATcc   >  1:2681555/1‑71 (MQ=255)                                                                   ttGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATcc   >  1:1753554/1‑71 (MQ=255)                                                                   ttGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATcc   >  1:3534799/1‑71 (MQ=255)                                                                    tGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATcc   <  1:866117/70‑1 (MQ=255)                                                                     gCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCct  >  1:3827448/1‑70 (MQ=255)                                                                      |                                                                     CTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCT  >  minE/282393‑282530

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑