Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,456,196	G→A	30.0%	A385V (GCG→GTG)	yfaQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,456,196	0	G	A	30.0%	56.1 / 28.9	40	A385V (GCG→GTG)	yfaQ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base A (5/7);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATC  >  minE/1456133‑1456263                                                                |                                                                    ccATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGt                                                              >  1:1372266/1‑71 (MQ=255) ccATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCAccccc                                                                <  1:288163/69‑1 (MQ=255)  cATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCAccccc                                                                <  1:1488925/68‑1 (MQ=255)     gATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTc                                                          >  1:1034710/1‑71 (MQ=255)      aTAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCg                                                         >  1:2696480/1‑71 (MQ=255)       tAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCg                                                         <  1:1700055/70‑1 (MQ=255)         aTGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCt                                                       <  1:3206561/70‑1 (MQ=255)          tGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTgg                                                     >  1:2698159/1‑71 (MQ=255)                 gATCGCCTGCATAGATGAGAACCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAg                                                >  1:589219/1‑69 (MQ=255)                   tCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCa                                              <  1:3176034/69‑1 (MQ=255)                   tCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAg                                             >  1:3291672/1‑70 (MQ=255)                    cGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                            >  1:544868/1‑70 (MQ=255)                     gCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGaaa                                          <  1:3707829/71‑1 (MQ=255)                     gCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAATCAGaaa                                          <  1:2313187/71‑1 (MQ=255)                     gCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAATCAGaaa                                          <  1:2317846/71‑1 (MQ=255)                     gCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGaaa                                          <  1:2849571/71‑1 (MQ=255)                          cATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCg                                                         >  1:3245271/1‑51 (MQ=255)                          cATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAacac                                       <  1:2302881/69‑1 (MQ=255)                              gATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCg                                 >  1:1753/1‑71 (MQ=255)                                tGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTg                               >  1:376724/1‑71 (MQ=255)                                     tCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGaaa                                          >  1:2264076/1‑55 (MQ=255)                                      ccTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTg                                                             >  1:233338/1‑35 (MQ=255)                                      ccTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGa                                            >  1:3541696/1‑52 (MQ=255)                                        tGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                       >  1:3384096/1‑71 (MQ=255)                                                ggCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGa               <  1:2760408/71‑1 (MQ=255)                                                ggCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGa               >  1:81875/1‑71 (MQ=255)                                                 gCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgct                          <  1:3394022/59‑1 (MQ=255)                                                  cGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAaca                                        >  1:2687241/1‑44 (MQ=255)                                                    gTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGt           <  1:1065438/71‑1 (MQ=255)                                                    gTCATCGTTCGCGCCCCCGTGGTCGCGGGCGAAGCAGAAACACGCTGCGTGGCgct                          <  1:247244/56‑1 (MQ=255)                                                     tCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTc          <  1:1115607/71‑1 (MQ=255)                                                      cATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGt           >  1:4064854/1‑69 (MQ=255)                                                       aTCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAg        <  1:3931030/71‑1 (MQ=255)                                                       aTCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAg        <  1:456724/71‑1 (MQ=255)                                                       aTCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAg        <  1:4136144/71‑1 (MQ=255)                                                       aTCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGa               <  1:1843789/64‑1 (MQ=255)                                                         cGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAgg       >  1:1057479/1‑70 (MQ=255)                                                          gTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAaca                                        >  1:1657885/1‑36 (MQ=255)                                                          gTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGc      >  1:3607065/1‑70 (MQ=255)                                                           ttCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAg        <  1:2727985/67‑1 (MQ=255)                                                             cgcgCCCCCGTGGTCGCTGGGGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATc  >  1:725821/1‑71 (MQ=255)                                                              gcgcCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATc  <  1:4122207/70‑1 (MQ=255)                                                              gcgcCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATc  <  1:1497945/70‑1 (MQ=255)                                                                |                                                                    CCATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATC  >  minE/1456133‑1456263

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