BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,739,214	T→C	9.4%	E70E (GAA→GAG) ‡	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,739,214	0	T	C	9.4%	464.9 / 27.3	159	E70E (GAA→GAG) ‡	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base T (89/55); new base C (9/6); total (98/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 160 total aligned reads displayed.

GTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGT  >  NC_002947/5739067‑5739363

gTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTc                                                                                                                                                     >  2:1178158/1‑150 (MQ=255)

gTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTc                                                                                                                                                     >  2:3597625/1‑150 (MQ=255)

 ttGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcg                                                                                                                                                    >  1:2100962/1‑150 (MQ=255)

 ttGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcg                                                                                                                                                    >  1:1913791/1‑150 (MQ=255)

     aGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcgacga                                                                                                                                                <  1:3537798/150‑1 (MQ=255)

       tgtTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATc                                                                                                                                              >  1:3388857/1‑150 (MQ=255)

         ttCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCtt                                                                                                                                            >  1:4753445/1‑150 (MQ=255)

           cccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGt                                                                                                                                          >  1:3243662/1‑150 (MQ=255)

           cccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGt                                                                                                                                          >  1:817006/1‑150 (MQ=255)

            ccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTc                                                                                                                                         <  2:828071/150‑1 (MQ=255)

            ccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTc                                                                                                                                         >  2:1811357/1‑150 (MQ=255)

             cAGTTGGCACAGTGTGGCACCAGGTTGGAGTTGGCAAACAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCg                                                                                                                                        >  2:4519170/1‑150 (MQ=255)

               gTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTa                                                                                                                                      <  1:502159/150‑1 (MQ=255)

                 tGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGc                                                                                                                                    <  1:4648891/150‑1 (MQ=255)

                         tgtgGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAg                                                                                                                            <  2:5111449/150‑1 (MQ=255)

                              cACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgcg                                                                                                                       >  1:1987684/1‑150 (MQ=255)

                               aCCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgcgc                                                                                                                      >  2:1074121/1‑150 (MQ=255)

                                ccAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgcgcg                                                                                                                     >  2:3989034/1‑150 (MQ=255)

                                 cAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgcgcgc                                                                                                                    >  2:2725404/1‑150 (MQ=255)

                                  aGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGcgcgcg                                                                                                                   <  2:4051533/150‑1 (MQ=255)

                                      tGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCGCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCgcgcggcgc                                                                                                               <  1:4900597/150‑1 (MQ=255)

                                          gTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGcc                                                                                                           <  2:4252836/150‑1 (MQ=255)

                                          gTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTAGTCGTAGCAGGCCCAGTTGCGCgcggcgcggca                                                                                                           >  2:3073852/1‑149 (MQ=255)

                                           ttGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCg                                                                                                          >  2:2348561/1‑150 (MQ=255)

                                               cAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGCCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATAc                                                                                                      >  1:3763281/1‑150 (MQ=255)

                                                gatgagaaCGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGcc                                                                                                     >  1:4398791/8‑150 (MQ=255)

                                                  tCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATAccgc                                                                                                   >  1:496919/1‑150 (MQ=255)

                                                   gagaaCGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATgccgcc                                                                                                  >  2:2408555/5‑150 (MQ=255)

                                                   cAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATAccgcc                                                                                                  <  1:4525609/150‑1 (MQ=255)

                                                       aCGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGtat                                                                                              >  1:4763714/1‑150 (MQ=255)

                                                                   tccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGctca                                                                                  <  2:3550851/146‑2 (MQ=255)

                                                                   tccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGctca                                                                                  >  2:3778979/5‑149 (MQ=255)

                                                                   tccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGctca                                                                                  <  2:4193489/146‑2 (MQ=255)

                                                                      cTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTc                                                                               <  1:3413091/150‑1 (MQ=255)

                                                                       tgtgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCg                                                                              <  2:1750856/150‑1 (MQ=255)

                                                                         tgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCggg                                                                            <  2:1913791/150‑1 (MQ=255)

                                                                         tgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCggg                                                                            >  1:1060647/1‑150 (MQ=255)

                                                                         tgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCggg                                                                            >  1:1060503/1‑150 (MQ=255)

                                                                           ggTCTTGAACACACCGACTGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCCGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGtt                                                                          >  2:3171790/1‑150 (MQ=255)

                                                                            gTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGttt                                                                         >  1:3716968/1‑150 (MQ=255)

                                                                             tCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGtttt                                                                        >  1:2830987/1‑150 (MQ=255)

                                                                               ttGAACACCCCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTcc                                                                      <  2:4944070/150‑1 (MQ=255)

                                                                               ttGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATAccgcc                                                                                                  >  1:4780134/1‑122 (MQ=255)

                                                                                tGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCcg                                                                     >  1:3114632/1‑150 (MQ=255)

                                                                                tGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCcg                                                                     <  2:2027558/150‑1 (MQ=255)

                                                                                 gAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCcgc                                                                    <  1:524060/150‑1 (MQ=255)

                                                                                     acacCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGAcc                                                                >  1:4717778/1‑150 (MQ=255)

                                                                                      cacCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCt                                                               <  2:3779392/150‑1 (MQ=255)

                                                                                       gcccaccgGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGCTCATGCGGGCGTTCGGCGAcatc                                                              >  1:3981306/5‑147 (MQ=255)

                                                                                       accgaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCtg                                                              >  2:2905181/1‑150 (MQ=255)

                                                                                        ccgaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCtgt                                                             <  2:2997311/150‑1 (MQ=255)

                                                                                         cgaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCtgtg                                                            >  1:3749360/1‑150 (MQ=255)

                                                                                         ccaccgGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGCTCATGCGGGCGTTCGGCGAcatccg                                                            <  2:1750502/148‑6 (MQ=255)

                                                                                          gaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTgg                                                           >  2:2659311/1‑150 (MQ=255)

                                                                                          gaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTgg                                                           >  1:4059926/1‑150 (MQ=255)

                                                                                          gaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTgg                                                           >  1:3330734/1‑150 (MQ=255)

                                                                                          caccgGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGCTCATGCGGGCGTTCGGCGAcatccgg                                                           <  2:1465455/149‑7 (MQ=255)

                                                                                               gCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGa                                                      <  1:418129/150‑1 (MQ=255)

                                                                                               gCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGa                                                      <  1:4146605/150‑1 (MQ=255)

                                                                                                cTTGCCCGACTGCCCCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAg                                                     <  1:64648/150‑1 (MQ=255)

                                                                                                cTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAg                                                     <  1:4355757/150‑1 (MQ=255)

                                                                                                cTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAg                                                     <  1:2348561/150‑1 (MQ=255)

                                                                                                cTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAg                                                     <  2:3491703/150‑1 (MQ=255)

                                                                                                cTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAg                                                     <  1:1178158/150‑1 (MQ=255)

                                                                                                 ttGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcga                                                                                                                                                   <  2:525658/55‑1 (MQ=255)

                                                                                                 ttGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcga                                                                                                                                                   >  1:525658/1‑55 (MQ=255)

                                                                                                  tGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgt                                                   <  1:5063158/150‑1 (MQ=255)

                                                                                                   gCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgtt                                                  <  1:3161486/150‑1 (MQ=255)

                                                                                                   gCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgtt                                                  <  2:709054/150‑1 (MQ=255)

                                                                                                    cccGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgttg                                                 <  2:215467/150‑1 (MQ=255)

                                                                                                    cccGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgttg                                                 >  1:1588068/1‑150 (MQ=255)

                                                                                                    cccGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgttg                                                 <  2:3243662/150‑1 (MQ=255)

                                                                                                     ccGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCCGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgttgt                                                <  2:3082961/150‑1 (MQ=255)

                                                                                                     ccGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgttgt                                                <  1:2995928/150‑1 (MQ=255)

                                                                                                     ccGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgttgt                                                >  1:1482175/1‑150 (MQ=255)

                                                                                                          tGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTtcatc                                           >  1:2070238/1‑150 (MQ=255)

                                                                                                           gCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTtcatca                                          >  1:3918649/1‑150 (MQ=255)

                                                                                                             aCCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGc                                        <  1:2665996/150‑1 (MQ=255)

                                                                                                             aCCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGc                                        >  1:1074269/1‑150 (MQ=255)

                                                                                                             aCCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGc                                        >  1:1074313/1‑150 (MQ=255)

                                                                                                             aCCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGc                                        <  2:249170/150‑1 (MQ=255)

                                                                                                              ccAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCa                                       <  2:2806614/150‑1 (MQ=255)

                                                                                                                agcagGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATg                                     >  1:4773914/1‑150 (MQ=255)

                                                                                                                 gcagGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTgg                                                           <  1:777185/127‑1 (MQ=255)

                                                                                                                 gcagGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTgg                                                           >  2:777185/1‑127 (MQ=255)

                                                                                                                       tttCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTg                              >  1:4996211/1‑150 (MQ=255)

                                                                                                                            tcgtcTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCt                         >  2:1562782/1‑150 (MQ=255)

                                                                                                                              gtcCTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTc                       >  1:1793058/1‑150 (MQ=255)

                                                                                                                                cTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCgg                     >  1:711350/1‑150 (MQ=255)

                                                                                                                                 ttCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGt                    <  1:3597625/150‑1 (MQ=255)

                                                                                                                                    cAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTcag                 >  2:2238202/1‑150 (MQ=255)

                                                                                                                                           gTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGATCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCtt          >  1:67290/1‑150 (MQ=255)

                                                                                                                                              aGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTg       >  2:1645912/1‑150 (MQ=255)

                                                                                                                                               gggTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTgg      >  1:5067454/1‑150 (MQ=255)

                                                                                                                                                ggTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGc     >  1:1466750/1‑150 (MQ=255)

                                                                                                                                                ggTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTGCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGc     >  2:3837057/1‑150 (MQ=255)

                                                                                                                                                 gTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCg    <  1:2774469/150‑1 (MQ=255)

                                                                                                                                                 gTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCg    >  2:4355429/1‑150 (MQ=255)

                                                                                                                                                 gTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCg    >  2:3346852/1‑150 (MQ=255)

                                                                                                                                                   ttCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGt  <  1:1728665/150‑1 (MQ=255)

GTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGT  >  NC_002947/5739067‑5739363

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