Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,211,683	A→G	30.1%	I465I (ATT→ATC)	rng ←	ribonuclease G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,211,683	0	A	G	30.1%	99.0 / 61.7	73	I465I (ATT→ATC)	rng	ribonuclease G
Reads supporting (aligned to +/- strand):  ref base A (29/22);  new base G (9/13);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

CATTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCA  >  minE/2211614‑2211749                                                                      |                                                                   catTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGaaaa                                                                   >  1:3430115/1‑71 (MQ=255)    tACGACGTCAAACTGCTCCTTGTTATAGAGCGGTTCAATTTGTACTTTATCCTGTTTGCCAACGAAAAtt                                                                 <  1:2953527/70‑1 (MQ=255)      cgacgTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTc                                                               <  1:3363420/70‑1 (MQ=255)      cgacgTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTc                                                               <  1:192059/70‑1 (MQ=255)      cgacgTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTc                                                               <  1:1887999/70‑1 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTcc                                                        >  1:2670330/1‑70 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCg                                                       >  1:1422034/1‑71 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCAc                                                            >  1:2437360/1‑66 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCg                                                       >  1:2143359/1‑71 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCg                                                       >  1:3202756/1‑71 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCg                                                       >  1:2697094/1‑71 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCg                                                       >  1:1857065/1‑71 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCg                                                       >  1:563058/1‑71 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCg                                                       >  1:1202212/1‑71 (MQ=255)             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCAAACGAAAATTTCCAc                                                            >  1:101870/1‑66 (MQ=255)                 tGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACt                                                           >  1:164476/1‑63 (MQ=255)                      cTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAgt                                              <  1:1549957/71‑1 (MQ=255)                         gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAgtgt                                            <  1:3415245/70‑1 (MQ=255)                         gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAgtgt                                            <  1:1999469/70‑1 (MQ=255)                         gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtgt                                            <  1:3395430/70‑1 (MQ=255)                         gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtgt                                            <  1:1737836/70‑1 (MQ=255)                          ttATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGa                                          <  1:191573/71‑1 (MQ=255)                          ttATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGa                                          <  1:551202/71‑1 (MQ=255)                          ttATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGa                                          <  1:817869/71‑1 (MQ=255)                          ttATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGa                                          <  1:3797660/71‑1 (MQ=255)                          ttATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGa                                          <  1:1503825/71‑1 (MQ=255)                           tataGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGAc                                         >  1:69036/1‑71 (MQ=255)                             taGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGActc                                       <  1:252970/71‑1 (MQ=255)                                  cGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGc                                  <  1:4059848/71‑1 (MQ=255)                                   ggTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGc                                  <  1:27116/70‑1 (MQ=255)                                    gTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCAc                                                            >  1:1926247/1‑43 (MQ=255)                                     ttCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTc                                    <  1:3171595/66‑1 (MQ=255)                                      tCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:2233714/1‑71 (MQ=255)                                      tCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:2504605/1‑71 (MQ=255)                                       cAATTTGTACTTTAACCTGTTTGGCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:2995046/1‑70 (MQ=255)                                       cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTc                                                               <  1:420357/37‑1 (MQ=255)                                       cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:3151792/1‑70 (MQ=255)                                       cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:3397673/1‑70 (MQ=255)                                       cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:167231/1‑70 (MQ=255)                                       cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:3469630/1‑70 (MQ=255)                                       cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:1276480/1‑70 (MQ=255)                                       cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:1844630/1‑70 (MQ=255)                                       cAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCAGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:2350655/1‑70 (MQ=255)                                       cAATTTGTACTTTAACATGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCttt                              >  1:3819617/1‑70 (MQ=255)                                        aaTTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTc                             >  1:2352695/1‑70 (MQ=255)                                        aaTTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTc                             >  1:3555646/1‑70 (MQ=255)                                        aaTTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCa                            <  1:783475/71‑1 (MQ=255)                                        aaTTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCa                            <  1:3048326/71‑1 (MQ=255)                                          tttGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAg                          <  1:2026746/71‑1 (MQ=255)                                          tttGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAg                          <  1:1953538/71‑1 (MQ=255)                                          tttGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAg                          <  1:1655095/71‑1 (MQ=255)                                          tttGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAg                          <  1:1113610/71‑1 (MQ=255)                                                  ttAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                  >  1:271526/1‑71 (MQ=255)                                                  ttAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                  >  1:161915/1‑71 (MQ=255)                                                  ttAACCTCTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                  >  1:721732/1‑71 (MQ=255)                                                   tAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                  <  1:880810/70‑1 (MQ=255)                                                   tAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                  >  1:3399300/1‑70 (MQ=255)                                                   tAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                  <  1:2674482/70‑1 (MQ=255)                                                   tAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                  >  1:2654493/1‑70 (MQ=255)                                                   tAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                  >  1:1166305/1‑70 (MQ=255)                                                     aCCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCt                 <  1:1769678/69‑1 (MQ=255)                                                          tttGCCAACGAAGATTTCCACTTCCGCCAGCGAGtgtg                                           <  1:3692679/38‑1 (MQ=255)                                                            tGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGa                                          >  1:1299917/1‑37 (MQ=255)                                                                aaCGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg    <  1:3175612/71‑1 (MQ=255)                                                                aaCGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg    <  1:3933790/71‑1 (MQ=255)                                                                 aCGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg    >  1:2465372/1‑70 (MQ=255)                                                                 aCGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg    >  1:2652383/1‑70 (MQ=255)                                                                 aCGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg    >  1:1727996/1‑70 (MQ=255)                                                                 aCGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg    >  1:1269678/1‑70 (MQ=255)                                                                 aCGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg    >  1:1212766/1‑70 (MQ=255)                                                                 aCGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGc   <  1:569410/71‑1 (MQ=255)                                                                 aCGAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGc   <  1:75753/71‑1 (MQ=255)                                                                  cGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa  <  1:2066187/71‑1 (MQ=255)                                                                  cGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa  <  1:3910212/71‑1 (MQ=255)                                                                    aaaaTTTCCACTTCCGCCAGCGAGTGTGACTCATCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa  >  1:45630/1‑69 (MQ=255)                                                                      |                                                                   CATTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCA  >  minE/2211614‑2211749

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