Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,042,502	G→T	8.6%	A147S (GCC→TCC) ‡	PP_5295 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,042,502	0	G	T	8.6%	181.1 / 6.5	58	A147S (GCC→TCC) ‡	PP_5295	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base T (2/3);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

AGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCAACGAACGCCTGGACAACATGGCCAG  >  NC_002947/6042360‑6042644                                                                                                                                               |                                                                                                                                                   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                                                      >  2:306706/1‑148 (MQ=255)                                                                                      aaGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGcgcg                                                      >  5:67608/1‑148 (MQ=255)                                                                                      aaGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGcgcg                                                      >  8:199339/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                       >  5:39157/1‑127 (MQ=255)                                                                                                                                          tggtggCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCaacga                       <  6:39157/127‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   AGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGA‑ACGCCGAAGCCATGGCGCTGTTCGACGAAGCGC‑‑TG‑GCCGAACTCAAGGCCGGGCGCCAGCGCGAAGGCCAGGAACTGGCCCGGCTGATCAACGAACGCCTGGACAACATGGCCAG  >  NC_002947/6042360‑6042644

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