Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,463,031	A→T	9.4%	intergenic (‑100/+42)	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,031	0	A	T	9.4%	178.1 / 10.3	64	intergenic (‑100/+42)	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 = 5.79e-01

GGAAGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCCAAC  >  NC_002947/1462883‑1463180                                                                                                                                                     |                                                                                                                                                      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                                                                                                          >  1:2214551/1‑116 (MQ=255)                                                                               cAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCAt                                                                                                          <  2:2214551/116‑1 (MQ=255)                                                                                   gCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTcc                                                                                                                                             <  1:215632/77‑1 (MQ=38)                                                                                   gCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTcc                                                                                                                                             >  2:215632/1‑77 (MQ=255)                                                                                             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                                                                                     <  2:393000/111‑1 (MQ=255)                                                                                                         cAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCggg                                                                                     >  1:393000/1‑111 (MQ=255)                                                                                                         cAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTg                                              <  2:113815/150‑1 (MQ=255)                                                                                                           gTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGtt                                            <  2:1026301/150‑1 (MQ=255)                                                                                                               aaGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAAc                                        <  1:2317178/150‑1 (MQ=255)                                                                                                                  gCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGc                                     >  1:668741/1‑150 (MQ=255)                                                                                                                  gCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGc                                     <  2:1515202/150‑1 (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                                                                                        >  1:1904981/1‑73 (MQ=255)                                                                                                                                            ggAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTc                                                                                        <  2:1904981/73‑1 (MQ=255)                                                                                                                                             gAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCgtg          >  1:797481/1‑150 (MQ=255)                                                                                                                                             gAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCgtg          >  2:2280011/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)                                                                                                                                                     |                                                                                                                                                      GGAAGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCCAAC  >  NC_002947/1462883‑1463180

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