Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,636	A→G	21.0%	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	21.0%	77.1 / 17.8	38	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (13/17);  new base G (6/2);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCA  >  minE/1763567‑1763704                                                                      |                                                                     gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                     >  1:814031/1‑71 (MQ=255)   cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                                   <  1:882162/71‑1 (MQ=255)     tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc                                                                 >  1:2131451/1‑71 (MQ=255)     tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc                                                                 >  1:2052093/1‑71 (MQ=255)      ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                                >  1:1995418/1‑71 (MQ=39)         gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                             >  1:1916668/1‑71 (MQ=39)           gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctct                                                           <  1:1765007/71‑1 (MQ=39)              gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg                                                        >  1:98179/1‑71 (MQ=255)              gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                        >  1:1356335/1‑71 (MQ=39)              gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                        >  1:1443316/1‑71 (MQ=39)                     caccaACAAGCTAATCCCATCTGGGCACATCCGATGTCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                                 >  1:1984392/1‑71 (MQ=255)                     caccaACAAGCTAATCCCATCTGGGCACATCCGATGTCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                                 >  1:1981769/1‑71 (MQ=255)                     caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                                 >  1:587150/1‑71 (MQ=39)                      accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                                 >  1:283909/1‑70 (MQ=255)                      accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                                 >  1:2062194/1‑70 (MQ=255)                      accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCg                                                <  1:1548833/71‑1 (MQ=39)                      accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCg                                                <  1:141445/71‑1 (MQ=39)                          acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAcgt                                            <  1:124535/71‑3 (MQ=9)                          acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAcgt                                            <  1:1265733/71‑3 (MQ=9)                          acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                            <  1:101692/71‑1 (MQ=38)                            aaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCttt                                                         >  1:1382007/1‑56 (MQ=9)                                  aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                     <  1:1106766/70‑1 (MQ=35)                                  aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                     <  1:572119/70‑1 (MQ=35)                                          tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                            >  1:2209534/1‑71 (MQ=34)                                             gCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCg                         <  1:943357/71‑1 (MQ=34)                                                  tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                                   >  1:1632594/1‑56 (MQ=35)                                                   ccGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                     >  1:1503505/1‑69 (MQ=35)                                                    cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                  <  1:386310/71‑1 (MQ=34)                                                    cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                  >  1:144758/1‑71 (MQ=34)                                                    cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                  >  1:1773656/1‑71 (MQ=34)                                                    cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                  <  1:745685/71‑1 (MQ=34)                                                    cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATCATGCGGTATTAGCTACCGTTTCCag                  <  1:1279711/71‑1 (MQ=21)                                                      aTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCt                             <  1:91607/58‑1 (MQ=35)                                                         gCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt             <  1:321172/71‑1 (MQ=34)                                                          cAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtagt              <  1:937075/69‑1 (MQ=35)                                                                   cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATc          <  1:1968302/64‑1 (MQ=35)                                                                    gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa  <  1:570195/71‑1 (MQ=34)                                                                    gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa  <  1:806990/71‑1 (MQ=34)                                                                     aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc        >  1:1197091/1‑64 (MQ=35)                                                                      |                                                                     GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCA  >  minE/1763567‑1763704

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