Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	527,917	A→G	16.2%	noncoding (1152/2902 nt)	PP_23SC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	527,917	0	A	G	16.2%	83.2 / 10.7	37	noncoding (1152/2902 nt)	PP_23SC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (13/18);  new base G (2/4);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

GGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAAC  >  NC_002947/527811‑528044                                                                                                           |                                                                                                                                   ggatgTTGGCTTAGAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTaga                                                                                                    >  1:327740/1‑136 (MQ=255) ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTaga                                                                                                    >  1:369383/1‑136 (MQ=17) ggaggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTaga                                                                                                    >  2:30437/1‑136 (MQ=32)  gaggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTagag                                                                                                   <  2:331453/136‑1 (MQ=32)    ggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAgagga                                                                                                 <  1:159849/136‑1 (MQ=32)    ggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAgagga                                                                                                 <  2:389212/136‑1 (MQ=32)    ggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAgagga                                                                                                 <  2:204747/136‑1 (MQ=32)         gCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGtt                                                                                            >  1:492982/1‑136 (MQ=32)           ttAGAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCt                                                                                          >  1:554196/1‑136 (MQ=12)                                   aagCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTga                                                                   >  1:442270/1‑136 (MQ=0)                                      cGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGaa                                                                >  1:427453/1‑136 (MQ=0)                                        tAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCt                                                             >  1:367397/1‑136 (MQ=32)                                          aTAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTg                                                           <  1:30437/136‑1 (MQ=32)                                          aTAGATCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCttt                                                           <  2:226648/136‑2 (MQ=17)                                                   tAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTa                                                  >  1:253086/1‑136 (MQ=32)                                                   tAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTa                                                  >  1:379649/1‑136 (MQ=32)                                                   tAGTCGAGTCGGCCTGCGCGGAAGAAGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATCCGGTAGAGGAGCGTTCTGAAAGCCTCTGAAGGTGAGTTGAGAAGCTTCCTGGAGGTa                                                  >  1:255302/1‑136 (MQ=11)                                                                   cgcgGATGATGTAACGGCGCTCAAACCAAGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGc                                  <  2:255302/136‑1 (MQ=12)                                                                   cgcgGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGc                                  <  2:253086/136‑1 (MQ=32)                                                                       gAAGATGTAACGGGGCTCAATCCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCTGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCTGAAGTGCGAATGCTGag                              <  2:136246/136‑2 (MQ=12)                                                                       gAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                              <  2:18543/136‑1 (MQ=32)                                                                       gAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                              <  2:379649/136‑1 (MQ=32)                                                                        aaGATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                              <  2:111498/136‑1 (MQ=0)                                                                         aGATGTAACGGGGCTCAATCCATGCACCTAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                            <  2:269160/136‑1 (MQ=17)                                                                         aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                            <  2:492982/136‑1 (MQ=32)                                                                          taTGTAACGGGGCTCAAACCATGCACCGAAGATACGTGTATCA‑TCTTA‑‑CGATGATTCTGTAGAGGATCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTTCGAATGCTGACATg                           <  2:535134/135‑1 (MQ=11)                                                                          gATGTAACGGGGCTCAAACCTTGCACCGAAGGTACTGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATg                           <  2:277657/136‑1 (MQ=12)                                                                          gATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATg                           <  2:554196/136‑1 (MQ=32)                                                                          gATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                            <  2:200825/136‑1 (MQ=0)                                                                          gATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                            <  2:333287/136‑1 (MQ=0)                                                                            tGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGa                          <  2:442270/136‑1 (MQ=0)                                                                                   gggCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGAc                  <  2:361113/136‑1 (MQ=32)                                                                                       tCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATg              >  2:413276/1‑136 (MQ=32)                                                                                       tCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATg              >  1:392871/1‑136 (MQ=32)                                                                                        cAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGc             >  2:445526/1‑136 (MQ=32)                                                                                             cATGCACCGATGCTACGGTTATCA‑TATTA‑‑CGATTATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCATAAGTGCTAATGCTGACATGAGTAACGACAATGCGAGTg        <  2:369383/136‑1 (MQ=11)                                                                                                   ccGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACAACAATGCGAGTGAAAAAc  >  1:491401/1‑136 (MQ=17)                                                                                                           |                                                                                                                                   GGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAAC  >  NC_002947/527811‑528044

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑