Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,463,032	A→T	9.4%	intergenic (‑101/+41)	algI ← / ← algL	alginate O‑acetylase AlgI/alginate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,463,032	0	A	T	9.4%	183.1 / 9.4	64	intergenic (‑101/+41)	algI/algL	alginate O‑acetylase AlgI/alginate lyase
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base T (4/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GGAAGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCCAACG  >  NC_002947/1462883‑1463181                                                                                                                                                      |                                                                                                                                                      ggAAGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCaa                                                                                                                                                       <  2:1740429/150‑1 (MQ=255)   aaGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCCTCCAGGGCAAAGCGCTCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGc                                                                                                                                                     <  1:869309/150‑1 (MQ=255)    aGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGcc                                                                                                                                                    <  2:1307245/150‑1 (MQ=255)    aGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGcc                                                                                                                                                    <  1:1399935/150‑1 (MQ=255)         ggCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTc                                                                                                                                               >  2:1507336/1‑150 (MQ=255)            aaGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCca                                                                                                                                            <  1:1076886/150‑1 (MQ=255)            aaGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCca                                                                                                                                            <  1:93547/150‑1 (MQ=255)                  aggaataggaacACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCGACAAACCCCCCAGTGCAAAGCTGTAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGtt                                                                                                                                      >  1:1106681/1‑150 (MQ=255)                          gaacaCGTTGTACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCca                                                                                                                                            <  1:1065699/136‑1 (MQ=255)                          gaacaCGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCca                                                                                                                                            >  2:1065699/1‑136 (MQ=255)                              acGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCaa                                                                                                                          <  2:1689656/150‑1 (MQ=255)                                gTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCaacg                                                                                                                        >  1:1851006/1‑148 (MQ=255)                                   ggACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATgggg                                                                                                                     <  1:2015520/150‑1 (MQ=255)                                        agaagaCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGAAATGGGTGTGaa                                                                                                                <  2:1625854/150‑1 (MQ=255)                                          aagaCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAACGGGGTTCAAtt                                                                                                              <  2:1851006/150‑1 (MQ=255)                                              ccATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCAtt                                                                                                          <  2:866174/150‑1 (MQ=255)                                               cATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTcc                                                                                                                                    <  2:314616/123‑1 (MQ=255)                                               cATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTcc                                                                                                                                    >  1:314616/1‑123 (MQ=255)                                                      gTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGctc                                                                                                  <  1:340532/150‑1 (MQ=255)                                                               aTTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTc                                                                                         >  2:394711/1‑150 (MQ=255)                                                               aTTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTc                                                                                         >  2:1713957/1‑150 (MQ=255)                                                                    ttCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTc                                                                                    >  1:1515202/1‑150 (MQ=255)                                                                      cacaGCATCCAGGGCAAAGCGCCCCACAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGt                                                                                  <  2:1106681/150‑1 (MQ=255)                                                                        caGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTc                                                                                                                                               <  1:2265068/87‑1 (MQ=255)                                                                        caGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTc                                                                                                                                               >  2:2265068/1‑87 (MQ=255)                                                                             tCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCAt                                                                                                           <  1:1217990/118‑1 (MQ=255)                                                                             tCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCAt                                                                                                           >  2:1217990/1‑118 (MQ=255)                                                                             tCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCcg                                                                           >  2:1480309/1‑150 (MQ=255)                                                                               cAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCAt                                                                                                           <  2:2214551/116‑1 (MQ=255)                                                                               cAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCAt                                                                                                           >  1:2214551/1‑116 (MQ=255)                                                                                   gCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTcc                                                                                                                                              >  2:215632/1‑77 (MQ=255)                                                                                   gCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTcc                                                                                                                                              <  1:215632/77‑1 (MQ=38)                                                                                             ccccAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGAAGCGGGCAAGCCAGCTCCAACGGGCTCCGCGCGCGCTATGGGGAGCAGGTGATGTCGTTCTCGCTGCCCGCTGTGGTTATACCCTCGTCACATCGCCGccc                                                           >  2:1653136/1‑150 (MQ=255)                                                                                              cccAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGCTACGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCa                                                          >  1:747262/1‑150 (MQ=255)                                                                                               ccAAACCCACCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAg                                                         >  1:1026301/1‑150 (MQ=255)                                                                                                 aaaCCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGcc                                                       >  1:1401809/1‑150 (MQ=255)                                                                                                     cccccAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGaa                                                   >  1:820532/1‑150 (MQ=255)                                                                                                         cAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCggg                                                                                      >  1:393000/1‑111 (MQ=255)                                                                                                         cAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCggg                                                                                      <  2:393000/111‑1 (MQ=255)                                                                                                         cAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTg                                               <  2:113815/150‑1 (MQ=255)                                                                                                           gTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGtt                                             <  2:1026301/150‑1 (MQ=255)                                                                                                               aaGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAAc                                         <  1:2317178/150‑1 (MQ=255)                                                                                                                  gCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGc                                      <  2:1515202/150‑1 (MQ=255)                                                                                                                  gCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGc                                      >  1:668741/1‑150 (MQ=255)                                                                                                                     ggAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGc                                   <  2:820532/150‑1 (MQ=255)                                                                                                                         ggCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCt                                     >  1:147509/1‑144 (MQ=255)                                                                                                                         ggCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCt                                     <  2:147509/144‑1 (MQ=255)                                                                                                                            aaCTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGt                            <  1:1507336/150‑1 (MQ=255)                                                                                                                                    tCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTgcgc                    <  1:637281/150‑1 (MQ=255)                                                                                                                                      ccTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTc                  >  1:1941238/1‑150 (MQ=255)                                                                                                                                        tgtgGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCt                                     <  1:2267628/129‑1 (MQ=255)                                                                                                                                        tgtgGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCt                                     >  2:2267628/1‑129 (MQ=255)                                                                                                                                          tgGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGc              <  1:394711/150‑1 (MQ=255)                                                                                                                                           gggAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGaa                                                   <  1:2232787/112‑1 (MQ=255)                                                                                                                                           gggAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGaa                                                   >  2:2232787/1‑112 (MQ=255)                                                                                                                                            ggAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAcgcgc            >  2:2004592/1‑149 (MQ=255)                                                                                                                                            ggAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCgt            >  2:1796309/1‑150 (MQ=255)                                                                                                                                            ggAGCGGGCAAGCGCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACTGCTGCATTTCGTGCTTGCTCTCAAGctt            >  2:135664/1‑148 (MQ=255)                                                                                                                                            ggAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTc                                                                                         <  2:1904981/73‑1 (MQ=255)                                                                                                                                            ggAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTc                                                                                         >  1:1904981/1‑73 (MQ=255)                                                                                                                                             gAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCgtg           >  2:2280011/1‑150 (MQ=255)                                                                                                                                             gAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCgtg           >  1:797481/1‑150 (MQ=255)                                                                                                                                                  ggCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGcc      >  1:1599172/1‑150 (MQ=255)                                                                                                                                                   gCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCCa     >  2:753270/1‑150 (MQ=255)                                                                                                                                                   gCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTACGCTCAAGCGTGTCGCCa     >  1:2056047/1‑150 (MQ=255)                                                                                                                                                     aaGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGTGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACACCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCCAAc   <  1:135664/150‑1 (MQ=255)                                                                                                                                                      aGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCCAACg  <  1:1713957/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GGAAGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCCAACG  >  NC_002947/1462883‑1463181

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑