Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	2,550,125	G→A	48.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	48.5%	17.8 / 49.9	33	noncoding (1436/1522 nt)	PP_16SF	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (9/8);  new base A (7/9);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

CCCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCA  >  AE015451/2550051‑2550354                                                                           |                                                                                                                                                                                                                                      cccGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAAcgcg                                                                      <  1:235331/236‑1 (MQ=11)        cTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGAccc                                                                  <  2:580287/233‑1 (MQ=14)             acacacCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAg                                                          >  1:6764/1‑236 (MQ=31)                    gCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGGATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                   >  1:595194/1‑236 (MQ=14)                     cccGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAACGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                   >  1:442663/1‑235 (MQ=14)                        gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTc                                               >  1:733326/1‑236 (MQ=31)                        gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTc                                               >  1:87487/1‑236 (MQ=14)                         tCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTACAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCt                                                >  2:717910/1‑234 (MQ=14)                               cATGGGAGTGGGTTGCACCAGCAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                     <  2:286530/207‑1 (MQ=11)                               cATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTgg                                         <  1:699180/235‑1 (MQ=31)                                aTGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCt                                                <  2:89294/227‑1 (MQ=31)                                   ggAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGt                                        >  2:645303/1‑232 (MQ=31)                                       tGGGTTGCACCAGGAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTACGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgc                                  <  2:405454/234‑1 (MQ=14)                                       tGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                <  1:244445/236‑1 (MQ=31)                                       tGGGTTGCACCAGAAGTAGCTAGTCTAACCGTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTgg                                         <  2:273704/227‑1 (MQ=14)                                         ggTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAc                              <  1:138301/236‑1 (MQ=31)                                          gTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgc                                  <  2:383981/231‑1 (MQ=31)                                          gTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAcc                             >  1:612882/1‑236 (MQ=31)                                            tGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAcccc                           >  1:603849/1‑236 (MQ=31)                                            tGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                     >  1:625353/1‑194 (MQ=14)                                                    aaGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                    <  2:529457/219‑1 (MQ=31)                                                     aGTAGCTAGTCGAACCGTCGGGGGGACGGTTACCAAGGTGGGAGTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                     <  2:678886/185‑1 (MQ=9)                                                      gTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAggg                  >  1:228881/1‑235 (MQ=31)                                                      gTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGt                 >  1:646217/1‑236 (MQ=14)                                                                 aaCCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTta                                                             >  2:486431/1‑181 (MQ=31)                                                                 aaCCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCa       >  1:245031/1‑235 (MQ=31)                                                                 aaCCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                >  2:675333/1‑210 (MQ=14)                                                                 aaCCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                <  1:675334/210‑1 (MQ=14)                                                                 aaCCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAg      >  1:115912/1‑236 (MQ=14)                                                                    cTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCt                          <  2:87487/213‑1 (MQ=14)                                                                     ttCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCa  <  1:527646/236‑1 (MQ=14)                                                                     tgcGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGGCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCGCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGt            <  2:131906/224‑1 (MQ=14)                                                                       cGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAc                              <  2:255033/206‑1 (MQ=31)                                                                        gggaggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                    <  2:577725/199‑1 (MQ=14)                                                                           |                                                                                                                                                                                                                                      CCCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCA  >  AE015451/2550051‑2550354

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