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

AGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGC  >  NC_002947/3951021‑3951301                                                                                                                                              |                                                                                                                                             aGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGAc                                                                                                                                       <  1:259154/148‑1 (MQ=255)    tCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGGAGCGAAAGGGCTGCGCAGCCGCCCCGGCAGTGGCAGC‑GCAAGCGCAGAAAGAGGGGCGCGCTGCGCGCGCCTTTCGcgccgc                                                                                                                                     >  4:261019/1‑147 (MQ=38)      gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGTGACGCAAgg                                                                                                                                 <  7:352013/149‑1 (MQ=255)             ccatgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGCGCGGCCTGGCGTCGCGAACGGGCTGCGCAGCACCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGCGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTc                                                                                                                           <  5:327972/148‑1 (MQ=38)                                 aGCAATAGCGCACGCCGTGTGGGCGCGGCCTGGGGGCGCTAACGGGCTGCGCCCCCGCCCGGGCAGTGCCAGC‑TCACGCGCACAAAGTGGGCCGCGCGGCGCGCGCCTTTCGCTACGCATGGCCGCTCCCACAAAAGCGGATGCCACag                                                                                                      <  8:315114/149‑1 (MQ=255)                                       atCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGAc                                                                                                                                       >  1:315206/1‑110 (MQ=255)                                       atCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGAc                                                                                                                                       <  2:315206/110‑1 (MQ=255)                                                           ggCCTTGCGTCGCTAACGGGCTGCGCAGCACCCCCGGCAGGGGCAGC‑TCAAGCGCAGAAAGCGGGGCGCGCTGCGCGCTCCTTGCGCCACGAAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                            <  7:126737/149‑1 (MQ=255)                                                            gCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                                                                            >  6:202328/1‑148 (MQ=255)                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATa                                                                          >  4:60320/1‑149 (MQ=255)                                                                 gCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCa                                                                      >  3:139724/1‑149 (MQ=255)                                                                  cGTCCCAAACGGTCTGCCCCGCACCCCCGGCTGTCCCCCA‑TCAAGCGCCAAAACCGGGGAGCGCGCCGCGCTCCTTCCGCGACTAAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                                                                     <  8:195647/149‑1 (MQ=255)                                                                      gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgac                                                                 >  1:190151/1‑149 (MQ=255)                                                                          aaGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg                                                              >  1:252511/1‑149 (MQ=255)                                                                             ggCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg                                                          >  7:273101/1‑149 (MQ=255)                                                                             ggCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg                                                          >  1:272308/1‑149 (MQ=255)                                                                              gCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc                                                         >  3:65695/1‑149 (MQ=255)                                                                                    cagcagCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg                                                   <  3:60320/149‑1 (MQ=255)                                                                                     agcagcCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTggg                                                  <  2:272308/149‑1 (MQ=255)                                                                                      gcagcCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                                 <  6:412863/149‑1 (MQ=255)                                                                                          ccccGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGACAGCGGCAGGGTCaa                                              >  8:374891/1‑148 (MQ=255)                                                                                              ggCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc                                                      <  8:95055/136‑1 (MQ=255)                                                                                              ggCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc                                                      >  7:95055/1‑136 (MQ=255)                                                                                                     gCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                                 >  5:8655/1‑134 (MQ=255)                                                                                                     gCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                                 <  6:8655/134‑1 (MQ=255)                                                                                                        gC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCt                                         <  3:19207/139‑1 (MQ=255)                                                                                                        gC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCt                                         >  4:19207/1‑139 (MQ=255)                                                                                                         c‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTgg                              <  2:252511/149‑1 (MQ=255)                                                                                                            cAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATa                                                                          >  5:240137/1‑103 (MQ=255)                                                                                                            cAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATa                                                                          <  6:240137/103‑1 (MQ=255)                                                                                                             aaGCGCAGAAAGAGTGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCa                                      <  2:134066/138‑1 (MQ=255)                                                                                                             aaGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCa                                      >  1:134066/1‑138 (MQ=255)                                                                                                             aaGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGc                            >  4:250602/1‑148 (MQ=255)                                                                                                              cgcgcAGAAACCGGGGCGCGCGGCGCGCTCCTTGCGCCACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCaa                          <  7:93569/148‑1 (MQ=255)                                                                                                               gcgcAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAAt                         >  1:190822/1‑149 (MQ=255)                                                                                                                  cAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGc                                       >  5:307249/1‑132 (MQ=255)                                                                                                                  cAGAAAGAGGGGAGCGCTGCGCGCTACTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGc                                       <  6:307249/132‑1 (MQ=255)                                                                                                                       cgaaaggAGCGCGCAGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCAGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTg                 >  7:357197/6‑149 (MQ=255)                                                                                                                         agGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGc               >  5:112379/1‑149 (MQ=255)                                                                                                                            ggAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGt            >  6:186030/1‑149 (MQ=255)                                                                                                                            ggAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGt            <  5:186030/149‑1 (MQ=255)                                                                                                                            ggAGCGCGCAGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCAGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGt            <  8:357197/149‑1 (MQ=255)                                                                                                                             gAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTg           <  4:65695/149‑1 (MQ=255)                                                                                                                                     cgcgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCgg    <  4:139724/148‑1 (MQ=255)                                                                                                                                      gcgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGc  <  3:250602/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                             AGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGC‑TCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGC  >  NC_002947/3951021‑3951301

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