Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,463,424	Δ1 bp	5.7%	coding (758/1116 nt)	algL ←	alginate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,463,424	0	T	.	5.7%	219.6 / 6.9	70	coding (758/1116 nt)	algL	alginate lyase
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base . (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

GGGTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCCCAGTAGCTG  >  NC_002947/1463277‑1463570                                                                                                                                                     |                                                                                                                                                   ggggCGCGGGCGCCGGCCTGCAGGCCTTCAGCCAGGCGCTGCAAGGCGTAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGTGCGAAGCTGGCATTCATCGGCAGCGGCGGCAAGTCGTAGTTGTGATAGGCCAGC‑GCGCgct                                                                                                                                                    <  8:259787/148‑1 (MQ=255) gggTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCgct                                                                                                                                                    <  5:246836/148‑1 (MQ=255)    tCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGctgctg                                                                                                                                                >  7:220509/1‑149 (MQ=255)    tCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGctgctg                                                                                                                                                >  5:343845/1‑149 (MQ=255)       cgGGCGCCGGCCAGCACGCCTTCACACAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGGGATAGGCCAGC‑GCGCGCTGCTGgcg                                                                                                                                             >  8:362235/1‑149 (MQ=255)        gggCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGCAGCGCGC‑GCTGgcct                                                                                                                                            <  4:247405/149‑3 (MQ=255)           cgcCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGCAGCGCGC‑GCTGgcctatc                                                                                                                                         <  8:94165/149‑6 (MQ=255)           cgcCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGATGGCGCTTg                                                                                                                                         >  5:260997/1‑149 (MQ=255)            gcCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGCAGCGCGC‑GCTGgcctatca                                                                                                                                        >  6:45524/1‑143 (MQ=255)            gcCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGa                                                                                                                                        <  1:184358/149‑1 (MQ=255)              cGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGAtt                                                                                                                                      >  8:27065/1‑149 (MQ=255)                  cAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGt                                                                                                                                  <  5:283618/149‑1 (MQ=255)                     cACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTgg                                                                                                                               >  2:133474/1‑149 (MQ=255)                      aCGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTggg                                                                                                                              <  7:220450/149‑1 (MQ=255)                        gCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCg                                                                                                                                   <  7:184913/142‑1 (MQ=255)                        gCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCg                                                                                                                                   >  8:184913/1‑142 (MQ=255)                        gCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTggg                                                                                                                              >  4:260304/1‑147 (MQ=255)                        gCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCa                                                                                                                            <  8:220509/149‑1 (MQ=255)                             cACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAg                                                                                                                       >  3:108915/1‑149 (MQ=255)                                ccAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGccc                                                                                                                    >  3:266726/1‑149 (MQ=255)                                ccAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGccc                                                                                                                    <  8:191996/149‑1 (MQ=255)                                ccAGGCGCAGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGccc                                                                                                                    >  5:82119/1‑149 (MQ=255)                                  aGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTg                                                                                                                  <  3:179818/149‑1 (MQ=255)                                  aGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTg                                                                                                                  <  1:33534/149‑1 (MQ=255)                                   ggCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTgg                                                                                                                 <  6:190086/149‑1 (MQ=255)                                   ggCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTgg                                                                                                                 <  5:126228/149‑1 (MQ=255)                                   ggCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTgg                                                                                                                 <  6:55206/149‑1 (MQ=255)                                      gcTGCAAGGCGAAGTTGTTCTCGGCGAGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCa                                                                                                              <  5:69239/149‑1 (MQ=255)                                          cAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATcc                                                                                                           <  7:312533/148‑1 (MQ=255)                                            aGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCAcc                                                                                                        <  6:343845/149‑1 (MQ=255)                                            aGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCAcc                                                                                                        >  5:214413/1‑149 (MQ=255)                                            aGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCAcc                                                                                                        <  1:261731/149‑1 (MQ=255)                                            aGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCAGCAAGGCGAAGTTGTGATAGGCCAGC‑GCGCTCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCCTCCAcc                                                                                                        >  7:258310/1‑149 (MQ=255)                                                  aGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTg                                                                                                  >  4:38000/1‑149 (MQ=255)                                                       ttCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGcc                                                                                             <  3:260304/149‑1 (MQ=255)                                                         ctcGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGa                                                                                           >  3:128821/1‑149 (MQ=255)                                                          tcGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGAc                                                                                          >  6:314156/1‑149 (MQ=255)                                                            ggCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGa                                                                                           <  7:33223/146‑1 (MQ=255)                                                            ggCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGAcc                                                                                         <  4:266726/148‑1 (MQ=255)                                                                   aGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGATTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTc                                                                                 <  2:266417/149‑1 (MQ=255)                                                                     gTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCt                                                                               >  1:233513/1‑149 (MQ=255)                                                                     aactacgccttgccgccgctggcgatgattgccagcttcgcccaggccaatggcgtAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTAGTATTCCt                                                                               >  7:210583‑M2/57‑149 (MQ=255)                                                                       ctacgccttgccgccgctggcgatgattgccagcttcgcccaggccaatggcgtAGTTGTGTTAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTgg                                                                                                                 <  5:275836‑M2/59‑1 (MQ=255)                                                                       ctacgccttgccgccgctggcgatgattgccagcttcgcccaggccaatggcgtAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTgg                                                                                                                 >  6:275836‑M2/55‑113 (MQ=255)                                                                             CATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGcc                                                                       >  5:314156‑M1/1‑149 (MQ=255)                                                                                      gggCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATc                                                                                                            >  6:210268/1‑103 (MQ=255)                                                                                      gggCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATc                                                                                                            <  5:210268/103‑1 (MQ=255)                                                                                       ggCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGaa                                                               >  8:107119/1‑147 (MQ=255)                                                                                       ggCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAAc                                                              <  7:27065/148‑1 (MQ=255)                                                                                         cGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAgg                                                           >  6:195660/1‑149 (MQ=255)                                                                                         cGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAgg                                                           >  7:77225/1‑149 (MQ=255)                                                                                             gCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTcgc                                                       >  7:336955/1‑149 (MQ=255)                                                                                                ggCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGAc                                                    >  1:347596/1‑149 (MQ=255)                                                                                                ggCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGAc                                                    >  3:184180/1‑149 (MQ=255)                                                                                                 gCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGa                                                     <  1:87348/147‑1 (MQ=255)                                                                                                 gCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGAcg                                                   >  3:43074/1‑149 (MQ=255)                                                                                                   aaTCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGAcggt                                                 <  1:133474/149‑1 (MQ=255)                                                                                                   aaTCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGAcggt                                                 <  4:128821/149‑1 (MQ=255)                                                                                                             aGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGGCGCGACGGTCGGTGGCGAc                                       >  6:270098/1‑149 (MQ=255)                                                                                                               cggcggCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGAcgg                                     <  7:107119/149‑1 (MQ=255)                                                                                                                 gcggcAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATc                                                                                                            <  1:165748/76‑1 (MQ=255)                                                                                                                 gcggcAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATc                                                                                                            >  2:165748/1‑76 (MQ=255)                                                                                                                  cggcAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGAcggcgg                                  <  5:270098/149‑1 (MQ=255)                                                                                                                    gcAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTa                                <  4:43074/149‑1 (MQ=255)                                                                                                                      aaGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGc                              >  7:143416/1‑149 (MQ=255)                                                                                                                      aaGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGc                              >  8:356394/1‑149 (MQ=255)                                                                                                                         gcgtAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCAt                           <  7:362235‑M2/145‑1 (MQ=255)                                                                                                                           gtAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGa                         >  2:346217‑M2/3‑149 (MQ=255)                                                                                                                                 gtgATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATgaccgac                    <  2:347596/148‑1 (MQ=255)                                                                                                                                   gATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGcc                             >  5:323408/1‑137 (MQ=255)                                                                                                                                   gATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGcc                             <  6:323408/137‑1 (MQ=255)                                                                                                                                      gcgCCAGCAGCGCGC‑GCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGgcc              >  3:149091/3‑149 (MQ=255)                                                                                                                                        gCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGgccgc            >  1:129647/1‑149 (MQ=255)                                                                                                                                                ggcgcTGCGGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCGGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGGCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCCCAGTAg     <  3:336486/148‑1 (MQ=255)                                                                                                                                                   gctgctGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCCCAGTAGCtg  <  4:184180/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGGTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCCCAGTAGCTG  >  NC_002947/1463277‑1463570

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