Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	2,655,897	C→T	P98P (CCG→CCA)	cysH ←	phosphoadenosine phosphosulfate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,655,897	0	C	T	100.0%	80.1 / NA	25	P98P (CCG→CCA)	cysH	phosphoadenosine phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/14);  total (11/14)

CGGGCTTTGGTCGCGGCGCTGGCCGGTGGCCCAGGCGCTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAACGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCAGGCGCCCAGTGTCGAGGCTGAATACCTTTACCTGCTTGTTCA  >  NC_002947/2655756‑2656033                                                                                                                                              |                                                                                                                                         cGGGCTTTGGTCGCGGCGCTGGCCGGTGGCCCAGGCGCTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGt                                                                                                                                      >  1:109705/1‑146 (MQ=255) cGGGCTTTGGTCGCGGCGCTGGCCGGTGGCCCAGGCGCTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGACGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGt                                                                                                                                      >  6:176262/1‑146 (MQ=255)         gggcgcgGCGCTGGGCGGTGGCCCAGGCGCTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTgg                                                                                                                             <  5:38606/144‑1 (MQ=255)                 cgcTGGCCGGTGGCCCAGGCGCTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCgg                                                                                                                    <  2:114760/148‑1 (MQ=255)                        cGGTGGCCCAGGCGCTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAg                                                                                                             <  8:41739/148‑1 (MQ=255)                               ccAGGCGCTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCg                                                                                                      >  1:155234/1‑148 (MQ=255)                               ccAGGCGCTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCg                                                                                                      >  5:149459/1‑148 (MQ=255)                                   gcgcTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATgg                                                                                                  >  1:18058/1‑148 (MQ=255)                                     gcTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATgg                                                                                                  >  1:103401/1‑146 (MQ=255)                                     gcTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGc                                                                                                <  4:92416/148‑1 (MQ=255)                                                  aGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGtt                                                                                   <  8:171026/148‑1 (MQ=255)                                                                       aTCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAAc                                                               <  5:176262/147‑1 (MQ=255)                                                                          ttGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCg                                                             >  8:110301/1‑146 (MQ=255)                                                                          tggCGAGTGCTGCGGCATTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCg                                                             >  2:14628/3‑146 (MQ=255)                                                                             cGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTTAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAg                                                         <  3:47888/147‑1 (MQ=255)                                                                             cGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGt                                                           <  2:66910/145‑1 (MQ=255)                                                                              gAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGt                                                       >  6:89998/1‑148 (MQ=255)                                                                                         caCTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCa                                            <  1:112527/148‑1 (MQ=255)                                                                                         caCTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCa                                            <  1:109481/148‑1 (MQ=255)                                                                                                 aTGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCAGGCGCCCa                                    <  7:110301/148‑1 (MQ=255)                                                                                                       gTCCTTGTAGAAGCTGAACAGGCCCGTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTCCTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCAGGCGCCCAgtgt                                >  7:35629/1‑146 (MQ=255)                                                                                                          cTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCAGGCGCCCAGTGTCGAgg                           >  4:90588/1‑148 (MQ=255)                                                                                                               aGAACCTGACCAGGCCCTTTTCCTTGACGACTGGGGCGCGCTGGGCGCGGGCGGGGCGGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCAGGCGCCCAGTGTCGAGGCTGaa                      <  2:72139/148‑1 (MQ=255)                                                                                                                                   tCCTTGACGAATGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCAGGCGCCCAGTGTCGAGGCTGAATACCTTTACCTGCTTGtt    <  6:56347/146‑1 (MQ=255)                                                                                                                                    ccTTGACGAATGGGGCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCAGGCGCCCAGTGTCGAGGCTGAATACCTTTACCTGCTTGTTCa  <  5:120269/147‑1 (MQ=255)                                                                                                                                              |                                                                                                                                         CGGGCTTTGGTCGCGGCGCTGGCCGGTGGCCCAGGCGCTGACGGTGGCCAGCTTGCGTCGCAGCGGCTCGATCTTGCGAATGCCGCAGCACTCGCCATGGCCGTCCTTGTAGAAGCTGAACAGGCCCTTTTCCTTGACGAACGGGTCGAGCTTGGCGCGGTCGGGGCTGAGGATCTCGATGGGCAGGTTGTACTGTTCGCGTACCTGGTCGATGAACCGGTAGGTCTCCGGGTGCAGGCGCCCAGTGTCGAGGCTGAATACCTTTACCTGCTTGTTCA  >  NC_002947/2655756‑2656033

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