Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,559	G→T	7.5%	noncoding (285/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,559	0	G	T	7.5%	253.4 / 5.6	79	noncoding (285/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base T (3/3);  total (37/42)
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.28e-01

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAA  >  minE/1763502‑1763626                                                          |                                                                    ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTg                                                        <  1:153424/71‑1 (MQ=11) ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                        >  1:1287326/1‑71 (MQ=32) ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                         <  1:2053499/70‑1 (MQ=14)    gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                      >  1:245794/1‑70 (MQ=35)       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                       >  1:2854884/1‑66 (MQ=35)       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                  >  1:2822562/1‑71 (MQ=34)       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                  >  1:1542170/1‑71 (MQ=21)         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 <  1:1533134/70‑1 (MQ=35)         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 <  1:36093/70‑1 (MQ=35)         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                   <  1:54745/68‑1 (MQ=25)          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               <  1:92239/71‑1 (MQ=34)          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               <  1:2291173/71‑1 (MQ=34)           ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                              <  1:322810/71‑1 (MQ=34)           ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               >  1:34137/1‑70 (MQ=35)           ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               >  1:725313/1‑70 (MQ=35)            cGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                             >  1:2064553/1‑71 (MQ=11)            cGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                             >  1:2234487/1‑71 (MQ=11)             gtgtCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAc                                            >  1:2812637/1‑71 (MQ=11)             gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGAACAGCCAGGGATCGTCGCCtt                                                         >  1:612039/1‑58 (MQ=11)               gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCATACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                            <  1:629313/69‑1 (MQ=255)               gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                          >  1:117713/1‑71 (MQ=34)                tctcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcctc                                         <  1:2114044/71‑3 (MQ=11)                tctcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcctc                                         <  1:1445443/71‑3 (MQ=11)                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               <  1:1206147/65‑1 (MQ=35)                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                          <  1:2317573/70‑1 (MQ=35)                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 >  1:1441911/1‑61 (MQ=35)                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 >  1:2419411/1‑61 (MQ=35)                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                       >  1:1439129/1‑71 (MQ=21)                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                     <  1:2998771/71‑1 (MQ=35)                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                     <  1:805551/71‑1 (MQ=35)                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                      >  1:2588559/1‑70 (MQ=35)                      ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                   <  1:1069538/71‑1 (MQ=39)                       tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                   >  1:553835/1‑70 (MQ=255)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 >  1:567127/1‑55 (MQ=35)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 >  1:1399471/1‑55 (MQ=35)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                  <  1:2970047/70‑1 (MQ=255)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                  <  1:1426133/70‑1 (MQ=255)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                   >  1:708595/1‑69 (MQ=40)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                   >  1:2829712/1‑69 (MQ=40)                        ccAGTGTGGCGGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 <  1:2338003/55‑1 (MQ=255)                         cAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                                <  1:2579144/71‑1 (MQ=255)                         cAGTGTGGCTGGTCATCCTCCCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                                 <  1:1201307/70‑1 (MQ=255)                          aGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                                <  1:1998922/70‑1 (MQ=255)                          aGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                                 <  1:2043415/69‑1 (MQ=255)                             gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                            >  1:2376979/1‑71 (MQ=255)                             gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGAGAGCCGTTACCCCACCAACAAGCTaa                            >  1:857008/1‑71 (MQ=255)                              tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                            <  1:972724/70‑1 (MQ=255)                                gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                         >  1:1797169/1‑71 (MQ=255)                                 cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                      <  1:996444/41‑1 (MQ=38)                                 cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                      <  1:2720448/41‑1 (MQ=38)                                 cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                         >  1:2771500/1‑70 (MQ=255)                                   ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                            <  1:1307396/65‑1 (MQ=255)                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                         <  1:517172/35‑1 (MQ=33)                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                         <  1:1603974/35‑1 (MQ=33)                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                         <  1:2014764/35‑1 (MQ=33)                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                     <  1:112475/71‑1 (MQ=255)                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                      <  1:766141/70‑1 (MQ=255)                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                      <  1:1798062/70‑1 (MQ=255)                                       aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                  >  1:236996/1‑71 (MQ=255)                                          ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca               >  1:2567421/1‑71 (MQ=255)                                          ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                             >  1:1407705/1‑41 (MQ=255)                                          ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                             >  1:2310381/1‑41 (MQ=255)                                           tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                 >  1:148330/1‑68 (MQ=255)                                           tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca               >  1:1789183/1‑70 (MQ=255)                                            ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                             >  1:1160361/1‑55 (MQ=255)                                            ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                             >  1:1535652/1‑55 (MQ=255)                                             tcAGACCAGCTAGGGATCGACGCCTTGGTGAGCCGTTa                                             >  1:2522269/1‑38 (MQ=255)                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc          <  1:677553/71‑1 (MQ=255)                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc          <  1:1231157/71‑1 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt            <  1:828174/69‑1 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc          <  1:1356616/71‑1 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc          >  1:2513018/1‑71 (MQ=255)                                                gACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg         <  1:521077/71‑1 (MQ=255)                                                 aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa        <  1:928956/71‑1 (MQ=255)                                                 aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa        <  1:915099/71‑1 (MQ=255)                                                   cAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                       >  1:1707202/1‑54 (MQ=255)                                                   cAGCTAGGGATCGTCCCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                       >  1:303607/1‑54 (MQ=255)                                                     gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc    <  1:388795/71‑1 (MQ=255)                                                       caGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa  <  1:1840125/70‑1 (MQ=255)                                                          ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg         <  1:1403563/61‑1 (MQ=255)                                                          ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg         <  1:2108278/61‑1 (MQ=255)                                                          |                                                                    GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAA  >  minE/1763502‑1763626

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