BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,471,402	G→T	5.1%	R478R (CGG→AGG)	alg8 ←	glycosyltransferase Alg8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,471,402	0	G	T	5.1%	1058.2 / 30.6	311	R478R (CGG→AGG)	alg8	glycosyltransferase Alg8
Reads supporting (aligned to +/- strand): ref base G (115/180); new base T (7/9); total (122/189)
Fisher's exact test for biased strand distribution p-value = 7.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 316 total aligned reads displayed.

GGCGCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCG  >  NC_002947/1471254‑1471551

ggCGCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCcgc                                                                                                                                                      <  2:2634168/150‑1 (MQ=255)

   gcTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAg                                                                                                                                                   >  1:1101215/1‑150 (MQ=255)

   gcTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAg                                                                                                                                                   <  1:5379188/150‑1 (MQ=255)

   gcTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAg                                                                                                                                                   >  2:7131222/1‑150 (MQ=255)

       ggCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGAcc                                                                                                                                               <  2:4778268/150‑1 (MQ=255)

        gCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCa                                                                                                                                              <  2:8657690/150‑1 (MQ=255)

        gCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCa                                                                                                                                              <  2:7331747/150‑1 (MQ=255)

         cTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAg                                                                                                                                             >  1:2140383/1‑150 (MQ=255)

          ttCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAgg                                                                                                                                            <  1:8886071/150‑1 (MQ=255)

          ttCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAgg                                                                                                                                            >  1:3113899/1‑150 (MQ=255)

           tCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGgt                                                                                                                                           >  1:6139076/1‑150 (MQ=255)

            cGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatgac                                                                                                                                          <  2:3934927/150‑5 (MQ=255)

             ggACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGgtgt                                                                                                                                         <  2:1797607/150‑1 (MQ=255)

                 tCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGaa                                                                                                                                     >  1:7165514/1‑150 (MQ=255)

                 tCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGaa                                                                                                                                     <  2:137172/150‑1 (MQ=255)

                 tCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGaa                                                                                                                                     <  1:4242649/150‑1 (MQ=255)

                   aTGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAAcc                                                                                                                                   >  1:8237713/1‑150 (MQ=255)

                   aTGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAAcc                                                                                                                                   >  1:4200084/1‑150 (MQ=255)

                     gCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAt                                                                                                                                 >  1:6746219/1‑150 (MQ=255)

                      cACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                                <  1:2469102/150‑1 (MQ=255)

                       aCGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttg                                                                                                                               >  1:3306843/1‑150 (MQ=255)

                       aCGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGGTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGCAGGTCATGGTCCGCGAGGACCAGGTGGTGTACCAttg                                                                                                                               >  1:7922796/1‑150 (MQ=255)

                        cGACATTGACGTTCACGGCGGTATTCATGGGGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttgt                                                                                                                              >  1:5208868/1‑150 (MQ=255)

                           cATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGa                                                                                                                           <  2:3329065/150‑1 (MQ=255)

                           cATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGa                                                                                                                           >  1:3876407/1‑150 (MQ=255)

                            aTTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCa                                                                                                                                              >  1:1276490/1‑130 (MQ=255)

                            aTTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGaa                                                                                                                          >  2:5826818/1‑150 (MQ=255)

                             ttGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGaaa                                                                                                                         <  1:3272087/150‑1 (MQ=255)

                              tGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGTACCATTGTTGAAAg                                                                                                                        >  1:8432277/1‑150 (MQ=255)

                                  gTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTgg                                                                                                                    >  1:3961290/1‑150 (MQ=255)

                                   ttCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGc                                                                                                                   >  1:8669803/1‑150 (MQ=255)

                                          cggTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCa                                                                                                            >  2:8692968/1‑150 (MQ=255)

                                           ggTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCAc                                                                                                           >  2:6070337/1‑150 (MQ=255)

                                               ttCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCtt                                                                                                       <  1:1754900/150‑1 (MQ=255)

                                                 gaTGGTGGGTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGTAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                     <  1:5183914/149‑1 (MQ=255)

                                                 cATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                     <  2:4046360/150‑1 (MQ=255)

                                                 cATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                     <  2:75461/150‑1 (MQ=255)

                                                 cATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                     <  1:3825811/150‑1 (MQ=255)

                                                 cATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                     <  1:1499387/150‑1 (MQ=255)

                                                   tggtggTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAgg                                                                                                   <  2:3874807/150‑1 (MQ=255)

                                                   tggtggTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAgg                                                                                                   >  1:3444185/1‑150 (MQ=255)

                                                   tggtggTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAgg                                                                                                   <  1:5031193/150‑1 (MQ=255)

                                                    ggtggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAggg                                                                                                  >  2:1463251/1‑150 (MQ=255)

                                                     gtggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGc                                                                                                 >  2:3166023/1‑150 (MQ=255)

                                                      tggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCa                                                                                                <  2:3780855/150‑1 (MQ=255)

                                                        gtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGt                                                                                              <  2:2867796/150‑1 (MQ=255)

                                                        gtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAACACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGt                                                                                              <  2:2766338/150‑1 (MQ=255)

                                                          ttGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTgg                                                                                            <  2:3538346/150‑1 (MQ=255)

                                                            gTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGc                                                                                          <  2:5228436/150‑1 (MQ=255)

                                                             ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCt                                                                                         >  2:1019431/1‑150 (MQ=255)

                                                                cTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAc                                                                                      <  2:2109992/150‑1 (MQ=255)

                                                                  gTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcgc                                                                                    >  2:6648955/1‑150 (MQ=255)

                                                                   ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcgcg                                                                                   >  2:5165900/1‑150 (MQ=255)

                                                                   ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcgcg                                                                                   >  2:3037679/1‑150 (MQ=255)

                                                                     aGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTc                                                                                 <  1:6520150/150‑1 (MQ=255)

                                                                      gatccgatccGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTcc                                                                                <  1:8788562/150‑1 (MQ=255)

                                                                      gatccgatccGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTcc                                                                                <  2:7690022/150‑1 (MQ=255)

                                                                       atccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGGCGCGTCCa                                                                               <  2:8137735/150‑1 (MQ=255)

                                                                          cgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGa                                                                            <  2:5318052/150‑1 (MQ=255)

                                                                           gatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGAc                                                                           >  1:2846317/1‑150 (MQ=255)

                                                                             tccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAATGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGATTGAGGGCAGTGGGCTGACGCGTCCAcaccaa                                                                         <  2:3744765/150‑6 (MQ=255)

                                                                                   ttGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTc                                                                   <  2:6922162/150‑1 (MQ=255)

                                                                                     gggCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGa                                                                 <  2:5695800/150‑1 (MQ=255)

                                                                                       gCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAgg                                                               <  2:1120812/150‑1 (MQ=255)

                                                                                       gCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAgg                                                               >  1:6342264/1‑150 (MQ=255)

                                                                                              cGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaaga                                                        <  1:292940/150‑1 (MQ=255)

                                                                                                aCCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAac                                                      <  1:1449860/150‑1 (MQ=255)

                                                                                                 ccATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAaca                                                     >  1:2324769/1‑150 (MQ=255)

                                                                                                   aTGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACg                                                   >  2:8375260/1‑150 (MQ=255)

                                                                                                   aTGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCGTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACg                                                   >  1:1025699/1‑150 (MQ=255)

                                                                                                     gAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTa                                                 <  2:1653385/150‑1 (MQ=255)

                                                                                                     gAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTa                                                 <  1:1607245/150‑1 (MQ=255)

                                                                                                          gcacagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCt                                            <  2:7796765/150‑1 (MQ=255)

                                                                                                            acagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                          <  2:236234/150‑1 (MQ=255)

                                                                                                            acagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                          <  2:4254818/150‑1 (MQ=255)

                                                                                                                caACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtcatca                                      <  1:4661042/150‑1 (MQ=255)

                                                                                                                caACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtcatca                                      <  1:4660938/150‑1 (MQ=255)

                                                                                                                caACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                          <  1:3502106/146‑1 (MQ=255)

                                                                                                                caACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                          >  2:3502105/1‑146 (MQ=255)

                                                                                                                    gAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAgcgc                                  >  1:7596123/1‑150 (MQ=255)

                                                                                                                            tGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCgatga                          >  2:7630031/1‑150 (MQ=255)

                                                                                                                             ggCAGCCGAGAAGGTCATGGTCCGCGAGGACCGGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCgatgat                         >  1:1610359/1‑150 (MQ=255)

                                                                                                                              gCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATc                        <  1:1867952/150‑1 (MQ=255)

                                                                                                                                aGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTg                      >  2:5355116/1‑150 (MQ=255)

                                                                                                                                   cGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGtt                   <  2:42390/150‑1 (MQ=255)

                                                                                                                                      gaAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTa                <  2:8323367/150‑1 (MQ=255)

                                                                                                                                       aaGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTaa               <  1:4393165/150‑1 (MQ=255)

                                                                                                                                       aaGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTaa               <  2:3616396/150‑1 (MQ=255)

                                                                                                                                        aGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAAt              <  2:1997769/150‑1 (MQ=255)

                                                                                                                                         ggTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATa             <  1:5747650/150‑1 (MQ=255)

                                                                                                                                         ggTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATa             <  1:7215927/150‑1 (MQ=255)

                                                                                                                                            cctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATagag          >  2:7987030/3‑150 (MQ=255)

                                                                                                                                             aTGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAgaga         <  2:6371708/150‑1 (MQ=255)

                                                                                                                                              tGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGaa        <  2:7195577/150‑1 (MQ=255)

                                                                                                                                               ggTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAt       >  2:6246288/1‑150 (MQ=255)

                                                                                                                                                gTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATc      <  1:7719866/150‑1 (MQ=255)

                                                                                                                                                 tCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCa     <  1:456865/150‑1 (MQ=255)

                                                                                                                                                 tCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCa     <  2:882880/150‑1 (MQ=255)

                                                                                                                                                   cgcgAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAGTAGAGAATCAGc   >  1:6439335/1‑150 (MQ=255)

                                                                                                                                                   cgcgAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCAGATGATCTGGTTGTAATAGAGACTCAGc   <  1:6070337/150‑1 (MQ=255)

GGCGCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCG  >  NC_002947/1471254‑1471551

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