Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,556	T→C	46.4%	noncoding (288/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,763,556	0	T	C	46.4%	17.1 / 125.3	84	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (21/24);  new base C (18/21);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCC  >  minE/1763502‑1763618                                                       |                                                               ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                 <  1:177386/70‑1 (MQ=31)  gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                 <  1:1992637/69‑1 (MQ=14)  gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                               <  1:2632081/71‑1 (MQ=21)       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                              >  1:913719/1‑67 (MQ=35)       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                          >  1:616403/1‑71 (MQ=34)         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAg                                           <  1:991844/68‑1 (MQ=35)         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                         <  1:1216035/70‑1 (MQ=35)         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCttt                                                <  1:1272175/63‑2 (MQ=9)         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                         <  1:1163532/70‑1 (MQ=25)         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCg                                        <  1:2450304/71‑1 (MQ=21)          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                       <  1:959387/71‑1 (MQ=21)          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                       >  1:376287/1‑71 (MQ=21)          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                       <  1:2486551/71‑1 (MQ=21)           ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                      <  1:2542348/71‑1 (MQ=34)            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                                      <  1:2776732/70‑1 (MQ=25)            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:1601925/1‑71 (MQ=21)              tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                            <  1:1694624/62‑1 (MQ=35)              tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                   >  1:54480/1‑71 (MQ=34)                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                    <  1:2356271/68‑1 (MQ=35)                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                  <  1:708469/70‑1 (MQ=25)                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:2787752/1‑65 (MQ=35)                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAcccaa                                >  1:2815467/1‑68 (MQ=255)                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                               >  1:2728964/1‑71 (MQ=21)                   cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                               >  1:818269/1‑70 (MQ=25)                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                              >  1:1147651/1‑70 (MQ=35)                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                               <  1:1394047/69‑1 (MQ=25)                      ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGGCtt                                                 <  1:2437610/49‑1 (MQ=11)                      ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                           <  1:2282424/71‑1 (MQ=255)                       tCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcc                                   >  1:301976/1‑62 (MQ=11)                       tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                           >  1:1286160/1‑70 (MQ=40)                       tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGTGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                           <  1:413915/70‑1 (MQ=255)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:2114655/1‑59 (MQ=35)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:2117972/1‑59 (MQ=35)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:662506/1‑59 (MQ=35)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                          <  1:1465701/70‑1 (MQ=255)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                           >  1:1067624/1‑69 (MQ=40)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                          <  1:2243526/70‑1 (MQ=255)                        acAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                 <  1:1315999/46‑1 (MQ=17)                         cAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                         <  1:610281/70‑1 (MQ=255)                         cAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                         <  1:312777/70‑1 (MQ=255)                           gtgtGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                      <  1:2049427/71‑1 (MQ=255)                           gtgtGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                      <  1:178920/71‑1 (MQ=255)                             gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                     >  1:883474/1‑70 (MQ=255)                              tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                   >  1:2003415/1‑71 (MQ=255)                              tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                <  1:2796566/42‑1 (MQ=255)                               ggCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                  >  1:1700434/1‑71 (MQ=255)                                gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                              <  1:2495918/42‑1 (MQ=38)                                gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                >  1:2382979/1‑40 (MQ=255)                                gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                              >  1:266420/1‑42 (MQ=255)                                gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGa                                            >  1:1647427/1‑44 (MQ=255)                                 cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                       <  1:2788880/48‑1 (MQ=38)                                  tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                         <  1:1596256/61‑1 (MQ=255)                                   ggTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGa                                            >  1:551120/1‑41 (MQ=14)                                   ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                            >  1:2331714/1‑41 (MQ=38)                                     tcgtCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTa                     <  1:511182/59‑1 (MQ=11)                                     tcgtCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTa                     <  1:83927/59‑1 (MQ=11)                                     tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt            <  1:2197121/71‑1 (MQ=255)                                      cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg           >  1:2559216/1‑71 (MQ=255)                                       aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg           <  1:2016289/70‑1 (MQ=255)                                       aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg           <  1:2178486/70‑1 (MQ=255)                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                       >  1:2258991/1‑42 (MQ=255)                                        tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg         >  1:517409/1‑71 (MQ=255)                                        tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                >  1:2208264/1‑64 (MQ=255)                                        tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt            <  1:1259189/68‑1 (MQ=255)                                        tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg         >  1:945292/1‑71 (MQ=255)                                         ccTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt              >  1:702954/1‑65 (MQ=255)                                         ccTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:248404/1‑42 (MQ=255)                                          ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca       <  1:198390/71‑1 (MQ=255)                                          ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca       <  1:1970787/71‑1 (MQ=255)                                          ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg          >  1:452041/1‑68 (MQ=255)                                            ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca     >  1:2199794/1‑71 (MQ=255)                                            ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                  >  1:352851/1‑58 (MQ=255)                                             tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                   <  1:2798503/56‑1 (MQ=255)                                             tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                   <  1:1531413/56‑1 (MQ=255)                                             tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                <  1:879338/59‑1 (MQ=255)                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTccc                                   >  1:2749497/1‑38 (MQ=255)                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc  <  1:2799744/71‑1 (MQ=255)                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc  >  1:2035393/1‑71 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt              >  1:480009/1‑59 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt              >  1:2868278/1‑59 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca     <  1:2834625/68‑1 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc  <  1:1216708/71‑1 (MQ=255)                                                gACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                <  1:1939678/56‑1 (MQ=255)                                                 aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg           >  1:1946926/1‑60 (MQ=255)                                                       tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                    >  1:2890464/1‑45 (MQ=255)                                                       |                                                               GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCC  >  minE/1763502‑1763618

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