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

CTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCG  >  NC_002947/3951018‑3951289                                                                                                                                                |                                                                                                                                  accccctCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc‑gc                                                                                                                              <  1:156694/142‑1 (MQ=255)      gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc‑gt                                                                                                                              >  2:156694/1‑142 (MQ=255)          aaGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGc                                                                                                                    <  6:339258/149‑1 (MQ=255)           tgTTGCAATGCAACCCACCCTGTTACGCCATACCCCACGCCGTGTGGGCGCGGCCTTGCGTCGCTAATGGGCTGCGCAGCCGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGcc                                                                                                                   <  6:66745/148‑1 (MQ=21)            gTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGAGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc‑gc                                                                                                                              <  8:216458/137‑1 (MQ=255)            gTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGAGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc‑gc                                                                                                                              >  7:216458/1‑137 (MQ=255)                           aCGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCATGCGTCGCGAAACGGCAGCGCAGCAGCCCCCGCAGTGGCAGCACAAGCGCAGAAAAAAGGGAGCGCTGCGCGCTCCATTA‑GCGACGCAAGGCCGCTCCCACaaaaacc                                                                                                    >  6:201522/1‑145 (MQ=37)                             gCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGACTCGCGCAAGCGCTGGGCCGCACCCCCCGCACGGGCACCACCACCGGCCAAAAAGGGGGGCGCAGGGCCCTCCATTA‑GGGCCGCAAGGCCCCGCCCCCAAAAGCgg                                                                                                   >  2:135662/1‑147 (MQ=255)                                      ccATATCGCAAGCCCTGTGGGAGCGGCCATTCGGAGCGAAAGGGCATGGCCCCAGCCCCGGCAGTGGCCGCGCCAGCCCACAAACAGGGGCGCCCTGCGCGGGCCGTTA‑GCGGCGCACGGCCGCGCCCACCAAACCGGCTGGCCCagac                                                                                        >  6:192600/1‑148 (MQ=255)                                                             cGGCGTTGCGTCGCTAAAGGGCTGCGCAGCAGCCCCGGCCGTGGCAGCTCAAGCGCAGAAAGAGGGGCGCGCTGCGCGCTCCTTTC‑GCGACGAAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATg                                                                 <  5:201522/149‑1 (MQ=38)                                                               gCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                >  5:310686/1‑148 (MQ=255)                                                                       tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGc                                                                                      >  6:88542/1‑118 (MQ=255)                                                                       tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGc                                                                                      <  5:88542/118‑1 (MQ=255)                                                                       tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCacc                                                        >  8:74757/1‑148 (MQ=255)                                                                         gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgac                                                     >  6:180696/1‑149 (MQ=255)                                                                          cgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAgg                                                                                       >  4:259522/1‑114 (MQ=255)                                                                          cgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAgg                                                                                       <  3:259522/114‑1 (MQ=255)                                                                           gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgaccg                                                   >  5:380952/1‑149 (MQ=255)                                                                            aaaGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg                                                  <  6:310686/149‑1 (MQ=255)                                                                             aaGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGc                                                 <  6:380952/149‑1 (MQ=255)                                                                               gggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGc                                                 >  8:4373/1‑147 (MQ=39)                                                                                  cTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCg                                            >  3:260677/1‑149 (MQ=255)                                                                                   tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCgg                                           >  5:65595/1‑149 (MQ=255)                                                                                   tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCgg                                           >  7:246217/1‑149 (MQ=255)                                                                                       cagcagCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGCGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg                                       >  2:127452/1‑149 (MQ=255)                                                                                         gcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                     <  4:260677/149‑1 (MQ=255)                                                                                            gcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGccc                                                                    >  8:56382/1‑115 (MQ=255)                                                                                            gcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGccc                                                                    <  7:56382/115‑1 (MQ=255)                                                                                            gcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                  >  3:399315/1‑149 (MQ=255)                                                                                                cGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGc                              >  6:117469/1‑149 (MQ=255)                                                                                                 ggCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGc                                                                                      >  1:151489/1‑92 (MQ=255)                                                                                                 ggCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGc                                                                                      <  2:151489/92‑1 (MQ=255)                                                                                                   cAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTg                            >  8:129742/1‑148 (MQ=255)                                                                                                     gTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCAt                         <  5:117469/149‑1 (MQ=255)                                                                                                        gCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcg                      >  4:181014/1‑149 (MQ=255)                                                                                                          aGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGCCAACAGCTGCATGCGCt                    >  7:343193/1‑149 (MQ=255)                                                                                                               aaGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGc                <  1:127452/148‑1 (MQ=255)                                                                                                                        aaGAGGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGCCAACAGCTGCATGCGCTGGGCAATCGCCGCt      <  8:343193/149‑1 (MQ=255)                                                                                                                         agagGGGAGCGCTGCGCGCTCCTTTCGGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGCGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCt      >  1:228808/1‑149 (MQ=255)                                                                                                                         agagGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCTACAGCTGCATGCGCTGGGCAATCGCCGCTg     >  6:395604/1‑149 (MQ=255)                                                                                                                          gagGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                <  6:352100/89‑1 (MQ=255)                                                                                                                          gagGGGAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                >  5:352100/1‑89 (MQ=255)                                                                                                                            ggggAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCg  >  7:42211/1‑149 (MQ=255)                                                                                                                            ggggAGCGCTGCGCGCTCCTTTC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCg  <  8:246217/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                  CTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATC‑GCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCG  >  NC_002947/3951018‑3951289

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