Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,313,324	C→G	13.8%	intergenic (‑595/+119)	PP_16SG ← / ← msrQ	16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,313,324	0	C	G	13.8%	75.8 / 6.1	29	intergenic (‑595/+119)	PP_16SG/msrQ	16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit
Reads supporting (aligned to +/- strand):  ref base C (16/9);  new base G (2/2);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

GTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGC  >  NC_002947/5313178‑5313470                                                                                                                                                   |                                                                                                                                                   gTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGccgcgcgg                                                                                                                                                 <  1:655450/150‑2 (MQ=255)  tCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCcgcg                                                                                                                                                >  1:586548/1‑150 (MQ=255)         aaGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                         <  2:586548/150‑1 (MQ=255)           gTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAACGCCGCCCGCGCGGCGCTCg                                                                                                                                       >  2:1383606/1‑150 (MQ=255)              gggCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                    >  2:1353318/1‑150 (MQ=255)                gcgcATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGACAGCGGCGAATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCATCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGAtctc                                                                                                                                  >  1:812033/1‑150 (MQ=14)                        taAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGTCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCATAGgcgc                                                                                                                          <  2:378145/150‑1 (MQ=11)                                      agGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAgg                                                                                                            <  2:227925/150‑1 (MQ=255)                                         gCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                         <  1:1353318/150‑1 (MQ=255)                                         gCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                         <  2:165346/150‑1 (MQ=255)                                                 cTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCACCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCtctc                                                                                                 >  1:1362694/1‑150 (MQ=255)                                                     taTTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGAATGCTCTGACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCcaca                                                                                             <  1:331361/149‑1 (MQ=21)                                                       tttCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAGACAAAGTCGGGCCCTCTCCACAga                                                                                           >  2:331361/1‑149 (MQ=255)                                                           aaGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCAtt                                                                                       >  2:456168/1‑150 (MQ=255)                                                                aGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCAAAAAAAGTTTCAGCGCTTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAt                                                                                  >  2:481313/1‑150 (MQ=255)                                                                   gAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtctc                                                                               <  2:390036/150‑1 (MQ=255)                                                                   gAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtctc                                                                               <  1:997227/150‑1 (MQ=255)                                                                    aaaTATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCa                                                                              >  1:62539/1‑150 (MQ=255)                                                                                                   ccAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCActct                                               >  1:302636/1‑150 (MQ=255)                                                                                                   ccAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCActct                                               >  2:1073736/1‑150 (MQ=255)                                                                                                         gtTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAAcc                                         <  1:456168/150‑1 (MQ=255)                                                                                                                  aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCcgcgc                                                                                                                                               >  1:1246162/1‑39 (MQ=255)                                                                                                                  aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCcgcgc                                                                                                                                               <  2:1246162/39‑1 (MQ=255)                                                                                                                  aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTAccc                                >  1:960241/1‑150 (MQ=255)                                                                                                                               tctgtGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCATAGGCGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCACGAGAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCaa                   >  1:375236/5‑150 (MQ=37)                                                                                                                                  aTGAGTTCGAGCGCCGCCCGCGAGGCGCTCGATCTCACAGAGACAAGAGAAACAAAGACGGGCCGTCTCCACAGCCATTGTTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCCAACCCCCCTACCCTACTACCTGGCAACCt                >  2:923749/1‑150 (MQ=255)                                                                                                                                       aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGa           >  2:828020/1‑150 (MQ=255)                                                                                                                                          gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTgg             >  1:896469/1‑145 (MQ=255)                                                                                                                                          gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTgg             <  2:896469/145‑1 (MQ=255)                                                                                                                                            gcgcCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgcgc      >  2:923570/1‑150 (MQ=255)                                                                                                                                                cgcccgcGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCcgc  >  1:107357/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGC  >  NC_002947/5313178‑5313470

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