BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,781,099	A→C	7.0%	K152Q (AAG→CAG)	sotB →	sugar efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,781,099	0	A	C	7.0%	1012.1 / 37.4	337	K152Q (AAG→CAG)	sotB	sugar efflux transporter
Reads supporting (aligned to +/- strand): ref base A (154/156); new base C (16/11); total (170/167)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

GTTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAA  >  NC_002947/2780951‑2781248

 ttCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAg                                                                                                                                                     >  1:6225355/1‑150 (MQ=255)

   cATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgc                                                                                                                                                   >  2:626632/1‑150 (MQ=255)

   cATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgc                                                                                                                                                   <  1:6103710/150‑1 (MQ=255)

    aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgcg                                                                                                                                                  >  2:3255597/1‑150 (MQ=255)

     ttGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCAGGCCAAGgcgc                                                                                                                                                 >  1:2274400/1‑150 (MQ=255)

       gTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTg                                                                                                                                               <  1:4090814/150‑1 (MQ=255)

       gTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCAGTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTg                                                                                                                                               <  2:2025247/150‑1 (MQ=255)

        tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                              <  2:8496469/150‑1 (MQ=255)

        tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                              <  2:6658743/150‑1 (MQ=255)

        acAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                              >  2:8328914/2‑150 (MQ=255)

            ccACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCct                                                                                                                                          <  2:8766378/150‑1 (MQ=255)

             cACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctg                                                                                                                                         <  1:3503741/150‑1 (MQ=255)

               cTTGCTGTCCTGGCTTTCGCAAAGCTGCGCCATGCTGCTGGTCAGCCGCCTCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                       <  1:6040943/150‑1 (MQ=255)

                ttGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgctg                                                                                                                                      <  2:7654241/150‑1 (MQ=255)

                ttGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgctg                                                                                                                                      <  1:6776053/150‑1 (MQ=255)

                gtGCTGTCCTGGCTTTCGCAACGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgctg                                                                                                                                      <  1:5938556/149‑1 (MQ=255)

                  gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGc                                                                                                                                    <  1:3247584/150‑1 (MQ=255)

                   cTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCCTGCTGGcc                                                                                                                                   >  1:904202/1‑150 (MQ=255)

                      tCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCcccc                                                                                                                                >  1:1425054/1‑150 (MQ=255)

                      tCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCAcc                                                                                                                                >  1:6352143/1‑150 (MQ=255)

                          ggCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGta                                                                                                                            >  1:5112241/1‑150 (MQ=255)

                             tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCCGGCCAAGGCGCGGGGCCGGCTGGCCACCCGtccaa                                                                                                                         >  1:767630/1‑146 (MQ=255)

                             tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtacta                                                                                                                         >  2:4667676/1‑150 (MQ=255)

                             tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtacta                                                                                                                         >  1:2842083/1‑150 (MQ=255)

                             tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtacaa                                                                                                                         >  1:6330791/1‑148 (MQ=255)

                             tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCACCGGtacta                                                                                                                         >  1:8159243/1‑150 (MQ=255)

                                 gCAAAGCTTCGCGATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCACGCCAAGGCCCGGGGCCTGCTGGCCACCCGTACTAaccg                                                                                                                     >  2:3866055/1‑149 (MQ=255)

                                 gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCt                                                                                                                     >  1:7314099/1‑150 (MQ=255)

                                 gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCt                                                                                                                     >  2:528818/1‑150 (MQ=255)

                                 gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCt                                                                                                                     >  2:5149324/1‑150 (MQ=255)

                                 gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCACGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCt                                                                                                                     >  2:7517689/1‑150 (MQ=255)

                                  cAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCACCGGTACTAaccgg                                                                                                                    >  1:496546/1‑148 (MQ=255)

                                  cAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCAGGGCCTGCTGGCCACCCGTACAACCCTg                                                                                                                    >  1:4891782/1‑150 (MQ=255)

                                  cAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCACTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTg                                                                                                                    <  2:4859539/150‑1 (MQ=255)

                                   aaaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAgg                                                                                                                                                    >  2:3645768/1‑117 (MQ=255)

                                   aaaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTgg                                                                                                                   >  2:4407226/1‑150 (MQ=255)

                                   aaaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTgg                                                                                                                   >  2:6783969/1‑150 (MQ=255)

                                   aaaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCCTGCTGGCCACCGGTACTACCCTgg                                                                                                                   >  2:67207/1‑150 (MQ=255)

                                       cTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCAt                                                                                                               <  1:151843/150‑1 (MQ=255)

                                        ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATg                                                                                                              >  2:4820251/1‑150 (MQ=255)

                                        ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAACGCGCTGGGCCGGCTGGCCACCCGTACTACCCAGGCCaag                                                                                                              >  1:1044080/1‑148 (MQ=255)

                                        ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCAGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATg                                                                                                              <  1:5536445/150‑1 (MQ=255)

                                           gCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTc                                                                                                           <  2:5423470/150‑1 (MQ=255)

                                             gaTGCGGCTGACCAGCAGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGGCAGCAGGCCCAGTCGCTGGGCCGGCTGTCCACCGGTAATACCCTGGCCATGGTGCt                                                                                                         >  1:5416686/2‑150 (MQ=255)

                                             cATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCt                                                                                                         <  2:7045750/150‑1 (MQ=255)

                                             cATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCAAGGTCCt                                                                                                         >  1:5670577/1‑150 (MQ=255)

                                              aTGCGGCTGACCAGCAGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTg                                                                                                        >  2:2067507/1‑150 (MQ=255)

                                                          aGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTg                                                                                            >  1:7732207/1‑150 (MQ=255)

                                                           gCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTgg                                                                                           <  1:7987510/150‑1 (MQ=255)

                                                            ccGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCCTGCTgg                                                                                                                                     <  2:2887568/107‑1 (MQ=38)

                                                              gCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGcc                                                                                        >  2:4478977/1‑150 (MQ=255)

                                                                  cGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCgtgt                                                                                    <  2:6920972/150‑1 (MQ=255)

                                                                          cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgctg                                                                                                                                      <  2:2609593/92‑1 (MQ=37)

                                                                          cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCg                                                                            >  1:3566351/1‑150 (MQ=255)

                                                                          cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCg                                                                            <  2:8606965/150‑1 (MQ=255)

                                                                           ccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                           <  2:1110960/150‑1 (MQ=255)

                                                                           ccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                           <  1:1279765/150‑1 (MQ=255)

                                                                            cTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAg                                                                          <  2:6400103/150‑1 (MQ=255)

                                                                             tGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGTCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAgg                                                                         <  2:8058248/150‑1 (MQ=255)

                                                                             tGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAgg                                                                         <  1:7018975/150‑1 (MQ=255)

                                                                                 gcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCACCCCGCCAAGGCGCAGGGCCGGCTGGCCACCCGGACTACCCAGGCCAGGGGCCGGGGCAACCCCCGGGGCCGGGTGGTGGGCGAGGCAcg                                                                     >  1:3866521/1‑149 (MQ=255)

                                                                                  cATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTg                                                                    <  1:2581119/150‑1 (MQ=255)

                                                                                    tGCGGTGTTCTGGTCCATCCCCGCCTCGCTGGCGGTGCGTGGGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTggg                                                                  <  1:2172240/150‑1 (MQ=255)

                                                                                       ggTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTg                                                               <  2:5609373/150‑1 (MQ=255)

                                                                                         tgtTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             >  2:7251711/1‑150 (MQ=255)

                                                                                         tgtTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             <  1:3662168/150‑1 (MQ=255)

                                                                                         tgtTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             <  2:3039618/150‑1 (MQ=255)

                                                                                         tgtTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             <  2:2762350/150‑1 (MQ=255)

                                                                                           ttCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGt                                                           <  1:8742282/150‑1 (MQ=255)

                                                                                              tGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCCGCTGGGCCGGCTGGCCACCCGGAATACCCAGGCCAACGGGCTGGGCAACCCGCTGGGCCGTGTGGTGGGCCAGCCACAGCGCGGGCCTATc                                                        >  1:1891399/1‑150 (MQ=255)

                                                                                              tGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGGACTACCCTGGCCATGGTCCTGGGCCTCCCGCTGGGCCGTGTGGTGGGCGAGGCCCTGGGCTGGCGTATc                                                        >  2:5548006/1‑150 (MQ=255)

                                                                                                     tCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCCGGCTGGCCACCGGTACTACCCGGGCCAAGGTCCGGGGCATCCCGCGGGGCCGGGTGGTGGGCGAGGCACTGGGCGGGCGGATCCCGGTcc                                                 >  1:8308213/1‑150 (MQ=255)

                                                                                                        ccGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                              <  1:6681425/150‑1 (MQ=255)

                                                                                                        ccGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                              <  2:3374961/150‑1 (MQ=255)

                                                                                                        ccGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                              <  2:135915/150‑1 (MQ=255)

                                                                                                                cTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGcc                                      >  2:4356127/1‑150 (MQ=255)

                                                                                                                 tGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCg                                     <  2:6308772/150‑1 (MQ=255)

                                                                                                                  ggCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgg                                    <  2:4630232/150‑1 (MQ=255)

                                                                                                                  ggCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgg                                    <  2:244055/150‑1 (MQ=255)

                                                                                                                   gCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCggg                                   <  1:8659283/150‑1 (MQ=255)

                                                                                                                   gCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACGACCCTGGCCAGGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCggg                                   <  2:6984487/150‑1 (MQ=255)

                                                                                                                       tGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAg                               <  1:3395328/150‑1 (MQ=255)

                                                                                                                        gCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgc                              <  1:1122297/150‑1 (MQ=255)

                                                                                                                            gtgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTa                          <  2:6095994/150‑1 (MQ=255)

                                                                                                                            gtgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCCCCGGTACTACCCAGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGGGGTGGGCGAGGCACTGGGCGGGCGTATCACGTTCCTGAGCCTTGCCGGGGTAGCGCTa                          >  2:2283285/1‑150 (MQ=255)

                                                                                                                               gcgcCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCa                                                       >  2:1290462/1‑118 (MQ=255)

                                                                                                                               gcgcCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGcc                       <  1:8470194/150‑1 (MQ=255)

                                                                                                                                  ccgccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacc                    <  2:8679942/150‑1 (MQ=255)

                                                                                                                                  ccgccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacc                    <  2:7130817/150‑1 (MQ=255)

                                                                                                                                    cggtggCCAGCAGGCCCAGCGCCTTGGCCTGCTGGCCCCCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAg                                                                          >  1:6492017/5‑94 (MQ=21)

                                                                                                                                         gCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCttt             <  1:578367/150‑1 (MQ=255)

                                                                                                                                           cagcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc           >  2:1763895/1‑150 (MQ=255)

                                                                                                                                            agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc          >  2:3514197/1‑150 (MQ=255)

                                                                                                                                            agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc          >  2:1627922/1‑150 (MQ=255)

                                                                                                                                             gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCt         <  1:5187949/150‑1 (MQ=255)

                                                                                                                                              cagGCAAAGGCGCTGGGCCTGCTGGCAACCGGTACTACCCTGGCCATGGTCCTGGTCAACCCGCTCGGCCGTGTGGTGGGCGAGGCACTGGGCTGTCGTATCACGTTCCTTAGAATTGACGGGGTAGCGCTAGCAACCATGCTTGGCATg        >  1:6618293/1‑150 (MQ=255)

                                                                                                                                               agGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGTGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGa       <  1:1955198/150‑1 (MQ=255)

                                                                                                                                               agGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGa       >  1:2761023/1‑150 (MQ=255)

                                                                                                                                                ggCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt      <  1:4372098/150‑1 (MQ=255)

                                                                                                                                                   cAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCaa   >  2:3457019/1‑150 (MQ=255)

GTTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAA  >  NC_002947/2780951‑2781248

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