BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,739,220	G→A	9.9%	I68I (ATC→ATT)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,739,220	0	G	A	9.9%	426.5 / 28.6	151	I68I (ATC→ATT)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base G (82/54); new base A (9/6); total (91/60)
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.64e-01

Only 100 of 163 total aligned reads displayed.

AGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAGCTT  >  NC_002947/5739072‑5739369

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

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

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

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

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

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

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

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

              gCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAg                                                                                                                                        >  1:1896356/1‑150 (MQ=255)

               cccAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTAGTAGCAgg                                                                                                                                       >  2:4171645/1‑150 (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)

                             aGGTTGGTGTTGGCAATCAGCTCGCGCGGTGCGTTGCTGTGGGTCTTGAACTCACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGcgcgcg                                                                                                                         <  2:4051489/150‑1 (MQ=255)

                                ttGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCgcgcggcg                                                                                                                      >  1:1523305/1‑150 (MQ=255)

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

                                     gTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGTGCGGcc                                                                                                                 <  1:3989034/150‑1 (MQ=255)

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

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

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

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

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

                                                gCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGt                                                                                                      >  1:1748174/1‑150 (MQ=255)

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

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

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

                                                           gcatccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGc                                                                                           <  1:516679/143‑1 (MQ=255)

                                                            cGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGACCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGct                                                                                          >  1:3624655/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)

                                                                 ttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGctcatgc                                                                                     >  2:4343873/2‑146 (MQ=255)

                                                                 cTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTc                                                                                     >  1:2806614/1‑150 (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)

                                                                        tCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCgggctca                                                                              >  2:2274915/1‑146 (MQ=255)

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

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

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

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

                                                                          ttGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTcc                                                                            <  2:3058406/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                    ccaccgGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGCTCATGCGGGCGTTCGGCGAcatccg                                                                  <  1:4997090/148‑6 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                           cTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAg                                                           <  1:1178158/150‑1 (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)

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

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

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

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

                                                                                                  gACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTtc                                                    <  1:3423730/150‑1 (MQ=255)

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

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

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

                                                                                                       cACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAg                                               >  1:3395126/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                   ttCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTgg                                   >  2:2395237/1‑150 (MQ=255)

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

                                                                                                                          tcTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCg                            <  1:722204/150‑1 (MQ=255)

                                                                                                                            ttCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGACGTATACCACCAGCTCTTTCGGGTTTTCAGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGt                          >  1:4926404/1‑150 (MQ=255)

                                                                                                                                      gTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCtt                >  1:67632/1‑150 (MQ=255)

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

                                                                                                                                       tCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCttt               <  1:53718/150‑1 (MQ=255)

                                                                                                                                        cAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCtttt              >  2:3616302/1‑150 (MQ=255)

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

                                                                                                                                          gggTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTgg            <  1:4171645/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                                 gacgaTCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAg     >  2:3882067/1‑150 (MQ=255)

                                                                                                                                                 gacgaTCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAg     >  2:1380273/1‑150 (MQ=255)

                                                                                                                                                 gacgaTCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAg     >  1:3351692/1‑150 (MQ=255)

                                                                                                                                                    gaTCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAGCtt  >  1:2338442/1‑150 (MQ=255)

AGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAGCTT  >  NC_002947/5739072‑5739369

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