Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	2,550,125	G→A	55.5%	noncoding (1436/1522 nt)	PP_16SF →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	2,550,125	0	G	A	55.5%	‑2.7 / 21.5	18	noncoding (1436/1522 nt)	PP_16SF	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base A (3/7);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

CCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCT  >  AE015451/2550052‑2550359                                                                          |                                                                                                                                                                                                                                           ccGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCGTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAAcgcg                                                                           <  1:270479/235‑1 (MQ=11) acGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                          <  1:325793/235‑1 (MQ=11)   gggCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGAcc                                                                        >  1:235381/1‑236 (MQ=31)               cacCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCCGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTg                                                                     <  2:170601/227‑1 (MQ=14)                     ccGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                        <  2:235377/234‑1 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTc                                                    >  1:509577/1‑236 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTc                                                    >  1:465758/1‑236 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTc                                                    >  1:385423/1‑236 (MQ=14)                        tCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCa                                                   >  1:441657/1‑236 (MQ=14)                                  ggAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  2:385419/236‑1 (MQ=14)                                      tGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                     <  1:172507/236‑1 (MQ=14)                                           tGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCta                                                    >  1:343703/1‑215 (MQ=14)                                             cACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcg                                      >  1:342352/1‑228 (MQ=14)                                             cACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcg                                      <  2:342348/228‑1 (MQ=14)                                              aCCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGttt                                               >  2:275458/1‑217 (MQ=14)                                               ccAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                        <  1:278117/208‑1 (MQ=14)                                                         cTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAc                                   <  1:153302/219‑1 (MQ=14)                                                                      cGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         <  2:169880/200‑1 (MQ=31)                                                                       ggggggACGGTTACCACGCGGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                        <  2:343699/184‑1 (MQ=14)                                                                         ggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                         >  1:264183/1‑197 (MQ=31)                                                                         ggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCt  >  1:40163/1‑236 (MQ=31)                                                                         ggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCt  <  1:372882/236‑1 (MQ=31)                                                                          |                                                                                                                                                                                                                                           CCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCT  >  AE015451/2550052‑2550359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑