Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,636	A→G	14.6%	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	14.6%	95.0 / 9.6	41	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (11/24);  new base G (2/4);  total (13/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCC  >  minE/1763569‑1763699                                                                    |                                                                cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCTAGGGt                                                              >  1:1028553/1‑71 (MQ=255) cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGt                                                              >  1:4804307/1‑71 (MQ=255)  ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                             >  1:1251940/1‑71 (MQ=255)  ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                             >  1:4484325/1‑71 (MQ=255)        aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCct                                                        <  1:304188/70‑1 (MQ=255)        aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                       <  1:2904722/71‑1 (MQ=39)        aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                       <  1:2615005/71‑1 (MQ=39)        aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                       <  1:2200549/71‑1 (MQ=39)             ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                  <  1:1494854/71‑1 (MQ=255)                ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt                                               <  1:1016439/71‑1 (MQ=255)                ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                               <  1:2987580/71‑1 (MQ=39)                ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                               <  1:4917521/71‑1 (MQ=39)                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                  <  1:478146/65‑1 (MQ=255)                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                            <  1:4913370/71‑1 (MQ=39)                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                 >  1:3341024/1‑48 (MQ=255)                      caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                          <  1:3641603/70‑1 (MQ=39)                      caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAAGCCCGAAGGTCCCCCTCTTTg                                                   <  1:3287513/61‑1 (MQ=255)                      caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAAGCCCGAAGGTCCCCCTCTTTg                                                   <  1:2732525/61‑1 (MQ=255)                      caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAAGCCCGAAGGTCCCCCTCTTTg                                                   <  1:4468679/61‑1 (MQ=255)                       aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                          <  1:3915921/69‑1 (MQ=37)                       aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAAGCCCGAAGGTCCCCCTCTTTg                                                   <  1:3690255/60‑1 (MQ=255)                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                              >  1:419489/1‑48 (MQ=38)                           aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:4272631/71‑1 (MQ=34)                           aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:3532497/71‑1 (MQ=34)                           aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:739830/71‑1 (MQ=34)                           aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:1495861/71‑1 (MQ=34)                           aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:1768344/71‑1 (MQ=34)                            gCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATg                                   >  1:3521586/1‑71 (MQ=34)                            gCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATg                                   <  1:2947079/71‑1 (MQ=34)                            gCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATg                                   >  1:1528437/1‑71 (MQ=34)                                  cccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                        >  1:488861/1‑60 (MQ=31)                                  cccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                        >  1:4744069/1‑60 (MQ=31)                                          ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                         <  1:2985546/67‑1 (MQ=35)                                                   gATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt            <  1:3109179/71‑1 (MQ=34)                                                         aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt      <  1:4078482/71‑1 (MQ=34)                                                            agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc   >  1:3509991/1‑71 (MQ=34)                                                            agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc   >  1:1060731/1‑71 (MQ=34)                                                            agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc   >  1:1997506/1‑71 (MQ=34)                                                            agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc   >  1:3879270/1‑71 (MQ=34)                                                             ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc  <  1:4318533/71‑1 (MQ=34)                                                             ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc  <  1:463664/71‑1 (MQ=34)                                                             ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc  <  1:1430355/71‑1 (MQ=34)                                                                    aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              >  1:2229688/1‑36 (MQ=40)                                                                    aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              >  1:1392609/1‑36 (MQ=40)                                                                    aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              >  1:1178561/1‑36 (MQ=40)                                                                    aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              >  1:3202234/1‑36 (MQ=40)                                                                    aGGTCCCCCTCTTTGGTCTTGCGACATTATGCAGt                               >  1:161583/1‑35 (MQ=255)                                                                    |                                                                CTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCC  >  minE/1763569‑1763699

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