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%	72.3 / NA	23	intergenic (‑125/+70)	PP_3486/PP_3487	cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (14/9);  total (14/9)

GGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGG  >  NC_002947/3951016‑3951299                                                                                                                                                  |                                                                                                                                           ggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcg                                                                                                                                         <  3:246269/149‑1 (MQ=255)      aGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACg                                                                                                                                    <  6:340612/149‑1 (MQ=255)                tgccatgccaAGCACGCAGGTAAGGCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGc                                                                                                                           <  7:178748/148‑1 (MQ=255)                        caAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGc                                                                                                                              <  2:90932/137‑1 (MQ=255)                        caAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCACGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGc                                                                                                                              >  1:90932/1‑137 (MQ=255)                        caAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGAAAGGGCTGCGCAGCCGCCCCCGCAGGGGCAGCGGAAACGGAGAAAGAGGGGAGCGCGGCGCGCGCCTTTAGCGCCGCAAGGCCGCTCCCACaa                                                                                                                   >  8:258008/1‑148 (MQ=255)                                                                  ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCCATTCCCCAGGCGGCCCATGATa                                                                        >  3:245057/1‑149 (MQ=255)                                                                            cgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgacc                                                              >  1:47438/1‑149 (MQ=255)                                                                             gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgaccg                                                             >  4:138644/1‑149 (MQ=255)                                                                              ccggggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg                                                            >  4:83959/4‑149 (MQ=39)                                                                                   gCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGAGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCAAAAACGCGGATGCCACAGAGGCCAATTCCACAGGCGGCCCCTGATACTCACCGACCGGCCAGc                                                       >  6:360955/1‑149 (MQ=255)                                                                                    cTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCg                                                      >  6:184472/1‑149 (MQ=255)                                                                                                   ggCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCt                                       <  5:184472/149‑1 (MQ=255)                                                                                                    gCAGTGGCACCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTg                                      <  7:258008/149‑1 (MQ=255)                                                                                                          gCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcg                                >  5:25662/1‑149 (MQ=255)                                                                                                           cAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcgc                               >  6:368684/1‑149 (MQ=255)                                                                                                                     gcAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcg                     >  4:338974/1‑149 (MQ=255)                                                                                                                             aaaggAGCGCGCAGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGc             >  3:88375/4‑149 (MQ=255)                                                                                                                              ggggAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCg            <  2:47438/149‑1 (MQ=255)                                                                                                                               gggAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCgg           >  2:334869/1‑149 (MQ=255)                                                                                                                               gggAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCgg           <  1:334869/149‑1 (MQ=255)                                                                                                                                       tGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCg   >  5:11327/1‑149 (MQ=255)                                                                                                                                        cagcgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCgg  <  4:88375/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                           GGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGG  >  NC_002947/3951016‑3951299

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