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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,763,556	0	T	C	61.5%	1.3 / 10.4	13	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (2/6);  total (4/9)
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.08e-01

TCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATG  >  Exported/1763506‑1763622                                                   |                                                                   tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                                <  1:177270/71‑1 (MQ=34)           gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                      >  1:118461/1‑71 (MQ=21)              tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                   >  1:272747/1‑71 (MQ=34)                                gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCTTTACCCCACCAACAAGCTAATCCCa                   <  1:297215/69‑1 (MQ=14)                                gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                   <  1:282817/69‑1 (MQ=255)                                 tCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                <  1:176852/71‑1 (MQ=255)                                 tCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                <  1:225392/71‑1 (MQ=255)                                 tCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                <  1:441271/71‑1 (MQ=255)                                   aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg              <  1:499991/71‑1 (MQ=255)                                         tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt        >  1:494258/1‑71 (MQ=255)                                          cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc       >  1:275034/1‑71 (MQ=255)                                           aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt        <  1:369285/69‑1 (MQ=255)                                               cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg  <  1:377373/71‑1 (MQ=255)                                                   |                                                                   TCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATG  >  Exported/1763506‑1763622

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