Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,636	A→G	32.3%	noncoding (208/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,636	0	A	G	32.3%	74.5 / 57.3	62	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (19/23);  new base G (13/7);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATC  >  minE/1763568‑1763706                                                                     |                                                                       ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAggg                                                                      >  1:2609408/1‑71 (MQ=255) ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAggg                                                                      <  1:3434099/71‑1 (MQ=255)   ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                        >  1:2561156/1‑67 (MQ=255)   ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTc                                                                    >  1:1120632/1‑71 (MQ=255)   ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTc                                                                    >  1:2343415/1‑71 (MQ=255)     ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccc                                                                  >  1:883299/1‑71 (MQ=255)       tGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccccc                                                                >  1:2579289/1‑71 (MQ=39)          gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctct                                                             <  1:1693343/71‑1 (MQ=255)           ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCtt                                                            <  1:2496027/71‑1 (MQ=255)            cGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCATCttt                                                           <  1:1542556/71‑1 (MQ=255)              ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                         >  1:1644727/1‑71 (MQ=255)              ttACCCCACCAACAAGCTAATCCCATCTGGCCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                         <  1:1668699/71‑1 (MQ=255)                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAggg                                                                      >  1:443954/1‑56 (MQ=255)                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                       >  1:2337939/1‑71 (MQ=255)                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                       >  1:1870451/1‑71 (MQ=255)                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc                                                       >  1:474206/1‑71 (MQ=39)                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc                                                       >  1:1344829/1‑71 (MQ=39)                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc                                                       >  1:4032991/1‑71 (MQ=39)                   ccaccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt                                                           <  1:769025/64‑1 (MQ=40)                    caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCct                                                               >  1:2928957/1‑59 (MQ=255)                    caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                                   >  1:1523939/1‑71 (MQ=255)                    caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                                   >  1:1332465/1‑71 (MQ=255)                    caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                                   <  1:2217629/71‑1 (MQ=39)                    caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                                   <  1:1268480/71‑1 (MQ=39)                     accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCg                                                  <  1:2210415/71‑1 (MQ=255)                     accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTCGTCTTGCg                                                  <  1:4016895/71‑1 (MQ=255)                     accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCg                                                  <  1:3427728/71‑1 (MQ=39)                     accaACAAGCTAATCCCATCTGGGCACATCAGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCg                                                  <  1:306879/71‑1 (MQ=14)                      ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCg                                                  <  1:673425/70‑1 (MQ=255)                      ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                 >  1:2144367/1‑71 (MQ=39)                      ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                 >  1:2353055/1‑71 (MQ=39)                      ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                 <  1:1799389/71‑1 (MQ=39)                      ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                 >  1:3558664/1‑71 (MQ=39)                      ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGACCGAGGGTCCCCCTCTTTGGt                                                        <  1:945652/64‑1 (MQ=255)                       caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                                   >  1:1024658/1‑68 (MQ=255)                       caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                                   >  1:3144472/1‑68 (MQ=255)                       caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                         <  1:1630678/62‑1 (MQ=39)                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAcg                                               >  1:137927/1‑70 (MQ=19)                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAcg                                               >  1:1563947/1‑70 (MQ=19)                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcccc                                                                 >  1:2275145/1‑53 (MQ=38)                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcccc                                                                 >  1:823108/1‑53 (MQ=38)                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcccc                                                                 >  1:1740630/1‑53 (MQ=38)                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                               <  1:2173764/71‑1 (MQ=38)                         acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                              <  1:2532351/71‑1 (MQ=38)                         acaaGCTAATCCCATCTGGGCACATCAGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                              >  1:473290/1‑71 (MQ=255)                            aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                           >  1:906101/1‑71 (MQ=34)                                    ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                    >  1:3704793/1‑70 (MQ=35)                                    ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                    <  1:876313/70‑1 (MQ=35)                                    ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                    >  1:3573391/1‑70 (MQ=35)                                     cATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                               <  1:1658123/58‑1 (MQ=31)                                     cATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                                  <  1:1332739/71‑1 (MQ=34)                                        cTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                                >  1:1464192/1‑70 (MQ=35)                                           ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                                <  1:1619574/67‑1 (MQ=35)                                           ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAcc                            <  1:2035421/71‑1 (MQ=34)                                           ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAc                             <  1:874701/70‑1 (MQ=35)                                            gCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATTCGGTATTAGCt                               <  1:3881550/67‑1 (MQ=25)                                                aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                       <  1:1616973/71‑1 (MQ=34)                                                aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                       >  1:1753656/1‑71 (MQ=34)                                                 tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                              >  1:1774394/1‑63 (MQ=35)                                                 tCCGATGGCAAGAGGCCCGAAGGTCACCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                       <  1:20491/70‑1 (MQ=255)                                                            gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcc           >  1:76361/1‑71 (MQ=34)                                                              ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc         >  1:1229870/1‑71 (MQ=34)                                                                 ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGt                                      <  1:1063853/39‑1 (MQ=40)                                                                 ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                                <  1:3693295/45‑1 (MQ=38)                                                                   gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa    <  1:1219268/71‑1 (MQ=34)                                                                   gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa    <  1:2764752/71‑1 (MQ=34)                                                                     aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATc  <  1:3826820/71‑1 (MQ=34)                                                                     aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATc  <  1:3495279/71‑1 (MQ=34)                                                                     |                                                                       CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATC  >  minE/1763568‑1763706

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