Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,042,506	A→C	8.5%	D148A (GAC→GCC)	PP_5295 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,042,506	0	A	C	8.5%	166.0 / 7.8	59	D148A (GAC→GCC)	PP_5295	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base C (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

AGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCAACGAACGCCTGGACAACATGGCCAGCGAAGTC  >  NC_002947/6042360‑6042651                                                                                                                                                   |                                                                                                                                                      aGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGaca                                                                                                                                                     >  2:162430/1‑149 (MQ=255)  gACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGacac                                                                                                                                                    <  3:212180/149‑1 (MQ=255)        aaTGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGTCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGTCGGCCACCAGcac                                                                                                                                              <  8:18278/149‑3 (MQ=255)        aaTGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGTCGGCCACCAGcac                                                                                                                                              <  3:171376/149‑3 (MQ=255)        aaTGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGa                                                                                                                                              >  2:6799/1‑149 (MQ=255)          tGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGAcc                                                                                                                                            >  7:225366/1‑149 (MQ=255)           ggCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGAccc                                                                                                                                           >  3:30978/1‑149 (MQ=255)            gCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGAcccc                                                                                                                                          >  8:237188/1‑149 (MQ=255)               aaCCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGCCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAg                                                                                                                                       <  7:114968/149‑1 (MQ=255)                     tGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACt                                                                                                                                  >  7:167278/1‑148 (MQ=255)                              aCCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGaa                                                                                                                       <  8:295934/149‑1 (MQ=255)                                cATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGc                                                                                                                     >  6:202556/1‑149 (MQ=255)                                         gccgcgcAACTGGTTGCCGCAGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCACTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCt                                                                                                            >  2:272816/1‑149 (MQ=255)                                             cgcAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTc                                                                                                        <  2:34007/149‑1 (MQ=255)                                              gcAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGcc                                                                                                                    >  1:283119/1‑136 (MQ=255)                                              gcAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGcc                                                                                                                    <  2:283119/136‑1 (MQ=255)                                                 aCTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTcgac                                                                                                     >  3:47510/1‑148 (MQ=255)                                                            ccgccgAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGACACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑Gc                                                                                       <  1:57693/149‑1 (MQ=255)                                                                 gAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACt                                                                                 >  1:60310/1‑149 (MQ=255)                                                                  aGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑Tg‑g                                                                                        <  5:73636/141‑1 (MQ=255)                                                                  aGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑Tg‑g                                                                                        >  6:73636/1‑141 (MQ=255)                                                                    gAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTc                                                                                >  7:297303/1‑147 (MQ=255)                                                                    gAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCaa                                                                              <  8:167278/149‑1 (MQ=255)                                                                            tggcCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGTCGGCCACCAGcac                                                                                                                                              >  1:66194/1‑79 (MQ=255)                                                                            tggcCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGTCGGCCACCAGcac                                                                                                                                              <  2:66194/81‑3 (MQ=255)                                                                              gcCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCTGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGgcgc                                                                    <  7:163392/149‑1 (MQ=255)                                                                                   aTCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAgc                                                                <  1:6799/148‑1 (MQ=255)                                                                                   aTCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAgc                                                                <  4:30978/148‑1 (MQ=255)                                                                                     cAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAgcgc                                                              >  4:216000/1‑148 (MQ=255)                                                                                      aaGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGcgcg                                                             >  8:199339/1‑148 (MQ=255)                                                                                      aaGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGcgcg                                                             >  5:67608/1‑148 (MQ=255)                                                                                      aaGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGcgcg                                                             >  2:306706/1‑148 (MQ=255)                                                                                       agcagcCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGcgcg                                                             >  1:161192/1‑147 (MQ=255)                                                                                        gcagcCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAg                                                          <  5:114841/149‑1 (MQ=255)                                                                                         cagcCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAgg                                                         <  2:166773/149‑1 (MQ=255)                                                                                            ccGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGc                                                        >  3:161931/1‑147 (MQ=255)                                                                                            ccGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGc                                                        >  6:230506/1‑147 (MQ=255)                                                                                             cGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAg                                                     <  5:202556/149‑1 (MQ=255)                                                                                                   cccTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTg                                               <  4:161931/149‑1 (MQ=255)                                                                                                    ccTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTgg                                              <  1:162430/149‑1 (MQ=255)                                                                                                    ccTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTgg                                              >  1:37620/1‑149 (MQ=255)                                                                                                     cTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGc                                             >  3:164343/1‑149 (MQ=255)                                                                                                           ccGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCt                                       >  8:305507/1‑149 (MQ=255)                                                                                                           ccGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCt                                       <  6:74124/149‑1 (MQ=255)                                                                                                            cGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTg                                      <  7:199339/149‑1 (MQ=255)                                                                                                             gCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCAGTGTAGCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGc                                        >  1:268628/1‑149 (MQ=255)                                                                                                             gCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGa                                     <  7:237188/149‑1 (MQ=255)                                                                                                              cTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGAt                                    >  4:178800/1‑149 (MQ=255)                                                                                                              cTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGAt                                    <  8:225366/149‑1 (MQ=255)                                                                                                              cTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGAt                                    <  6:67608/149‑1 (MQ=255)                                                                                                                 gAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCaa                                 <  4:47510/149‑1 (MQ=255)                                                                                                                  aaGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCaac                                >  3:11452/1‑149 (MQ=255)                                                                                                                       cTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCaacgaacg                           >  3:48008/1‑149 (MQ=255)                                                                                                                            aTGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCGACGAACGCCTg                       <  6:142308/148‑1 (MQ=255)                                                                                                                               gCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCAACGAACGCCTGGac                    >  5:29911/1‑148 (MQ=255)                                                                                                                                ccAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCAACGAACGCCTGGacaa                  <  4:11452/149‑1 (MQ=255)                                                                                                                                    gCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCAACGAACGCCTGGACAACATg              <  4:164343/149‑1 (MQ=255)                                                                                                                                         cTGGTGTCGGCCACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCAACGAACGCCTGGACAACATGGCCAg         >  8:188136/1‑149 (MQ=255)                                                                                                                                          tggtggCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCaacga                              <  6:39157/127‑1 (MQ=255)                                                                                                                                          tggtggCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCaacga                              >  5:39157/1‑127 (MQ=255)                                                                                                                                                ccGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGTCCAGGAACTGGCCCGGCTGATCAACGAACGCCTGGACAACATGGCCAGCGAAGTc  <  5:230506/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      AGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCAACGAACGCCTGGACAACATGGCCAGCGAAGTC  >  NC_002947/6042360‑6042651

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