BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,469,125	T→G	7.6%	M5L (ATG→CTG)	algL ←	alginate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,469,125	0	T	G	7.6%	862.1 / 10.1	294	M5L (ATG→CTG)	algL	alginate lyase
Reads supporting (aligned to +/- strand): ref base T (146/118); new base G (2/28); total (148/146)
Fisher's exact test for biased strand distribution p-value = 1.12e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

CGCTTCGCAGCGGAAATTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTT  >  NC_002947/1468976‑1469273

cgcTTCGCAGCGGAAATTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCAt                                                                                                                                                       <  1:1182769/150‑1 (MQ=255)

  cTTCGCAGCGGATATTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATgg                                                                                                                                                     <  2:1516008/150‑1 (MQ=255)

  cTTCGCAGCGGAAATTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATgg                                                                                                                                                     <  1:1630950/150‑1 (MQ=255)

   ttCGCAGCGGAAATTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATggg                                                                                                                                                    <  1:5997843/150‑1 (MQ=255)

      gCAGCGGAAATTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCg                                                                                                                                                 >  1:2945972/1‑150 (MQ=255)

           ggAAATTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCa                                                                                                                                            <  1:3349692/150‑1 (MQ=255)

             aaaTTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTGTGAAAAGGGTCATGGGTCGCCTCACt                                                                                                                                          >  2:373557/1‑150 (MQ=255)

               aTTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTgg                                                                                                                                        <  2:5898525/150‑1 (MQ=255)

               aTTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTgg                                                                                                                                        <  2:4415259/150‑1 (MQ=255)

                ttACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGgc                                                                                                                                       >  1:5026003/1‑150 (MQ=255)

                ttACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGgc                                                                                                                                       <  1:5223689/150‑1 (MQ=255)

                    cATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGgcc‑tgc                                                                                                                                  <  1:943852/150‑1 (MQ=255)

                    cATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGgcc‑tgc                                                                                                                                  >  2:3508681/1‑150 (MQ=255)

                    cATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGgcc‑tgc                                                                                                                                  >  1:3983133/1‑150 (MQ=255)

                     aTCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGgcc‑tgcc                                                                                                                                 <  1:1841465/150‑1 (MQ=255)

                      tCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGgcc‑tgcct                                                                                                                                >  1:4735498/1‑150 (MQ=255)

                          cGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTgcgc                                                                                                                            >  2:132174/1‑150 (MQ=255)

                               ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCt                                                                                                                       >  2:7295746/1‑150 (MQ=255)

                                tGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTg                                                                                                                      >  2:2404971/1‑150 (MQ=255)

                                 gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTgc                                                                                                                     >  2:5328698/1‑150 (MQ=255)

                                 gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTgc                                                                                                                     >  1:3413797/1‑150 (MQ=255)

                                 gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTgc                                                                                                                     <  1:2293638/150‑1 (MQ=255)

                                   gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTgcgc                                                                                                                   <  1:4243082/150‑1 (MQ=255)

                                    cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAACGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTAATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCa                                                                                                                  <  2:7126555/150‑1 (MQ=255)

                                     ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAg                                                                                                                 >  2:1731548/1‑150 (MQ=255)

                                     ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAg                                                                                                                 >  2:7192416/1‑150 (MQ=255)

                                     ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAg                                                                                                                 >  2:7040886/1‑150 (MQ=255)

                                        tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCa                                                                                                              <  2:6201952/150‑1 (MQ=255)

                                        tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCAGACCGCGCAGGCa                                                                                                              <  2:6110965/150‑1 (MQ=255)

                                          gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCGTGCCTGCGCAGACCGCGCAGGCag                                                                                                             <  1:4128998/150‑2 (MQ=255)

                                          gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCGTGCCTGCGCAGACCGCGCAGGCag                                                                                                             <  2:5612959/150‑2 (MQ=255)

                                           aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGGGCGCGGCACCACCAACCACGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                           <  2:2818496/150‑1 (MQ=255)

                                            tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGc                                                                                                          >  2:3364705/1‑150 (MQ=255)

                                             ccccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCAGACCGCGCAGGCaggcca                                                                                                         <  2:5946753/150‑6 (MQ=255)

                                               ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                       >  1:6102510/1‑150 (MQ=255)

                                               ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                       >  2:1009454/1‑150 (MQ=255)

                                                cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCCCGGCGAAg                                                                                                      >  1:5765704/1‑150 (MQ=255)

                                                cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                      >  2:3912391/1‑150 (MQ=255)

                                                cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                      >  1:4338884/1‑150 (MQ=255)

                                                  ggTAGTAACCCTGCGGCGGTACCAGCGGGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGCGCCGGGGGGGGCTGTAATTTGTTTTGAAATGGGTCTTGGGTCGCCTCACTGGCC‑TGCCTGCGGGGTCTGCGCAGGCACGGCGAAgtt                                                                                                    <  2:1812217/150‑1 (MQ=255)

                                                    tagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCCGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCAGGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgt                                                                                                  >  2:5237667/1‑150 (MQ=255)

                                                     tgtaCCCCGGGGGCGGTCCCCGCGCGGCGTGCGCGGCACCACCAAACACGGCCACCGCGAGCCGGGGGGGGCTGTAAATTCTTTTGAAATGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgtt                                                                                                 <  2:1203840/149‑1 (MQ=255)

                                                         aCCCTGCGGGGGTACCAGCGCGGCGTGCGCGGCACAACAAAACAAGCCAAACGCGAGCAGGGGGGGGCTGAAAATTCTTTTTAAAAGGGTATTGGGTCGCCTCACTGGCC‑TGCCTGCGCGTTCTGCGCAGGCACGGCGAAGTTGTTGCGt                                                                                             <  2:3952112/150‑1 (MQ=255)

                                                          cccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGGCTGCGCACGCACGGCGAAGGTGTTGCgat                                                                                            >  1:3733044/1‑148 (MQ=255)

                                                           ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGCGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                           >  2:2896076/1‑150 (MQ=255)

                                                            cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTg                                                                                          >  2:7184781/1‑150 (MQ=255)

                                                             tGCGGCGGTACCCGCGCGGCGTGCGCGGCACCACCAACCACGGCCAACGCGAGCTGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑GGCCTGCGCGTTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGc                                                                                         <  2:7432966/150‑1 (MQ=255)

                                                             tGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGc                                                                                         >  2:2960163/1‑150 (MQ=255)

                                                             tGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGc                                                                                         >  1:2574539/1‑150 (MQ=255)

                                                             gggggcggTACCCGCGCGGGGTGCGCGCCACCACAAACCCCGCCCACCGCGCGCTGGGCGGGGCGGCAAATTCTTTTCAAATGGGTCTGGGGTCGCCTCCGTGGCC‑GGCGTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGc                                                                                         <  2:1947109/147‑1 (MQ=255)

                                                              gcggcgGTACCAGCGCGGCGTGCGCGGCACCACCAAACACGGCCAACGCGAGCCGGGGGGGGCTGGAAATTCGTTTGAAAAGGGTCTTGGGTCGCCGCACGGGCC‑GGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCa                                                                                        <  2:5987666/150‑1 (MQ=255)

                                                               ggggggTACCAGGGGGGGGGGCGCGGCACCACCAAACACGGCAAACGCTTGCTGGGCGGGGCTGTAAATTGTTTTTAAATGGGTCTTGGGTCGCCTCACTGGCC‑TGCCTGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCag                                                                                       <  1:5508129/146‑1 (MQ=255)

                                                                    gTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTTCAGAGCtt                                                                                  >  1:2484551/1‑150 (MQ=255)

                                                                      aCCTGCGCGGCGTGCGCGGCACCACAAACCATGGCCACCGCGAGCTGGGCGGGGCTGGAAATTCTTTTTAAATGGGTCTGGGGCCGCCTCACTGGCG‑TGCCTGCGCGTTCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTgg                                                                                <  1:4926696/150‑1 (MQ=255)

                                                                       cccgggcggggGGCGCGGACCCACCAAACACGCCAACCGGTCGCGGGGGGGGGCTGTAATTTGGTTTCAATTGGGTCTTGGGTCGCCTCAGTGGCG‑GGCCTGGGGGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGc                                                                               <  2:3983139/145‑1 (MQ=255)

                                                                       cccgcgcGGCGTGCGCGGCACCCCCAAACACGGCCACCGCGCGCCGGGCGGGGCTGCAAATTCTTTTGAAACGGGTCTTGGGTCGCCTCACTGGCC‑TGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGc                                                                               <  2:1085581/147‑1 (MQ=255)

                                                                        cAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCt                                                                              >  1:3165567/1‑150 (MQ=255)

                                                                          gcgcGGGGGGCGCGGCACCACCAAACACGGCCACCGCGCGCGGGGGGGGGCTGGAAATTGTTTTGAAATGGGTCAGGGGTCGCCTCCCGGGCC‑GGCCTGCGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                            <  2:1599123/150‑1 (MQ=255)

                                                                            gcggcgTGCGCGCCCCCACCAAACAAGGCCAACGCGTGCCGGGGGGGGCTGTAATTTCGTTTGAAATGGGTCATGGGTCGCCTCACTGGCC‑GGCGTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGa                                                                          <  2:4562671/150‑1 (MQ=255)

                                                                             cggcgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCGGCGCGGTCTGCGCAGGCACGGCGCAGGTGTTGCGTTTGCAGAGCTTGGCTTCGAc                                                                         >  1:3542277/1‑150 (MQ=255)

                                                                             cggcgTCCGCGCCACCACCAAACACGGCAACCCCGAGCAGGGCGGGGCTGCAAATTCGTTTTAAAAGGGTCAGGGGTCGCCTCACTGGCC‑GGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGAc                                                                         <  2:4901824/150‑1 (MQ=255)

                                                                                cgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCtt                                                                      >  2:6708426/1‑150 (MQ=255)

                                                                                    cgcgGCACCCCAAAACAAGGCCAACGCGAGCTGGGCGGGGCTGTAAATTCGTTTGAAAAGGGTCTTGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGGGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                  <  2:5537958/150‑1 (MQ=255)

                                                                                    cgcgGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                  >  2:2927825/1‑150 (MQ=255)

                                                                                     gcgGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCCGAGCTTGGCTTCGACCTTCTggg                                                                 >  1:1950360/1‑150 (MQ=255)

                                                                                         caccacCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGGTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGcc                                                             >  2:2149972/1‑150 (MQ=255)

                                                                                          accaccAAACCAGGCCACCGCGAGCCGGGCGGGGCTGTAATTTCGTTTGAAACGGGTCATGGGTCGCCTCCCTGGCC‑GGCGTGCGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCg                                                            <  2:6333830/150‑1 (MQ=255)

                                                                                          accacaAAACACGCCCACCGCGCGCGGGGGGGGGCTGTAATTTCTTTTCAAATGGGACATGGGTCGCCTCACGGGCC‑GGCCGGCGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCg                                                            <  2:2584614/150‑1 (MQ=255)

                                                                                            caccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGACACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCt                                                          >  1:1481051/1‑150 (MQ=255)

                                                                                               caAACATGGCCAACGCGAGCCGGGCGGGGCTGGAAATTCTTTTTAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTc                                                       <  2:4633794/150‑1 (MQ=255)

                                                                                                aaaCACGGCCACCGCTAGCCGGGGGGGGCTGTAATTTCGTTTTAAATGGGCCTTGGGCCGCATCACGGGCC‑GGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                      <  1:3020791/150‑1 (MQ=255)

                                                                                                aaaCAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGGTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                      >  2:6856196/1‑150 (MQ=255)

                                                                                                    aaGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGAAGGCACTGAGAAGTTGTTGCTTTTGCAGATCTTGGCTTCGAACTTCTGGGTGACGCTTTCCGGcc                                                  >  2:6410232/1‑150 (MQ=255)

                                                                                                     aGGCCAACGCGCGCCGGGCGGGGCTGTAACTTCGTTTTAAACGGGTCCTGGGCCGCCTCACGGGCC‑TGCCGGCGGGGTGTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                 <  2:1266567/150‑1 (MQ=255)

                                                                                                       gCCAACGCGAGCCGGGGGGGGCTGGAAATTCGTTTGAAATGGGTCCTGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCt                                               <  1:2169367/150‑1 (MQ=255)

                                                                                                         caccgcgcgcgGGGGGGGGCTGTAAATTCGTTTGAAACGGGCCTTGGGTCGCCGCAGGGGCC‑GGCGTGCGCGGTGTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGt                                             <  2:3111851/147‑1 (MQ=255)

                                                                                                          aaCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTa                                            >  1:6986081/1‑150 (MQ=255)

                                                                                                           aCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAc                                           >  2:360774/1‑150 (MQ=255)

                                                                                                           aCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAc                                           >  1:1645929/1‑150 (MQ=255)

                                                                                                              cgcgctggggggggCTTGAAATTCGTTTTAAATGGGTATTGGGTCGCCTCACGGGCC‑TGCCTGCGCGGTGTGGGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                        <  2:5742267/144‑1 (MQ=255)

                                                                                                              cgAGCTGGGCGGGGCGGTAAATGCGTTTGAAAAGGGTCTTGGGGCGCCTCACGGGCC‑GGCCTGCGCGGTCGGCTCTGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                        <  2:4133715/150‑1 (MQ=255)

                                                                                                                  cgggggggggCTGGAATTGCGTTTTAAATGGGTCTGGGGGCGCCTCCCGGGGG‑GGCGGGGGGGGTGGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAc                                    <  1:767698/148‑1 (MQ=255)

                                                                                                                   tgggCGGGGCTGGAATTGCGTTTGAAACGGGTCTGGGGGCGCCTCACTGGCC‑TGCCTGCGCGGTCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcg                                   <  2:5035063/149‑1 (MQ=255)

                                                                                                                      ggggggGTGGAATTTCGTTTGAAATGGGTCATGGGGCGCCTCACTGGCC‑GGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                <  2:5739472/148‑1 (MQ=255)

                                                                                                                      gCGGGGCGGTAAATTCTTTTTAAAAGGGTCTTGGGTCGCCTCACGGGCC‑TGCGTGCGGGGTCTGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                <  1:1093393/150‑1 (MQ=255)

                                                                                                                        ggggCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGCAGGTGTTGCGTTTGCCGAGCTTGGGTTCGACCGTCTGGGGGCCGCGTTCCGGCCCCGGTACCGCGAACGCGAGc                              >  1:6838120/1‑150 (MQ=255)

                                                                                                                             tGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGGTGTTGCGTTTGCAGAGCTTGGCTTCGACCTGCTGGGTGCCGCTTTCCGGCCCCGGTACCACGAACGCGAGCAGgcg                         >  2:1891015/1‑150 (MQ=255)

                                                                                                                              ggAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCAGCGCGGGCTGCGCAGGCACGGCGAAGTTGTTTCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCGTTCCGGCCCCAGTAACTCGAACGCGAGCAGgcgc                        >  1:3836460/1‑150 (MQ=255)

                                                                                                                                       tttGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTTCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTACGGCCCCTGTAACTCGAACGCGAGCAGGCGCAGGCTGGcc               >  1:6456269/1‑150 (MQ=255)

                                                                                                                                       tttGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCGTCTGGGTGCCGCTTTCCGGCCCCTGTACCACGAAAGCGAGCAGGCGCTGGCTGGcc               >  1:4413023/1‑150 (MQ=255)

                                                                                                                                       tttGAAAAGGGTAATGGGTCTACTCACTGGCC‑TGCAAGGGCGGTCAGAGAAGGAAAGGACAATTTGTTTGGTTTGAAGATCTTTGATTCGACCTTCTGGGTGCCGCTTTAAGGACCCTGTACCTCGAAAGCGAGAAGGAgctggcagtcc               >  2:6271735/1‑145 (MQ=255)

                                                                                                                                         ttaattgggTCTGGGGCCGCCTCCGTGGCC‑GGCCGGCGGGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa             <  1:67988/144‑1 (MQ=255)

                                                                                                                                         ttaaatGGGTCCGGGGCCGCCCCCCGGGCC‑TGCGGGCGGGGCCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa             <  1:1731548/148‑1 (MQ=255)

                                                                                                                                         tgaactgggTCTGGGGCCGCCACCCTGGCC‑GGCGGGGGGGGTGTGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa             <  1:4075278/144‑1 (MQ=255)

                                                                                                                                          gAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAg            >  2:5466076/1‑150 (MQ=255)

                                                                                                                                             aaGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGGAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCAGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTcc         >  1:6908721/1‑150 (MQ=255)

                                                                                                                                             aaGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTcc         >  1:6203385/1‑150 (MQ=255)

                                                                                                                                             aaGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCAGACCGCGCAGGCag                                                                                                             <  2:6613954/50‑2 (MQ=21)

                                                                                                                                             aaGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCAGACCGCGCAGGCag                                                                                                             >  1:6613954/1‑49 (MQ=21)

                                                                                                                                               gggTCTTGGGCCGCCCCCCGGGCC‑GGCCTGCGGGGTCTGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc       <  1:3508681/150‑1 (MQ=255)

CGCTTCGCAGCGGAAATTACCATCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCC‑TGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTT  >  NC_002947/1468976‑1469273

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