BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,739,220	G→A	6.5%	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	6.5%	1139.7 / 44.4	367	I68I (ATC→ATT)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base G (202/141); new base A (15/9); total (217/150)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 386 total aligned reads displayed.

GTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAGCTT  >  NC_002947/5739073‑5739369

gtgtTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGAt                                                                                                                                                     >  2:6427351/1‑150 (MQ=255)

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

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

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

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

    tCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTg                                                                                                                                                 >  1:794098/1‑150 (MQ=255)

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

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

       cAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCg                                                                                                                                              <  2:2414238/150‑1 (MQ=255)

       cAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCg                                                                                                                                              >  1:2457919/1‑150 (MQ=255)

         gTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTa                                                                                                                                            >  1:4695669/1‑150 (MQ=255)

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

               ccAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACTGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGc                                                                                                                                      <  2:4002670/150‑1 (MQ=255)

                cAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGcc                                                                                                                                     <  2:8672557/150‑1 (MQ=255)

                 aGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGccc                                                                                                                                    <  2:1375345/150‑1 (MQ=255)

                  gtgtGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCa                                                                                                                                   <  1:5230450/150‑1 (MQ=255)

                    gtgGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGt                                                                                                                                 >  1:310328/1‑150 (MQ=255)

                       gCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgc                                                                                                                              >  2:1168644/1‑150 (MQ=255)

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

                              gTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCgcgcggc                                                                                                                       >  1:1331510/1‑150 (MQ=255)

                              gTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCgcgcggc                                                                                                                       >  2:8397183/1‑150 (MQ=255)

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

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

                                        gCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATa                                                                                                             >  1:697356/1‑150 (MQ=255)

                                          aaTCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATAcc                                                                                                           >  1:4834120/1‑150 (MQ=255)

                                           aTCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATAccg                                                                                                          <  2:5447980/150‑1 (MQ=255)

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

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

                                                 aCGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGACCCAGTTGCGCGCGGCGCGGCCGATACCGCCGtat                                                                                                    <  2:1753748/150‑1 (MQ=255)

                                                   gcgcgGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATac                                                                                                  <  1:1170062/150‑1 (MQ=255)

                                                   gcgcgGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATac                                                                                                  <  1:1653951/150‑1 (MQ=255)

                                                       ggTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATaccacc                                                                                              <  2:3013242/150‑1 (MQ=255)

                                                        gTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATAccacca                                                                                             <  1:2324482/150‑1 (MQ=255)

                                                           catccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGct                                                                                          >  2:5157741/7‑150 (MQ=255)

                                                             ttGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGctct                                                                                        >  2:2029339/1‑150 (MQ=255)

                                                              ccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGctcat                                                                                       <  2:2964733/147‑3 (MQ=255)

                                                                cTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTc                                                                                     >  1:3007672/1‑150 (MQ=255)

                                                                cTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTc                                                                                     >  1:990807/1‑150 (MQ=255)

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

                                                                  gtgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCgg                                                                                   >  1:3374370/1‑150 (MQ=255)

                                                                  gtgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCgg                                                                                   <  2:1874383/150‑1 (MQ=255)

                                                                  gtgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCgg                                                                                   <  1:197474/150‑1 (MQ=255)

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

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

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

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

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

                                                                           gAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCcgc                                                                          >  2:8098342/1‑150 (MQ=255)

                                                                           gAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCcgc                                                                          >  1:1592019/1‑150 (MQ=255)

                                                                             acacgcCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGCTCATGCGGGCGTTCGGCGa                                                                        <  1:430571/150‑1 (MQ=255)

                                                                             acacacCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGa                                                                        >  2:5665415/1‑150 (MQ=255)

                                                                             acacacCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGa                                                                        >  1:8456245/1‑150 (MQ=255)

                                                                              cacacCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGAc                                                                       <  2:5253025/150‑1 (MQ=255)

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

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

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

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

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

                                                                                   cgaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGt                                                                                                      >  1:2924023/1‑114 (MQ=255)

                                                                                    caccgGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGCTCATGCGGGCGTTCGGCGAcatccgg                                                                 <  1:7560403/149‑7 (MQ=255)

                                                                                         gCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGAc                                                                       >  1:1077899/1‑139 (MQ=255)

                                                                                         gCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGa                                                            >  1:5522650/1‑150 (MQ=255)

                                                                                         gCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGa                                                            >  1:1441616/1‑150 (MQ=255)

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

                                                                                           ttGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAgg                                                          >  1:8257225/1‑150 (MQ=255)

                                                                                           ttGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAgg                                                          >  1:6808195/1‑150 (MQ=255)

                                                                                            tGCCCGACTGCACCAGCTGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGTCGATCTTGTCGTAGCAGTCCCAGTTGCGCGCTGCGCGGGCTTTACCGCCGTATACCACCAGCTCTTTCTGGTTTTCCGCGACCTGTGGATCGAGgt                                                         <  1:1295232/150‑1 (MQ=255)

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

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

                                                                                               ccGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTAAACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGgttgt                                                      >  2:7572083/1‑150 (MQ=255)

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

                                                                                                   cTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTtcat                                                  <  2:307936/150‑1 (MQ=255)

                                                                                                   cTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTtcat                                                  <  2:4575982/150‑1 (MQ=255)

                                                                                                    tGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTtcatc                                                 <  1:7842037/150‑1 (MQ=255)

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

                                                                                                            cagGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATgcg                                         >  2:7991531/1‑150 (MQ=255)

                                                                                                                     gtcgtcTTCCGGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGgcgc                                >  1:5408559/1‑150 (MQ=255)

                                                                                                                       cgtcTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCtt                              <  1:3638758/150‑1 (MQ=255)

                                                                                                                        gtcTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTc                             >  1:8142395/1‑150 (MQ=255)

                                                                                                                        gtcTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTc                             <  2:77712/150‑1 (MQ=255)

                                                                                                                        gtcTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTc                             >  1:6335923/1‑150 (MQ=255)

                                                                                                                         tcTTCCAGGCGGGTCAGGGTTTCGACGATCTTTTCGTAGAAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGAGCAGTGGCGCTTCg                            >  1:3477098/1‑150 (MQ=255)

                                                                                                                         tcTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGGGGTTGTTCATCAGCATGCGCAGTGGCGCTTCg                            <  2:2697894/150‑1 (MQ=255)

                                                                                                                         tcTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCg                            >  1:7276383/1‑150 (MQ=255)

                                                                                                                         tcTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCg                            >  2:6967409/1‑150 (MQ=255)

                                                                                                                           ttCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGt                          >  1:6720386/1‑150 (MQ=255)

                                                                                                                             ccAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTca                        <  2:3776200/150‑1 (MQ=255)

                                                                                                                                  cGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTcagccag                   <  1:8093987/150‑1 (MQ=255)

                                                                                                                                       cAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCtttt              <  1:4658415/150‑1 (MQ=255)

                                                                                                                                        aGGGTATCGACGATCATGTCTTAGCAGGCCCAGATGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCAGAGGATCGAGGATGATCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTATg             >  2:6871502/1‑150 (MQ=255)

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

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

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

                                                                                                                                              tCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTc       <  2:6158326/150‑1 (MQ=255)

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

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

                                                                                                                                                gacgaTCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTCGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAg     >  2:2191450/1‑150 (MQ=255)

                                                                                                                                                 acgaTCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAGc    >  1:1752803/1‑150 (MQ=255)

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

                                                                                                                                                   gaTCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGAGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAGCtt  >  1:7705834/1‑150 (MQ=255)

GTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTGTGGATCGAGGTTGTTCATCAGCATGCGCAGTGGCGCTTCGGTCAGCCAGCTTTTGGCGGTCAGCTT  >  NC_002947/5739073‑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: `‑`