Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,951,161	A→T	100%	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%	93.2 / NA	29	intergenic (‑125/+70)	PP_3486/PP_3487	cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (20/9);  total (20/9)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑TTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGC  >  NC_002947/3951026‑3951301                                                                                                                                          |                                                                                                                                             gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAgg                                                                                                                                 >  1:205268/1‑149 (MQ=255)  aaGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGc                                                                                                                                <  2:205268/149‑1 (MQ=255)                                 tatCGCACGCCCGGTGGGAGGGGCCTTGCGTCGCTAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTcg                                                                                                                                           <  6:25989/107‑1 (MQ=37)                                 tatCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTcg                                                                                                                                           >  5:25989/1‑107 (MQ=255)                                                      ggCCTTGCGTCCCCAAAGGCCTGCGATCCAGCCCCGGCAGTGGCAGCTCAAGCCCACAAAGAGGGTCGCGGTGCGCGCGC‑‑CTTGCGCCACCCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                            <  7:168283/149‑1 (MQ=255)                                                                     aaGGGCTGCGCAGCAGCCCCGGCAGTGGCGGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGc                                                             >  3:218630/1‑149 (MQ=255)                                                                        ggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg                                                          >  8:250335/1‑149 (MQ=255)                                                                          cTGCGCAGCAGCCCCGTCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGAGCTCCGTTTTCCCGACGCAAGGCCGCACCCACAAAACCGGATTCCACAGAGGCCAATTCCACAGGCGGCCCATGATACTCACCGACCGGCCAg                                                          >  1:290689/1‑149 (MQ=255)                                                                            gcgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                                                                     <  6:316174/134‑1 (MQ=255)                                                                            gcgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                                                                     >  5:316174/1‑134 (MQ=255)                                                                                  cagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                            <  1:130516/121‑1 (MQ=255)                                                                                  cagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                            >  2:130516/1‑121 (MQ=255)                                                                                    gcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGGCaa                                              >  6:183838/1‑149 (MQ=255)                                                                                          gCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTg                                        >  4:360614/1‑149 (MQ=255)                                                                                              tGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATg                                    >  7:241989/1‑149 (MQ=255)                                                                                                gCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcg                                  >  5:186571/1‑149 (MQ=255)                                                                                                  aGCTCAAGTGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCt                                >  4:258954/1‑149 (MQ=255)                                                                                                    cTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcgc                                 <  6:59717/146‑1 (MQ=255)                                                                                                    cTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcgc                                 >  5:59717/1‑146 (MQ=255)                                                                                                      cAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTggg                             >  8:3756/1‑148 (MQ=255)                                                                                                       aaGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTggg                             <  4:218630/147‑1 (MQ=255)                                                                                                       aaGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGc                            >  5:290797/1‑148 (MQ=255)                                                                                                          cgcAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATc                        >  3:103862/1‑149 (MQ=255)                                                                                                           gcAGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATTCTCACAGACCGGCCAGCGGCTGGTTCAACAGCTGCATGCGCTGGGCAATcg                       >  6:106584/1‑149 (MQ=39)                                                                                                             aGAAAGAGGGGAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgcc                     >  6:130906/1‑149 (MQ=255)                                                                                                                     gggAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCgg             <  5:35371/149‑1 (MQ=255)                                                                                                                     gggAGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCgg             >  6:35371/1‑149 (MQ=255)                                                                                                                        aGCGCTGCGCGCTC‑‑CTTTCGCGACGCAAGGCCGCTCACACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTg           >  1:247228/1‑148 (MQ=255)                                                                                                                                gcgcTC‑‑CTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGc  <  3:360614/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                             GAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTC‑‑TTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGC  >  NC_002947/3951026‑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: ‑