Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,211,696	G→A	49.2%	A461V (GCG→GTG)	rng ←	ribonuclease G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,211,696	0	G	A	49.2%	3.5 / 96.4	67	A461V (GCG→GTG)	rng	ribonuclease G
Reads supporting (aligned to +/- strand):  ref base G (15/19);  new base A (17/16);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGACCAGGAAACG  >  minE/2211626‑2211763                                                                       |                                                                    aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCg                                                                     >  1:591088/1‑71 (MQ=255) aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCg                                                                     <  1:3735008/71‑1 (MQ=255)       cTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCg                                                               <  1:1043219/71‑1 (MQ=255)          cTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAgt                                                            >  1:1119730/1‑71 (MQ=255)          cTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAg                                                             <  1:2584879/70‑1 (MQ=255)          cTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAg                                                             <  1:1032270/70‑1 (MQ=255)            ggTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtg                                                           >  1:4081081/1‑70 (MQ=255)             gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGa                                                              >  1:3935574/1‑66 (MQ=255)             gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtgt                                                          <  1:3620277/70‑1 (MQ=255)             gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAgtgt                                                          <  1:416797/70‑1 (MQ=255)             gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAgtgt                                                          <  1:2127194/70‑1 (MQ=255)              ttATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGa                                                        <  1:3540786/71‑1 (MQ=255)                   gagCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGActct                                                    <  1:1128966/70‑1 (MQ=255)                         ttCAATTTGTACTTTAACCTGTTTCCCAACGAAAATTTCCACTTCCGCCAGCGAg                                                             <  1:1039493/55‑1 (MQ=255)                          tCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCttt                                            >  1:1579715/1‑71 (MQ=255)                           cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGa                                                              >  1:3640704/1‑52 (MQ=255)                           cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGa                                                        >  1:3806084/1‑58 (MQ=255)                            aaTTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCa                                          <  1:956214/71‑1 (MQ=255)                            aaTTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTc                                           <  1:2204031/70‑1 (MQ=255)                             aTTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCaa                                         <  1:3434934/71‑1 (MQ=255)                              tttGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAg                                        <  1:2116756/71‑1 (MQ=255)                              tttGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCaa                                         <  1:2370883/70‑1 (MQ=255)                                tGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGc                                      >  1:2921434/1‑71 (MQ=255)                                tGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGc                                      >  1:3705801/1‑71 (MQ=255)                                tGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGc                                      >  1:271771/1‑71 (MQ=255)                                 gTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGc                                      <  1:1540953/70‑1 (MQ=255)                                  tACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCt                                     <  1:1687181/70‑1 (MQ=255)                                  tACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCt                                     <  1:2301317/70‑1 (MQ=255)                                    cTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGa                                                              <  1:344302/43‑1 (MQ=255)                                    cTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAgtgt                                                          <  1:865703/47‑1 (MQ=255)                                    cTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAgtgt                                                          <  1:96599/47‑1 (MQ=255)                                      ttAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                                >  1:2192593/1‑71 (MQ=255)                                       tAACCTGTTTTCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGa                                                        >  1:408553/1‑46 (MQ=39)                                       tAACCTGTTTGCCAACGAAAATTTCCACTTCCGcc                                                                   >  1:43150/1‑35 (MQ=255)                                       tAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGa                                                        >  1:2150680/1‑46 (MQ=255)                                       tAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGa                                                        >  1:3568634/1‑46 (MQ=255)                                       tAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                                <  1:3123230/70‑1 (MQ=255)                                       tAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                                >  1:1913984/1‑70 (MQ=255)                                       tAACCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTACGCCTTTCAAGGCTTCAGc                                >  1:4042986/1‑70 (MQ=255)                                        aaCCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                                <  1:3398539/69‑1 (MQ=255)                                         aCCTGTTTGCCAACGAAAATTTCCACTTCCACCAGCg                                                               <  1:2551887/37‑1 (MQ=255)                                           cTGTTTGCCAACGAAAATTTCCACTTCCACCAGCGAg                                                             >  1:2566258/1‑37 (MQ=255)                                             gTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                                <  1:2499719/64‑1 (MQ=255)                                              tttGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCg                        >  1:3130779/1‑71 (MQ=255)                                              tttGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCg                        >  1:1984506/1‑71 (MQ=255)                                                tGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCt                                     <  1:2724555/56‑1 (MQ=255)                                                tGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCg                        >  1:3661870/1‑69 (MQ=255)                                                    aaCGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGgaga                    <  1:3379739/69‑1 (MQ=255)                                                    aaCGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg                  >  1:2362376/1‑71 (MQ=255)                                                    aaCGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg                  <  1:1618658/71‑1 (MQ=255)                                                    aaCGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg                  <  1:489813/71‑1 (MQ=255)                                                    aaCGAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg                  <  1:1386518/71‑1 (MQ=255)                                                      cGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa                >  1:1600328/1‑71 (MQ=255)                                                      cGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa                <  1:1207308/71‑1 (MQ=255)                                                       gAAAATTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa                >  1:1248744/1‑70 (MQ=255)                                                         aaaTTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCaa                                         >  1:1139870/1‑43 (MQ=255)                                                         aaaTTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa                >  1:333771/1‑68 (MQ=255)                                                          aaTTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGa            >  1:3880962/1‑71 (MQ=255)                                                           aTTTCCACTTCCGCCAGCTAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGAc           <  1:2644716/71‑1 (MQ=255)                                                           aTTTCCACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGAc           >  1:3242071/1‑71 (MQ=255)                                                            tttCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGAcc          <  1:3843078/71‑1 (MQ=255)                                                              tCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGACCAg        >  1:1688126/1‑71 (MQ=255)                                                               ccACTTCCGCCAGCGAGTGTGAGTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGACCAg        <  1:867098/70‑1 (MQ=255)                                                               ccACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGACCAg        <  1:2073138/70‑1 (MQ=255)                                                               ccACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGACCAg        >  1:3740334/1‑70 (MQ=255)                                                               ccACTTCCACCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGACCAg        >  1:3613923/1‑70 (MQ=255)                                                                 aCTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCaa                                         >  1:2592607/1‑35 (MQ=255)                                                                   ttCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                                <  1:497404/42‑1 (MQ=255)                                                                    tCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGACCAGGAAACg  >  1:1532017/1‑71 (MQ=255)                                                                       |                                                                    AAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATAGACCAGGAAACG  >  minE/2211626‑2211763

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