Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,837,281	A→G	11.3%	S430G (AGC→GGC)	ccoN‑I →	cbb3‑type cytochrome c oxidase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,837,281	0	A	G	11.3%	155.9 / 10.2	62	S430G (AGC→GGC)	ccoN‑I	cbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base G (3/4);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

ATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAG  >  NC_002947/4837134‑4837422                                                                                                                                                     |                                                                                                                                              aTCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACCTCGTTCGTGGAAACCTTGGTGgcca                                                                                                                                               <  1:3724/149‑1 (MQ=255) aTCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGgccag                                                                                                                                              >  8:328237/1‑149 (MQ=255)    aaTGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAGCCTGGGTGgctggcca                                                                                                                                           >  7:231265/1‑149 (MQ=255)        cgcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCa                                                                                                                                       <  2:21528/149‑1 (MQ=255)            cTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCt                                                                                                                                   <  1:75210/149‑1 (MQ=255)                  gctggcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAGCCTGGGTGGCTGGCCACCAAGGTTtcca                                                                                                                               <  8:231265/147‑3 (MQ=18)                  gctggcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgtc                                                                                                                              <  7:105723/148‑1 (MQ=255)                        caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCt                                                                                                                       >  6:291868/1‑149 (MQ=255)                                     gTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                          <  7:328237/149‑1 (MQ=255)                                     gTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                          >  5:63929/1‑149 (MQ=255)                                               tCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCTGGCCAccaa                                                                                                                                       <  4:341961/110‑3 (MQ=255)                                               tCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCTGGCCAccaa                                                                                                                                       >  3:341961/1‑108 (MQ=255)                                                  ccTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCATGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                             >  4:218027/1‑149 (MQ=255)                                                    tCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                   >  3:69784/1‑141 (MQ=255)                                                    tCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                   <  4:69784/141‑1 (MQ=255)                                                      gATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCa                                                                                         <  5:291868/149‑1 (MQ=255)                                                       aTGTTGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCAt                                                                                        <  2:124172/149‑1 (MQ=255)                                                           gggTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTc                                                                                    <  3:218027/149‑1 (MQ=255)                                                             gTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                   >  1:30697/1‑132 (MQ=255)                                                             gTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                   <  2:30697/132‑1 (MQ=255)                                                              tCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTccag                                                                                 >  5:121234/1‑147 (MQ=255)                                                                      aTCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGACTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTg                                                                         <  3:258712/149‑1 (MQ=255)                                                                         gCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                   >  8:166020/1‑120 (MQ=255)                                                                          ccgaGGGCCTGATGTGGGGCGCGGTCTTCTGCGACGGCACGCTCTCGTA‑CTCGTTCGTGGATACCTTGGTGGCCAGCGTCCGAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGGATGTTCCTGAt                                                                               <  5:221092/136‑1 (MQ=255)                                                                          cccACGGCCTGATGTGGCGCGCGCACAACAGCGACGGCACGCTCACATA‑CTCGAACGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTACGTCCGCGGCGCCATCTTCCTCACCCCCATGTTCCTGAt                                                                               >  6:221092/1‑139 (MQ=255)                                                                                 ccTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                              >  1:80576/1‑149 (MQ=255)                                                                                    gATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCa                                                           >  3:261538/1‑149 (MQ=255)                                                                                     aTGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCAc                                                          >  1:105023/1‑149 (MQ=255)                                                                                     aTGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCAc                                                          >  1:101043/1‑149 (MQ=255)                                                                                       gtgGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCg                                                                                            <  1:152072/113‑1 (MQ=255)                                                                                       gtgGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCg                                                                                            >  2:152072/1‑113 (MQ=255)                                                                                       gtgGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCg                                                        >  4:304650/1‑149 (MQ=255)                                                                                        tgGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCg                                                        <  4:59781/148‑1 (MQ=255)                                                                                        tgGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGt                                                       <  2:101043/149‑1 (MQ=255)                                                                                        tgGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGt                                                       <  2:105023/149‑1 (MQ=255)                                                                                          gcgcgcGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                             <  3:329119/109‑1 (MQ=255)                                                                                          gcgcgcGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                             >  4:329119/1‑109 (MQ=255)                                                                                            gcgcgGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGt                                                   >  3:176462/1‑149 (MQ=255)                                                                                                ggTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGc                                               >  8:138549/1‑149 (MQ=255)                                                                                                  tCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCa                                             <  6:63929/149‑1 (MQ=255)                                                                                                       aGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCt                                        >  7:307944/1‑149 (MQ=255)                                                                                                        gCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTg                                                      <  2:52437/134‑1 (MQ=255)                                                                                                        gCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTg                                                      >  1:52437/1‑134 (MQ=255)                                                                                                             ggCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGt                                  >  6:227533/1‑149 (MQ=255)                                                                                                                aCGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGc                               <  5:227533/149‑1 (MQ=255)                                                                                                                    tCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTg                                      >  3:16768/1‑138 (MQ=255)                                                                                                                    tCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTg                                      <  4:16768/138‑1 (MQ=255)                                                                                                                      aCCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGc                                                                                          >  4:269887/1‑84 (MQ=255)                                                                                                                       ccTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGc                         >  3:60224/1‑148 (MQ=255)                                                                                                                         tA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCATCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGaa                      >  3:234986/1‑149 (MQ=255)                                                                                                                               gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                  <  4:176462/148‑1 (MQ=255)                                                                                                                                ttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTAGTCGTCCGCTTCGTCGGCGGCGCGAAATTCCACACCGGCATGTTCCTGATGGcc                                                                           >  7:81781/1‑90 (MQ=255)                                                                                                                                   gtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTg                          >  3:114678/1‑136 (MQ=255)                                                                                                                                   gtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTg                          <  4:114678/136‑1 (MQ=255)                                                                                                                                   gtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt             >  3:170871/1‑149 (MQ=255)                                                                                                                                        aaCCTTGGTGGCCAGCCAGCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt        <  4:234986/149‑1 (MQ=255)                                                                                                                                        aaCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt        >  1:97408/1‑149 (MQ=255)                                                                                                                                         gccTGGGTGGCTGGCCACCAAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTg                                      <  2:139927/117‑1 (MQ=255)                                                                                                                                         aCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt        <  8:4459/148‑1 (MQ=255)                                                                                                                                         aCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt        >  7:4459/1‑148 (MQ=255)                                                                                                                                           ctgggtggCTGGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag     >  8:108695/4‑149 (MQ=255)                                                                                                                                            ttGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag     >  5:206027/1‑148 (MQ=255)                                                                                                                                            tgggtggCTGGCCACCAAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg           >  7:228707/3‑142 (MQ=255)                                                                                                                                            tgggtggCTGGCCACCAAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg           <  8:228707/140‑1 (MQ=255)                                                                                                                                               gtggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgaggag  <  6:206027/148‑1 (MQ=255)                                                                                                                                                  gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGcccc                            >  7:313790/1‑119 (MQ=255)                                                                                                                                                  gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGcccc                            <  8:313790/119‑1 (MQ=255)                                                                                                                                                    cagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTg       <  6:337430/138‑1 (MQ=255)                                                                                                                                                    cagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTg       >  5:337430/1‑138 (MQ=255)                                                                                                                                                     |                                                                                                                                              ATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTA‑CTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAG  >  NC_002947/4837134‑4837422

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