Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	412,090	A→G	36.5%	V89A (GTG→GCG)	leuS ←	leucyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	412,090	0	A	G	36.5%	54.7 / 60.7	65	V89A (GTG→GCG)	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (19/22);  new base G (13/10);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

AGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGC  >  minE/412021‑412158                                                                      |                                                                     aGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCAc                                                                     >  1:1947093/1‑71 (MQ=255)  gCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCcgcg                                                                    <  1:3619298/71‑1 (MQ=255)   cTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCa                                                                      >  1:2443163/1‑68 (MQ=255)   cTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGc                                                                   >  1:2409186/1‑71 (MQ=255)      gTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCtt                                                                <  1:1701191/71‑1 (MQ=255)      gTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGcc                                                                  <  1:3262018/69‑1 (MQ=255)      gTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCtt                                                                >  1:3058816/1‑71 (MQ=255)       tttttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCt                                                                 <  1:3425971/69‑1 (MQ=255)        ttttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTc                                                               >  1:412271/1‑70 (MQ=255)        ttttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTc                                                               >  1:582363/1‑70 (MQ=255)          ttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCg                                                             <  1:147967/70‑1 (MQ=255)          ttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCg                                                             <  1:877183/70‑1 (MQ=255)          ttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCg                                                             >  1:1136185/1‑70 (MQ=255)          ttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCATTCCg                                                             >  1:196851/1‑70 (MQ=255)           tCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCg                                                             >  1:1995615/1‑69 (MQ=255)           tCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCg                                                             >  1:2990356/1‑69 (MQ=255)           tCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCg                                                             >  1:3636103/1‑69 (MQ=255)           tCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTT‑TTTTTCACCGCCACGCCTTCCgca                                                           <  1:3245233/70‑1 (MQ=255)            cATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCgca                                                           <  1:977271/70‑1 (MQ=255)             atatACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcagg                                                         >  1:2658167/1‑71 (MQ=255)             atatACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcagg                                                         <  1:14699/71‑1 (MQ=255)               ataCGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCCCCTTCCgca                                                           >  1:3312806/1‑67 (MQ=255)                taCGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcaggca                                                       <  1:118241/70‑1 (MQ=255)                taCGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcaggca                                                       <  1:122943/70‑1 (MQ=255)                taCGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcaggc                                                        <  1:1693500/69‑1 (MQ=255)                taCGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcaggaa                                                       <  1:1618542/70‑3 (MQ=255)                    cgATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGAcc                                                   >  1:1550515/1‑70 (MQ=255)                    cgATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTACGCAGGCAGACCa                                                  <  1:1093389/71‑1 (MQ=255)                          tgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAAcgcg                                            >  1:2922159/1‑71 (MQ=255)                                aCGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAg                                      >  1:3001176/1‑71 (MQ=255)                                 cGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGc                                     >  1:288896/1‑71 (MQ=255)                                   tCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCg                                   <  1:1775571/71‑1 (MQ=255)                                   tCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCg                                   <  1:3508307/71‑1 (MQ=255)                                       cGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCgc                                  >  1:1896432/1‑67 (MQ=255)                                        ggTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCgg                              <  1:3476465/71‑1 (MQ=255)                                              ggAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAg                        <  1:1799867/71‑1 (MQ=255)                                              ggAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAg                        <  1:1549/71‑1 (MQ=255)                                               gAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAg                        >  1:1664633/1‑70 (MQ=255)                                               gAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCACGCAGAc                                                    >  1:1577394/1‑42 (MQ=38)                                                  cGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCgg                              >  1:701595/1‑61 (MQ=255)                                                   ggTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGAt                                 <  1:3864020/57‑1 (MQ=255)                                                   ggTGTTGTTTTTCACCGCC‑CGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGtt                  <  1:3576836/71‑1 (MQ=255)                                                       ttgttTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGttttt               >  1:2136638/1‑71 (MQ=255)                                                       ttgttTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGttttt               <  1:3119772/71‑1 (MQ=255)                                                       ttgttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGttttt               >  1:494800/1‑71 (MQ=255)                                                        tgttTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTg              >  1:665308/1‑71 (MQ=255)                                                        tgttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCt                            <  1:2767687/57‑1 (MQ=255)                                                        tgttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCt                            <  1:1696605/57‑1 (MQ=255)                                                         gttTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGcc                                    >  1:226049/1‑48 (MQ=255)                                                         gttTTTCACCGCCACGCCTTCCGCAGGCAGACCaaa                                                >  1:3910124/1‑36 (MQ=255)                                                         gttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTg              <  1:258823/70‑1 (MQ=255)                                                         gttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTg              <  1:3747345/70‑1 (MQ=255)                                                          tttttCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGcc            <  1:1173713/71‑1 (MQ=255)                                                            tttCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGccc           <  1:3895304/70‑1 (MQ=255)                                                            tttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGttttt               >  1:2146074/1‑66 (MQ=255)                                                            tttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGcc            >  1:1584778/1‑69 (MQ=255)                                                             ttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGa                                  >  1:3387978/1‑46 (MQ=255)                                                             ttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCa          <  1:500099/70‑1 (MQ=255)                                                                aCCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCa       >  1:950108/1‑70 (MQ=255)                                                                aCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCa       <  1:2444326/70‑1 (MQ=255)                                                                aCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCa       >  1:2085899/1‑70 (MQ=255)                                                                aCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCAt      <  1:2588065/71‑1 (MQ=255)                                                                 ccgccgCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATa     >  1:2843549/1‑71 (MQ=255)                                                                  cgccgcGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCg                                   <  1:3290175/40‑1 (MQ=255)                                                                  cgccgcGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATa     <  1:749804/70‑1 (MQ=255)                                                                  cGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAg         <  1:2346599/66‑1 (MQ=255)                                                                  cGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATAc    >  1:1009014/1‑71 (MQ=255)                                                                   gCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGtt                  <  1:1996108/56‑1 (MQ=255)                                                                    ccACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGc  >  1:1766290/1‑71 (MQ=255)                                                                    ccACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGc  <  1:109626/71‑1 (MQ=255)                                                                      |                                                                     AGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGC  >  minE/412021‑412158

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