Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,463,417:1	+A	8.5%	coding (765/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,417	1	.	A	8.5%	223.1 / 15.7	59	coding (765/1116 nt)	algL	alginate lyase
Reads supporting (aligned to +/- strand):  ref base . (24/30);  new base A (2/3);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCTGGGGTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCC  >  NC_002947/1463270‑1463561                                                                                                                                                     |                                                                                                                                                 gTGGCTGGGGTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTGGTGATAGGCCAgc‑g                                                                                                                                                 >  3:197635/1‑149 (MQ=255)    gCTGGGGTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAgc‑gcgc                                                                                                                                              >  1:297548/1‑149 (MQ=255)      tGGGGTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGcagcgcg                                                                                                                                             >  4:75906/1‑149 (MQ=255)        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                                                                                                                                 <  5:42637/133‑5 (MQ=255)                                 cTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGCAGCGCGC‑GCTGgcctat                                                                                                                                 >  6:42637/1‑129 (MQ=255)                                      cccAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGcc                                                                                                            <  4:27946/149‑1 (MQ=255)                                       ccAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGccc                                                                                                           >  5:33187/1‑149 (MQ=255)                                       ccAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGccc                                                                                                           >  3:297511/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                                                                             <  7:86530/139‑1 (MQ=255)                                                                               gACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTg                                                                             >  8:86530/1‑139 (MQ=255)                                                                               gACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCa                                                                   <  7:26676/149‑1 (MQ=255)                                                                               gACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCa                                                                   >  8:161663/1‑149 (MQ=255)                                                                                     aTTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCACC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACggac                                                              >  3:302929/1‑146 (MQ=255)                                                                                                cGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAgg                                                  <  4:197635/149‑1 (MQ=255)                                                                                                        gCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTa                                                                           >  3:153978/1‑116 (MQ=255)                                                                                                        gCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTa                                                                           <  4:153978/116‑1 (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                                                                                                                        <  7:11424/64‑1 (MQ=255)                                                                                                               tcGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGtt                                                                                                                        >  8:11424/1‑64 (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)                                                                                                                                                     |                                                                                                                                                 GTGGCTGGGGTCGCGGGCGCCGGCCAGCACGCCTTCACCCAGGCGCTGCAAGGCGAAGTTGTTCTCGGCGCGCAGGTCGACGCCATTGGCCTGGGCGAAGCTGGCAATCATCGCCAGCGGCGGCAAGGCGTAGTTGTGATAGGCCAGC‑GCGCGCTGCTGGCGCTTGATTTCGTTGGGCAGGAAGCCCTGGTCATCCACCTGGTTGGCGCCGACCTTGTATTCCTTCACGGCCCAGTCGAACAGGTCGCGACGGTCGGTGGCGACGGCGGTAGCCATGACCGACCAGGCCGCC  >  NC_002947/1463270‑1463561

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