Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,581,079	C→T	27.4%	G47S (GGT→AGT)	ypeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,581,079	0	C	T	27.4%	78.8 / 31.0	51	G47S (GGT→AGT)	ypeB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/20);  new base T (8/6);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

CTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGATTAAG  >  minE/1581011‑1581149                                                                     |                                                                       cTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGa                                                                      <  1:1793656/71‑1 (MQ=255)     aGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCtt                                                                  >  1:1514260/1‑71 (MQ=255)      ggggCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTc                                                                 >  1:934649/1‑71 (MQ=255)       gggCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTcc                                                                <  1:3597845/71‑1 (MQ=255)           tGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGc                                                            >  1:2946125/1‑71 (MQ=255)               ccAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCggg                                                        <  1:3255290/71‑1 (MQ=255)               ccAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCgg                                                         >  1:3690802/1‑70 (MQ=255)               ccAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAACTCACTGAGCTTCCCTGCCggg                                                        <  1:3244764/71‑1 (MQ=255)                 aGTTGCGTCAGCCCCTCTGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGc                                                            >  1:1436316/1‑65 (MQ=255)                 aGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGc                                                                    <  1:1405855/57‑1 (MQ=255)                 aGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGc                                                            >  1:2697202/1‑65 (MQ=255)                   ttGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGc                                                            <  1:666450/63‑1 (MQ=255)                   ttGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAAc                                                     >  1:311127/1‑70 (MQ=255)                   ttGCGTCAGCACCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGc                                                            <  1:2090036/63‑1 (MQ=255)                     gCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCa                                                  >  1:1575433/1‑71 (MQ=255)                     gCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTc                                                   >  1:1330712/1‑70 (MQ=255)                        tCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCAcc                                                                        >  1:2226669/1‑46 (MQ=255)                         cAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCtt                                              >  1:2442642/1‑71 (MQ=255)                          aGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTc                                                   <  1:17448/65‑1 (MQ=255)                          aGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCttt                                             >  1:2849829/1‑71 (MQ=255)                          aGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCttt                                             <  1:2367233/71‑1 (MQ=255)                           gCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGa                                                                      >  1:84574/1‑45 (MQ=255)                           gCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTg                                            <  1:1379140/71‑1 (MQ=255)                            ccccTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGc                                           <  1:3338578/71‑1 (MQ=255)                             cccTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTg                                                             >  1:83958/1‑52 (MQ=255)                             cccTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCg                                                          >  1:2738055/1‑55 (MQ=255)                              ccTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCgg                                         <  1:885503/71‑1 (MQ=255)                               cTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCa                                                  >  1:3652493/1‑61 (MQ=255)                                 cGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGc                                       >  1:616796/1‑70 (MQ=255)                                  ggTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                      <  1:2957742/70‑1 (MQ=255)                                     tttAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCAt                                                 <  1:553031/56‑1 (MQ=255)                                       tAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACca                                 <  1:679652/70‑1 (MQ=255)                                            agcagcGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcag                              >  1:3733809/1‑68 (MQ=255)                                            agcagcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcagcag                           <  1:2145777/71‑1 (MQ=255)                                             gcagcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGa                          <  1:796976/71‑1 (MQ=255)                                              cagcGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGAt                         <  1:2628654/71‑1 (MQ=255)                                               agcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCg                                          >  1:2058711/1‑53 (MQ=255)                                               agcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGAt                         <  1:263299/70‑1 (MQ=255)                                                   tGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATActcc                    >  1:1010259/1‑71 (MQ=255)                                                     aTTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATActcctc                  <  1:125325/71‑1 (MQ=255)                                                      tttGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACtcctcc                 >  1:464592/1‑71 (MQ=255)                                                        tGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTg                                            >  1:694000/1‑42 (MQ=255)                                                         gCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCgg               <  1:123658/70‑1 (MQ=255)                                                         gCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCgg               <  1:131227/70‑1 (MQ=255)                                                           aTTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCa             <  1:3241105/70‑1 (MQ=255)                                                            ttAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                      >  1:3428682/1‑44 (MQ=255)                                                             tAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGt          >  1:2264073/1‑71 (MQ=255)                                                              aGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTc         >  1:3330924/1‑71 (MQ=255)                                                               gCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCg        <  1:1526817/71‑1 (MQ=255)                                                                 tCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATAc                       <  1:3196625/54‑1 (MQ=255)                                                                  cACTGAGCTTCCCTGCCGGGAACTCATCTTTGCggg                                        >  1:3496314/1‑36 (MQ=255)                                                                   aCCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGATTa    <  1:101183/71‑1 (MQ=255)                                                                     cGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGATTAAg  <  1:3545447/71‑1 (MQ=255)                                                                     |                                                                       CTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGATTAAG  >  minE/1581011‑1581149

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