Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,951,159	T→C	intergenic (‑123/+72)	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,159	0	T	C	100.0%	83.5 / NA	25	intergenic (‑123/+72)	PP_3486/PP_3487	cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/9);  total (16/9)

GCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGC  >  NC_002947/3951046‑3951301                                                                                                                  |                                                                                                                                                gCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcg                                                                                                                                            <  2:75296/119‑1 (MQ=255)                   aaGCCCTGTGGGAGCGGCCGTTGGGAGGGAAAGGGCGGGGCCGCCGGCCCGGCAGGGGCAGCGGAATCGCAGCAAGAGGGGGGCGCGGCCCCCGGCGTTAGCGGCGCAAGGGCGCGACCACAAAAAAGGAATGCCCCAAAGGCCAAtaa                                                                                            >  8:238209/1‑147 (MQ=255)                                    ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                            <  1:37072/147‑1 (MQ=255)                                    ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                            >  4:38317/1‑147 (MQ=255)                                    ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                            >  2:37072/1‑147 (MQ=255)                                    ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                            <  3:38317/147‑1 (MQ=255)                                               gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCGCCAGGCGGCCCATGATACTCaccgaccg                                                               >  5:3069/1‑149 (MQ=255)                                               gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgaccg                                                               >  6:328197/1‑149 (MQ=255)                                                           cagcagCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg                                                   <  5:328197/149‑1 (MQ=255)                                                              cagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                                 >  2:339192/1‑148 (MQ=255)                                                                 ccccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                              >  7:136606/1‑148 (MQ=255)                                                                      tcAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTAGCGACGCAAGGCCGCTCCCACAAAAACGC‑ATGCCACAGAGCCCAAATCCCCAGGCGGCACATGATACTCACCGCCCGGCCAGCGCCTgg                                                   >  8:136292/2‑138 (MQ=37)                                                                      gCAGTGGCAGCTCAAGCGCAGAAAGAGGTGAGCGCTGCGCTCTCCTTTCGCGACGCAATGCCGCTCCCACAAAAGCGG‑ATGCCAGAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg                                                   <  7:136292/138‑1 (MQ=255)                                                                      gCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTg                                        <  1:339192/149‑1 (MQ=255)                                                                        aGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCa                                      >  1:85861/1‑149 (MQ=255)                                                                          tGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATg                                    >  5:271240/1‑149 (MQ=255)                                                                                               agGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgccgc                   <  8:95146/145‑1 (MQ=255)                                                                                               agGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgccgc                   >  7:95146/1‑145 (MQ=255)                                                                                                      cgcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg                                                   <  6:75262/106‑1 (MQ=255)                                                                                                      cgcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg                                                   >  5:75262/1‑106 (MQ=255)                                                                                                         tGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATc                        >  4:280129/1‑130 (MQ=255)                                                                                                         tGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATc                        <  3:280129/130‑1 (MQ=255)                                                                                                         tGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCg     >  5:57216/1‑149 (MQ=255)                                                                                                           cgcgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCggg   >  3:54132/1‑149 (MQ=255)                                                                                                            gcgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTACATCCGGGc  >  8:200499/1‑149 (MQ=255)                                                                                                                  |                                                                                                                                                GCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGG‑ATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGC  >  NC_002947/3951046‑3951301

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: ‑