BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,739,076	T→A	19.0%	E116D (GAA→GAT)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,739,076	0	T	A	19.0%	428.2 / 96.9	190	E116D (GAA→GAT)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base T (73/80); new base A (20/16); total (93/97)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

Only 100 of 194 total aligned reads displayed.

TCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTT  >  NC_002947/5738928‑5739225

tCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGtt                                                                                                                                                      >  2:3477750/1‑150 (MQ=255)

 gaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATc                                                                                                                                                     >  2:5300102/2‑150 (MQ=255)

  aGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTcc                                                                                                                                                    <  2:2090542/150‑1 (MQ=255)

    gCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCa                                                                                                                                                  <  1:904652/150‑1 (MQ=255)

    gCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCa                                                                                                                                                  >  2:2500683/1‑150 (MQ=255)

             gTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACCTCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCa                                                                                                                                         <  2:3014643/150‑1 (MQ=255)

              taatGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCaa                                                                                                                                        >  2:4241346/4‑149 (MQ=255)

                 tGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAgtgt                                                                                                                                     <  2:309770/150‑1 (MQ=255)

                 tGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAatgc                                                                                                                                     >  2:3504344/1‑149 (MQ=255)

                     ggCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCa                                                                                                                                 >  1:979187/1‑150 (MQ=255)

                     ggCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCa                                                                                                                                 <  1:1620008/150‑1 (MQ=255)

                       cGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGCAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTGCCCAGTTGGCCCAGTGTGGCAcc                                                                                                                               <  1:137770/150‑1 (MQ=255)

                        gACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCa                                                                                                                              >  1:1058861/1‑150 (MQ=255)

                            ggCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGtt                                                                                                                          <  2:4098003/150‑1 (MQ=255)

                             gCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCATCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTg                                                                                                                         <  1:4364817/150‑1 (MQ=255)

                             gCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTg                                                                                                                         >  2:1270436/1‑150 (MQ=255)

                             gCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTg                                                                                                                         <  1:4364886/150‑1 (MQ=255)

                               ttCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGa                                                                                                                       <  1:171685/150‑1 (MQ=255)

                               ttCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGa                                                                                                                       <  1:3932345/150‑1 (MQ=255)

                                  gacgaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGtt                                                                                                                    >  2:1896529/1‑150 (MQ=255)

                                  gacgaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGtt                                                                                                                    <  1:4210459/150‑1 (MQ=255)

                                  cacgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAAtt                                                                                                                    <  1:3634626/149‑1 (MQ=255)

                                   acgaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTg                                                                                                                   <  1:2845741/150‑1 (MQ=255)

                                   acgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTg                                                                                                                   <  2:5323374/150‑1 (MQ=255)

                                   acgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTg                                                                                                                   >  2:2680646/1‑150 (MQ=255)

                                    cgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTgg                                                                                                                  <  1:1170298/150‑1 (MQ=255)

                                    cgaAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGACATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTgg                                                                                                                  >  1:648506/1‑150 (MQ=255)

                                     gaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGc                                                                                                                 >  1:3492593/1‑150 (MQ=255)

                                     gaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGc                                                                                                                 >  1:2811767/1‑150 (MQ=255)

                                       aGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCaa                                                                                                               <  2:1711976/150‑1 (MQ=255)

                                         gTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATc                                                                                                             <  1:1270436/150‑1 (MQ=255)

                                          tttCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCa                                                                                                            <  2:4573303/150‑1 (MQ=255)

                                                aGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCAcgc                                                                                                      <  1:2657113/150‑1 (MQ=255)

                                                   tGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCAcgcgcg                                                                                                   >  2:4860935/1‑150 (MQ=255)

                                                   tGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAAcgcgcg                                                                                                   <  2:4173301/150‑1 (MQ=255)

                                                    gCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCa                                                                                                         >  2:2253135/1‑143 (MQ=255)

                                                    gCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCa                                                                                                         <  1:2253135/143‑1 (MQ=255)

                                                    gCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCa                                                                                                         >  2:4032880/1‑143 (MQ=255)

                                                    gCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCa                                                                                                         <  1:4032880/143‑1 (MQ=255)

                                                     cccTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGt                                                                                                 <  2:2173986/150‑1 (MQ=255)

                                                     cccTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGt                                                                                                 <  2:5189510/150‑1 (MQ=255)

                                                       cTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGcc                                                                                                                                           >  2:702416/1‑106 (MQ=255)

                                                            cGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGc                                                                                          <  1:639685/150‑1 (MQ=255)

                                                             gATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCt                                                                                         >  1:809932/1‑150 (MQ=255)

                                                                gCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCtgtg                                                                                      <  1:4078441/150‑1 (MQ=255)

                                                                 cccTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTgg                                                                                     <  1:5171402/150‑1 (MQ=255)

                                                                    tgtgcaccGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGtg                                                                                  >  1:732237/7‑149 (MQ=255)

                                                                    tgtgcaccGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGtg                                                                                  >  1:3466391/7‑149 (MQ=255)

                                                                      gCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCtt                                                                                <  2:1841496/150‑1 (MQ=255)

                                                                        tGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGa                                                                              >  1:3386475/1‑150 (MQ=255)

                                                                         gCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGaa                                                                             <  1:895972/150‑1 (MQ=255)

                                                                         accGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGtgcggaa                                                                             >  2:5037294/2‑144 (MQ=255)

                                                                         accGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGtgcggaa                                                                             >  2:2343858/2‑144 (MQ=255)

                                                                            gATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAacac                                                                          >  2:2445859/1‑150 (MQ=255)

                                                                              tGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAacacac                                                                        >  2:4255213/1‑150 (MQ=255)

                                                                                tAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccg                                                                      >  2:46970/1‑150 (MQ=255)

                                                                                  gATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgac                                                                    >  2:299327/1‑150 (MQ=255)

                                                                                  gATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgac                                                                    <  1:3865902/150‑1 (MQ=255)

                                                                                  gATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgac                                                                    >  1:4911098/1‑150 (MQ=255)

                                                                                   aTCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgacc                                                                   >  1:2890039/1‑150 (MQ=255)

                                                                                    tCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCGTGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgaccg                                                                  <  2:757615/150‑1 (MQ=255)

                                                                                       aGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCt                                                               >  1:3747166/1‑150 (MQ=255)

                                                                                       aGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCt                                                               >  1:1731115/1‑150 (MQ=255)

                                                                                       aGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCt                                                               >  2:3780449/1‑150 (MQ=255)

                                                                                        gCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCtt                                                              <  2:897714/150‑1 (MQ=255)

                                                                                            ccGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGccc                                                          <  2:2811767/150‑1 (MQ=255)

                                                                                              ggcggTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGa                                                        >  2:1353906/1‑150 (MQ=255)

                                                                                              ggcggTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGGTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAAGCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGa                                                        >  2:3795185/1‑150 (MQ=255)

                                                                                                   tCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCa                                                   <  1:4860935/150‑1 (MQ=255)

                                                                                                    ggccTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCAc                                                  >  1:3831481/4‑150 (MQ=255)

                                                                                                      tCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTggg                                                                                    <  1:2354406/114‑1 (MQ=255)

                                                                                                        tGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACcagc                                              >  2:2749632/1‑150 (MQ=255)

                                                                                                           ccATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                           <  2:4303909/150‑1 (MQ=255)

                                                                                                           ccATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                           <  2:4188475/150‑1 (MQ=255)

                                                                                                            cATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGAGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAggg                                          >  2:984874/1‑150 (MQ=255)

                                                                                                             aTACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGt                                         >  1:2037452/1‑150 (MQ=255)

                                                                                                              tACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTc                                                                                  <  2:3861139/108‑1 (MQ=255)

                                                                                                               aCATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGttt                                       >  2:382347/1‑150 (MQ=255)

                                                                                                                 aTCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcg                                     >  2:2368139/1‑150 (MQ=255)

                                                                                                                  tCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGTTCACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgt                                    >  2:2060353/1‑150 (MQ=255)

                                                                                                                  tCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                           <  2:3087799/143‑1 (MQ=255)

                                                                                                                  tCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                           >  1:3087736/1‑143 (MQ=255)

                                                                                                                  tCGCCAGGCCCTTGGCGACCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcga                                    >  2:3309363/1‑149 (MQ=255)

                                                                                                                   cGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgtc                                   <  1:4275491/150‑1 (MQ=255)

                                                                                                                   cGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgtc                                   <  1:5031752/150‑1 (MQ=255)

                                                                                                                     ccAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcg                                     >  2:2963993/1‑146 (MQ=255)

                                                                                                                      cAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCgtcgtc                                <  2:2890039/150‑1 (MQ=255)

                                                                                                                       aGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCt                               >  1:5286535/1‑150 (MQ=255)

                                                                                                                                  gTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGt                    <  1:140532/150‑1 (MQ=255)

                                                                                                                                    ccAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAacaca                                                                         <  1:1899852/95‑1 (MQ=255)

                                                                                                                                     cAGTTCGTTGACGTGGTCCCAGTTGGCCCAGTGTGGCACCAGGCTGGAGTTGGCAATCAGCACGCGCGGTGCGCTGCTGTGGGTCTTGAACCCACCGACCGGCTTGCCCGACGGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAg                 <  2:5394533/150‑1 (MQ=255)

                                                                                                                                      aGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAgg                <  1:1573876/150‑1 (MQ=255)

                                                                                                                                       gTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAggg               >  2:1145412/1‑150 (MQ=255)

                                                                                                                                        ttCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGt              >  1:5359378/1‑150 (MQ=255)

                                                                                                                                             tGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcga         <  1:3115074/150‑1 (MQ=255)

                                                                                                                                             tGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcga         >  1:4958575/1‑150 (MQ=255)

                                                                                                                                              gAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcgac        <  1:1010360/150‑1 (MQ=255)

                                                                                                                                               aaGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcgacg       >  2:4489980/1‑150 (MQ=255)

                                                                                                                                                aGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                           >  2:2576207/1‑113 (MQ=255)

                                                                                                                                                aGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                           <  1:2576207/113‑1 (MQ=255)

TCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTT  >  NC_002947/5738928‑5739225

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