Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,312,050	A→T	50.5%	noncoding (510/2902 nt)	PP_23SG ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,312,050	0	A	T	50.5%	‑5.2 / 18.1	14	noncoding (510/2902 nt)	PP_23SG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base T (3/4);  total (5/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 = 1.02e-01

CTATTAAGACTCGCTTTCGCTACGCCT‑CCCCTATTCGGTTAAGCTCGCCACTGAATATAAGTCGCTGACCCATTATACAAAAGGTACGCAGTCACC‑‑TAACAAA‑‑GTAGGCTCCCACTGCTTGTACGCATACGGTTTCAGGTTCTATTTCACTCCCCTCTCCGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGTCAGTAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTC  >  NC_002947/5311922‑5312169                                                                                                                                      |                                                                                                                        ctgttcaGACTCGCTTTCGCTGCGCCTACACCTACCGGCTTAAGCTTGCTCGTCAGACTAAGTCGCTGACCCATTATACAAAAGGTACGCTGTCACCCAGAACGAATCTTGGGCTCCAACTGTTTGTAGGCATTCGGTTTCAGGTACTAt                                                                                                         <  1:1637909/144‑1 (MQ=0) ctgttcaGACTCGCTTTCGCTGCGCCTACACCTACCGGCTTAAGCTTGCTCGTCAGACTAAGTCGCTGACCCATTATACAAAAGGTACGCTGTCACCCAGAACGAATCTTGGGCTCCAACTGTTTGTAGGCATTCGGTTTCAGGTACTAt                                                                                                         >  2:6486792/7‑150 (MQ=0)  tgttcaGACTCGCTTTCGCTGCGCCTACACCTACCGGCTTAAGCTTGCTCGTCAGACTAAGTCGCTGACCCATTATACAAAAGGTACGCTGTCACCCAGAACGAATCTTGGGCTCCAACTGTTTGTAGGCATTCGGTTTCAGGTACTAtt                                                                                                        >  2:3269537/6‑150 (MQ=0)  tgttcaGACTCGCTTTCGCTGCGCCTACACCTACCGGCTTAAGCTTGCTCGTCAGACTAAGTCGCTGACCCATTATACAAAAGGTACGCTGTCACCCAGAACGAATCTTGGGCTCCAACTGTTTGTAGGCATTCGGTTTCAGGTACTAtt                                                                                                        <  2:8742954/145‑1 (MQ=0)   gttcaGACTCGCTTTCGCTGCGCCTACACCTACCGGCTTAAGCTTGCTCGTCAGACTAAGTCGCTGACCCATTATACAAAAGGTACGCTGTCACCCAGAACGAATCTTGGGCTCCAACTGTTTGTAGGCATTCGGTTTCAGGTACTAttt                                                                                                       >  1:4982995/5‑150 (MQ=0)   gttcaGACTCGCTTTCGCTGCGCCTACACCTACCGGCTTAAGCTTGCTCGTCAGACTAAGTCGCTGACCCATTATACAAAAGGTACGCTGTCACCCAGAACGAATCTTGGGCTCCAACTGTTTGTAGGCATTCGGTTTCAGGTACTAttt                                                                                                       <  1:7640622/146‑1 (MQ=0)    ttcaGACTCGCTTTCGCTGCGCCTACACCTACCGGCTTAAGCTTGCTCGTCAGACTAAGTCGCTGACCCATTATACAAAAGGTACGCTGTCACCCAGAACGAATCTTGGGCTCCAACTGTTTGTAGGCATTCGGTTTCAGGTACTATTTc                                                                                                      <  1:8289862/147‑1 (MQ=0)    ttcaGACTCGCTTTCGCTGCGCCTACACCTACCGGCATAAGCTTGCTCGTCAGACTAAGTCGCTGACCCATTATACAAAAGGTACGCTGTCACCCACAACGAATCTTGGGCTCCCACTGTTTGTAGGCATACGGTTTCAGGTACTAttta                                                                                                      >  1:6874133/4‑149 (MQ=0)                    ctaagacT‑CCCATATTCGGTTAAGCTCGCCACTGAATATAAGTCGCTGACCCATTATACAAAAGGTACGCAGTCTCC‑‑TAACAAA‑‑GTAGGCTCCCACTGCTTGTACGCATACGGTTTCAGGTTCTATTTCACTCCCCTCTCCGGGTTTCtt                                                                                 <  2:1255648/144‑1 (MQ=1)                             ccacTATTCGGTTAAGCTCGCCACTGAATATAAGTCGCTGACCCATTATACAAAAGGTACGCAGTCACC‑‑TAACAAA‑‑GTAGGCTCCCACTGCTTGTACTCATACGGTTTCAGGTTCTATTTCACTCCCCTCTCCGGGGTTCTTTTCGCCtt                                                                         <  2:8272597/147‑1 (MQ=1)                                      tgtaaatcTCGCAACTGAATATAAGTCGCTTACCCATTATACAAAAGGTACGCAGAAACC‑‑TAAAAAA‑‑GTAGGCTCCAACTTCTAGTACTCATACGGTTTAAGGTTCTATTACACTCCCCTCTCCGGGGTTCTTTTCGACTTTCACTCACg                                                                >  1:6202015/8‑150 (MQ=1)                                                                                 aacggTACGCAGTCACC‑‑TAACAAA‑‑GTAGGCTCCCACTGCTTGTACGCATACGGTTTCAGGTTCTATTTCACTCCCCTCTCCGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGTCAGTAGTATTTAGCCTTgga                     <  1:7232030/147‑1 (MQ=1)                                                                                         ccactCACC‑‑TAACAAA‑‑GTAGGCTCCCACTGCTTGTACGCATACGGTTTCAGGTTCTATTTCACTCCCCTCTCCGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGTCAGTAGTATTTAGCCATGGAGGATGGTc             <  2:1978787/146‑1 (MQ=1)                                                                                                      acgaa‑‑GTAGGCTCCCACTGCTTGTACGCATACGGTTTCAGGTTCTATTTCACTCCCCTCTCCGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGTCAGTAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTc  <  1:2506122/147‑1 (MQ=1)                                                                                                                                      |                                                                                                                        CTATTAAGACTCGCTTTCGCTACGCCT‑CCCCTATTCGGTTAAGCTCGCCACTGAATATAAGTCGCTGACCCATTATACAAAAGGTACGCAGTCACC‑‑TAACAAA‑‑GTAGGCTCCCACTGCTTGTACGCATACGGTTTCAGGTTCTATTTCACTCCCCTCTCCGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGTCAGTAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTC  >  NC_002947/5311922‑5312169

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