Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,951,161	A→T	intergenic (‑125/+70)	PP_3486 ← / ← PP_3487	cytochrome c/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,951,161	0	A	T	100.0%	83.5 / NA	26	intergenic (‑125/+70)	PP_3486/PP_3487	cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/8);  total (18/8)

GCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGC  >  NC_002947/3951032‑3951301                                                                                                                                  |                                                                                                                                              gccatgccaAGCACGCAGGTAAGCCATATCGCAAGCCGTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAgg                                                                                                                                  <  4:166552/143‑1 (MQ=255)                                           gcaGCGGCCTTGCGTCGCTAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCATAAAGAGGGGCGCGCTGCGCGCTCCTTTCGCTACGAAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGcc                                                                                 <  1:39012/147‑1 (MQ=38)                                             aGCGGCCTTGCGTCGCGAAATGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGATAGAGTGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TACCCATGCGGCCCa                                                                               >  2:32839/1‑149 (MQ=38)                                                  ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATa                                                                          >  2:11076/1‑149 (MQ=255)                                                      gCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCa                                                                      >  5:379812/1‑149 (MQ=255)                                                             gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCaccgaccg                                                               >  5:289621/1‑149 (MQ=255)                                                             gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCaccgaccg                                                               >  3:167586/1‑149 (MQ=255)                                                                ggggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGcc                                                            >  6:36613/2‑149 (MQ=39)                                                                       cgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTc                                                                       >  6:148138/1‑131 (MQ=255)                                                                       cgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTc                                                                       <  5:148138/131‑1 (MQ=255)                                                                               ccccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                              >  5:294450/1‑148 (MQ=255)                                                                               ccccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAAc                                             <  4:167586/149‑1 (MQ=255)                                                                                cccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                              >  8:62238/1‑147 (MQ=255)                                                                                 ccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAg                                           <  2:59205/149‑1 (MQ=255)                                                                                 ccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAg                                           >  1:59205/1‑149 (MQ=255)                                                                                  cGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGc                                          >  8:8568/1‑149 (MQ=255)                                                                                     cAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGc                                       <  7:8568/149‑1 (MQ=255)                                                                                           cAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcgc                                 >  6:250414/1‑149 (MQ=255)                                                                                                   gcgcAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCGCCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAAt                         >  4:252050/1‑149 (MQ=255)                                                                                                       aGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgcc                     >  4:235792/1‑149 (MQ=255)                                                                                                       aGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgcc                     >  7:175067/1‑149 (MQ=255)                                                                                                         aaaGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgccgc                   >  4:355287/1‑149 (MQ=255)                                                                                                           agagGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTg                 <  3:355287/149‑1 (MQ=255)                                                                                                               gggAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACGGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCgg             >  4:302829/1‑149 (MQ=255)                                                                                                                       tGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATAACCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGATGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATcc      >  1:114601/1‑149 (MQ=255)                                                                                                                          gcgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGc  <  1:11076/149‑1 (MQ=255)                                                                                                                                  |                                                                                                                                              GCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAAT‑TCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGC  >  NC_002947/3951032‑3951301

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