Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,463,424	Δ1 bp	8.8%	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	.	8.8%	170.8 / 10.8	57	coding (758/1116 nt)	algL	alginate lyase
Reads supporting (aligned to +/- strand):  ref base T (21/31);  new base . (2/3);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GGGTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCCCAGTAGCTGTG  >  NC_002947/1463277‑1463572                                                                                                                                                     |                                                                                                                                                     gggTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCgctg                                                                                                                                                     >  7:192131/1‑149 (MQ=255)  ggTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCgctgc                                                                                                                                                    >  5:78028/1‑149 (MQ=255)   gTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCgctgct                                                                                                                                                   >  8:120754/1‑149 (MQ=255)      gcgGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGgc                                                                                                                                                >  1:76778/1‑149 (MQ=255)           cgcCGGCCAGCACGCCTTCACCCATGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTAGAAGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTg                                                                                                                                           <  6:267299/149‑1 (MQ=255)           cgcCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGCAGCGCGC‑GCTGgcctatc                                                                                                                                           <  3:75906/149‑6 (MQ=255)            gcCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGa                                                                                                                                          <  7:120754/149‑1 (MQ=255)              cGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGAtt                                                                                                                                        <  4:176668/149‑1 (MQ=255)              cGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGAtt                                                                                                                                        <  5:23438/149‑1 (MQ=255)                  cAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGt                                                                                                                                    <  6:229885/149‑1 (MQ=255)                   aGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGtt                                                                                                                                   <  8:241842/149‑1 (MQ=255)                     cACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTgg                                                                                                                                 <  2:297548/149‑1 (MQ=255)                      aCGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTggg                                                                                                                                <  7:187651/149‑1 (MQ=255)                        gCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTggg                                                                                                                                >  1:276372/1‑147 (MQ=255)                          cTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGCAGCGCGC‑GCTGgcctat                                                                                                                                            >  6:42637/1‑129 (MQ=255)                          cTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGCAGCGCGC‑GCTGgcctat                                                                                                                                            <  5:42637/133‑5 (MQ=255)                               cccAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGcc                                                                                                                       <  4:27946/149‑1 (MQ=255)                                ccAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGccc                                                                                                                      >  3:297511/1‑149 (MQ=255)                                ccAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGccc                                                                                                                      >  5:33187/1‑149 (MQ=255)                                ccAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGccc                                                                                                                      >  8:26676/1‑149 (MQ=255)                                 cAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCt                                                                                                                     <  6:78028/149‑1 (MQ=255)                                       cTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCAt                                                                                                               <  2:276372/149‑1 (MQ=255)                                        tGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGcc                                                                                                                       <  6:239085/140‑1 (MQ=255)                                        tGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGcc                                                                                                                       >  5:239085/1‑140 (MQ=255)                                        tGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATc                                                                                                              <  7:114393/149‑1 (MQ=255)                                                gAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGCAGCGCGC‑GCTGgcctatc                                                                                                                                           <  8:88249/112‑6 (MQ=255)                                                gAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGCAGCGCGC‑GCTGgcctatc                                                                                                                                           >  7:88249/1‑107 (MQ=255)                                                gAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGt                                                                                                      >  8:282029/1‑149 (MQ=255)                                                 aaGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGtt                                                                                                     >  4:255804/1‑149 (MQ=255)                                                 aaGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGtt                                                                                                     <  8:5534/149‑1 (MQ=255)                                                        tctcGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGgcgc                                                                                                <  1:112003/147‑1 (MQ=255)                                                        tctcGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGgcgc                                                                                                >  2:112003/1‑147 (MQ=255)                                                               gcgcAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGt                                                                                       <  7:282029/149‑1 (MQ=255)                                                                     ggcGACGCCATTGGCCTGGGCGAAGCGGGCAAGCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCt                                                                                 <  7:265805/147‑1 (MQ=255)                                                                     gTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCt                                                                                 <  7:261063/149‑1 (MQ=255)                                                                        gACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTg                                                                                        >  8:86530/1‑139 (MQ=255)                                                                        gACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTg                                                                                        <  7:86530/139‑1 (MQ=255)                                                                        gACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCa                                                                              >  8:161663/1‑149 (MQ=255)                                                                        gACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCa                                                                              <  7:26676/149‑1 (MQ=255)                                                                              aTTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCACC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACggac                                                                         >  3:302929/1‑146 (MQ=255)                                                                                         cGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAgg                                                             <  4:197635/149‑1 (MQ=255)                                                                                                 gCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTa                                                                                      <  4:153978/116‑1 (MQ=255)                                                                                                 gCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTa                                                                                      >  3:153978/1‑116 (MQ=255)                                                                                                 gCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGAcg                                                     >  8:56909/1‑149 (MQ=255)                                                                                                 gCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGAc                                                      <  8:10434/148‑1 (MQ=255)                                                                                                    atcatcGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGAcggtc                                                  <  7:161663/149‑1 (MQ=255)                                                                                                        tcGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGtt                                                                                                                                   >  8:11424/1‑64 (MQ=255)                                                                                                        tcGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGtt                                                                                                                                   <  7:11424/64‑1 (MQ=255)                                                                                                           ccAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGc                                            >  6:176373/1‑148 (MQ=255)                                                                                                           ccAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCg                                           >  2:151540/1‑149 (MQ=255)                                                                                                             aGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGAc                                         <  3:255804/149‑1 (MQ=255)                                                                                                             aGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCCACCGACCGTCGCGAc                                         <  2:76778/149‑1 (MQ=255)                                                                                                                    gcAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGa                                                       >  1:214961/1‑128 (MQ=255)                                                                                                                    gcAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGa                                                       <  2:214961/128‑1 (MQ=255)                                                                                                                    gcAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTa                                  <  7:56909/149‑1 (MQ=255)                                                                                                                                    aTAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAAAAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAgg                  >  7:205284/1‑149 (MQ=255)                                                                                                                                        gCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGgccgc              >  7:281288/1‑149 (MQ=255)                                                                                                                                         ccAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGgccgcc             <  1:151540/149‑1 (MQ=255)                                                                                                                                         ccAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGgccgcc             >  4:182005/1‑149 (MQ=255)                                                                                                                                                  cgcTGCGGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCCCAGTAGCt     <  8:192131/149‑1 (MQ=255)                                                                                                                                                     tgctgGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCCCAGTAGCtgtg  >  4:11537/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGGTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCCCAGTAGCTGTG  >  NC_002947/1463277‑1463572

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