Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,313,359	A→C	48.3%	noncoding (1293/1522 nt)	PP_16SG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,313,359	0	A	C	48.3%	‑4.4 / 36.6	29	noncoding (1293/1522 nt)	PP_16SG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (4/11);  new base C (9/5);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 6.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTCTGGTGCAACCCACTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGACATTCTGATTCGCGATTACTAGCGATTCCGACTTCACGCAGTCGAGTTGCAGACTGCGATCCGGACTACGATCGGTTTTGTGAGATTAGCTCCACCTCGCG‑GCTTGGCAACCCTCTGTACCGACCATTGTAGCACGTGTGTAGCCCAGGCCGTAAGGGCCATGATGACTTGACGTCATCCCCACCTTCC  >  NC_002947/5313236‑5313495                                                                                                                            |                                                                                                                                          acgtCTGGTGCAACCCACTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGACATTCTGATTCGCGATTACTAGCGATTCCGACTTCACGCAGTCGAGTTGCAGACTGCGATCCGGACTACGATCGGtttt                                                                                                                 <  1:6660546/147‑1 (MQ=1)       ggcagaaCCAACTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGAt                                                                                                          <  2:6622353/145‑1 (MQ=1)       ggcagaaCCAACTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGAt                                                                                                          <  2:1392014/145‑1 (MQ=1)              ccaaCTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GATATTAGCTca                                                                                                   >  2:6004546/4‑149 (MQ=1)              ccaaCTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTca                                                                                                   >  1:4844414/4‑149 (MQ=1)              ccaaCTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTca                                                                                                   <  1:6120392/147‑2 (MQ=1)              ccaaCTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTca                                                                                                   >  1:6471878/4‑149 (MQ=1)                           tgagACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGACATTCTGATTCGCGATTACTAGCGATTCCGACTTCACGCAGTCGAGTTGCAGACTGCGATCCGGACTACGATCGGTTTTGTGAGATTAGCTCCACCTCGCG‑GCtt                                                                                      <  1:8477859/147‑1 (MQ=1)                                                                     cggcATGCTGATCCGCGATTACTAGCGATTCCAGCTTCATGTAGGCGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCa                                            <  2:2040187/147‑1 (MQ=0)                                                                         atgcTGATCCGCGATTACTAGCGATTCCAGCTTCATGTAGGCGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAggtc                                        <  1:1518647/147‑3 (MQ=0)                                                                         atgcTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCg                                       >  1:8378354/4‑150 (MQ=0)                                                                               atccgcgATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAggg                                 <  1:7371031/147‑1 (MQ=0)                                                                               atccgcgATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAggg                                 <  1:5638953/147‑1 (MQ=0)                                                                               atccgcgATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAggg                                 <  1:1969153/147‑1 (MQ=0)                                                                                                   ccaaCTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAGGGCCATGAGGACTTGACGTCAt             <  1:8356669/147‑1 (MQ=0)                                                                                                   ccaaCTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAGGGCCATGAGGACTTGACGTCAt             <  2:771242/147‑1 (MQ=0)                                                                                                           atgtaggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCAcc      <  1:21772/143‑1 (MQ=0)                                                                                                           atgtaggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCAcc      <  1:2016106/143‑1 (MQ=0)                                                                                                            tgtaggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCt     >  2:1976148/7‑150 (MQ=0)                                                                                                              taggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc   >  2:6350059/5‑150 (MQ=0)                                                                                                              taggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc   >  2:6350083/5‑150 (MQ=0)                                                                                                              taggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc   >  2:6873893/5‑150 (MQ=0)                                                                                                              taggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc   <  2:8834504/146‑1 (MQ=0)                                                                                                              cactCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAGGGCCATGAGGACTTGACGTCATCCCCACCTtcc  <  2:3005396/147‑1 (MQ=0)                                                                                                              aaggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc   >  2:6873868/5‑150 (MQ=0)                                                                                                               aggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtcc  >  2:2977507/4‑150 (MQ=0)                                                                                                               aggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtcc  >  2:7438776/4‑150 (MQ=0)                                                                                                               aggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtcc  >  2:7438887/4‑150 (MQ=0)                                                                                                               aggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtcc  >  2:2819442/4‑150 (MQ=0)                                                                                                                            |                                                                                                                                          ACTTCTGGTGCAACCCACTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGACATTCTGATTCGCGATTACTAGCGATTCCGACTTCACGCAGTCGAGTTGCAGACTGCGATCCGGACTACGATCGGTTTTGTGAGATTAGCTCCACCTCGCG‑GCTTGGCAACCCTCTGTACCGACCATTGTAGCACGTGTGTAGCCCAGGCCGTAAGGGCCATGATGACTTGACGTCATCCCCACCTTCC  >  NC_002947/5313236‑5313495

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