BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,738,862	G→T	5.7%	L188M (CTG→ATG)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,738,862	0	G	T	5.7%	591.2 / 17.6	176	L188M (CTG→ATG)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base G (79/87); new base T (5/5); total (84/92)
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.76e-01

Only 100 of 188 total aligned reads displayed.

GGAGATGGCCTTGCCTTCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAG  >  NC_002947/5738713‑5739010

ggAGATGGCCTTGCCTTCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCccagccag                                                                                                                                                      >  2:2359987/1‑150 (MQ=255)

 gagaTGGCCTTGCCTTCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTTGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAgg                                                                                                                                                     <  1:257949/150‑1 (MQ=255)

  agaTGGCCTTGCCTTCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAggg                                                                                                                                                    >  2:19551/1‑150 (MQ=255)

      ggCCTTGCCTTCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTcgc                                                                                                                                                >  1:1361872/1‑150 (MQ=255)

        ccttgccttCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCgg                                                                                                                                              <  2:357117/150‑1 (MQ=255)

         cttgccttCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGc                                                                                                                                             >  2:2685843/1‑150 (MQ=255)

          ttgccttCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGcc                                                                                                                                            <  2:978776/150‑1 (MQ=255)

             ccttCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGACATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGt                                                                                                                                         <  1:2079406/150‑1 (MQ=255)

              cttCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTg                                                                                                                                        <  2:3798429/150‑1 (MQ=255)

              cttCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTg                                                                                                                                        <  2:2894098/150‑1 (MQ=255)

               ttCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTgg                                                                                                                                       >  1:510673/1‑150 (MQ=255)

                  ggcggTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTg                                                                                                                                    <  2:1391958/150‑1 (MQ=255)

                   gcggTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTgg                                                                                                                                   <  1:2165696/150‑1 (MQ=255)

                     ggTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGgc                                                                                                                                 <  1:822246/150‑1 (MQ=255)

                     ggTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGgc                                                                                                                                 <  2:4737312/150‑1 (MQ=255)

                     ggTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGgc                                                                                                                                 <  2:1235487/150‑1 (MQ=255)

                            ttGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGccc                                                                                                                          >  2:1438845/1‑150 (MQ=255)

                            ttGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGccc                                                                                                                          >  2:1942232/1‑150 (MQ=255)

                               gCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGTAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATg                                                                                                                       <  1:4871650/150‑1 (MQ=255)

                               gCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATg                                                                                                                       <  1:4871705/150‑1 (MQ=255)

                               gCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATg                                                                                                                       <  1:1592526/150‑1 (MQ=255)

                                 gATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATgcc                                                                                                                     <  1:5072048/150‑1 (MQ=255)

                                 gATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATgcc                                                                                                                     >  2:4157929/1‑150 (MQ=255)

                                       tgcgAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAg                                                                                                               <  2:4923483/150‑1 (MQ=255)

                                        gcgAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAgg                                                                                                              <  1:2784266/150‑1 (MQ=255)

                                         cgAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGc                                                                                                             <  2:936734/150‑1 (MQ=255)

                                         cgAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGc                                                                                                             <  1:4684539/150‑1 (MQ=255)

                                                gtcgtcGAGGTCAGTGGCCTGCTCGTCGACGTAGGGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGTTGTTCAGCGAGCATGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCTGCGGt                                                                                                      <  2:367725/150‑1 (MQ=255)

                                                  cgtcgAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGa                                                                                                    <  1:3557617/150‑1 (MQ=255)

                                                  cgtcgAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGa                                                                                                    <  1:1001820/150‑1 (MQ=255)

                                                   gtcgAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAg                                                                                                   >  2:4029690/1‑150 (MQ=255)

                                                    tcgAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGt                                                                                                  >  2:4065613/1‑150 (MQ=255)

                                                    tcgAGGTAAGTGGCCTGCTCGTCGAAGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGt                                                                                                  >  2:4588285/1‑150 (MQ=255)

                                                     cgAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTa                                                                                                 <  2:2871894/150‑1 (MQ=255)

                                                         gTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCa                                                                                             <  1:2117374/150‑1 (MQ=255)

                                                               ggCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGcc                                                                                       <  2:5007034/150‑1 (MQ=255)

                                                                 ccTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCtt                                                                                     >  2:717014/1‑150 (MQ=255)

                                                                   tGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTc                                                                                   <  2:3246185/150‑1 (MQ=255)

                                                                   tGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTc                                                                                   >  1:1719107/1‑150 (MQ=255)

                                                                    gCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCa                                                                                  <  2:176430/150‑1 (MQ=255)

                                                                    gCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCa                                                                                  <  1:5028511/150‑1 (MQ=255)

                                                                    gCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCa                                                                                  <  2:4678771/150‑1 (MQ=255)

                                                                    gCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCa                                                                                  <  1:3359563/150‑1 (MQ=255)

                                                                        gtcgACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCt                                                                              >  2:2915949/1‑150 (MQ=255)

                                                                                 gCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTa                                                                     <  2:619097/150‑1 (MQ=255)

                                                                                   gggTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGt                                                                   >  1:4462383/1‑150 (MQ=255)

                                                                                   gggTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGt                                                                   >  1:3721391/1‑150 (MQ=255)

                                                                                     gTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGc                                                                 >  1:1905393/1‑150 (MQ=255)

                                                                                       ttCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTg                                                               <  2:3007534/150‑1 (MQ=255)

                                                                                        tCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTgg                                                              >  1:121230/1‑150 (MQ=255)

                                                                                        tCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTgg                                                              >  1:1712995/1‑150 (MQ=255)

                                                                                          cAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCg                                                            >  2:4227187/1‑150 (MQ=255)

                                                                                              cGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCg                                                        <  2:1896239/150‑1 (MQ=255)

                                                                                                   gTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTc                                                   <  1:1093075/150‑1 (MQ=255)

                                                                                                    tCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcg                                                  >  1:3892917/1‑150 (MQ=255)

                                                                                                      gATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcgac                                                >  1:4703853/1‑150 (MQ=255)

                                                                                                       aTACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcgacg                                               >  2:4392013/1‑150 (MQ=255)

                                                                                                         aCGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGaa                                             <  1:3593322/150‑1 (MQ=255)

                                                                                                         aCGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGaa                                             >  2:2524943/1‑150 (MQ=255)

                                                                                                         aCGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGaa                                             <  1:581652/150‑1 (MQ=255)

                                                                                                          cGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAg                                            <  2:938329/150‑1 (MQ=255)

                                                                                                          cGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAg                                            <  1:2996710/150‑1 (MQ=255)

                                                                                                          cGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAg                                            >  1:3941424/1‑150 (MQ=255)

                                                                                                           ggCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAgg                                           >  1:203577/1‑150 (MQ=255)

                                                                                                           ggCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAgg                                           >  2:3458089/1‑150 (MQ=255)

                                                                                                            gCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGt                                          >  2:5034326/1‑150 (MQ=255)

                                                                                                                tgctgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCa                                      <  1:636986/150‑1 (MQ=255)

                                                                                                                  ctgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATa                                    <  1:4074827/150‑1 (MQ=255)

                                                                                                                   tgGCATTCGTTGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                   <  1:2279824/150‑1 (MQ=255)

                                                                                                                      cATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCTGCTTCGACGAAGGTTTCATAGGTg                                <  1:4588285/150‑1 (MQ=255)

                                                                                                                        ttCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGcc                              >  2:4521018/1‑150 (MQ=255)

                                                                                                                          cGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCt                            <  1:1870469/150‑1 (MQ=255)

                                                                                                                           gATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAg                                                                                                   >  1:4323960/1‑78 (MQ=255)

                                                                                                                           gATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAg                                                                                                   <  2:4323960/78‑1 (MQ=255)

                                                                                                                           gATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTg                           >  1:4091437/1‑150 (MQ=255)

                                                                                                                               ttCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACg                       >  2:1963768/1‑150 (MQ=255)

                                                                                                                                tCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACACACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGa                      <  2:1716899/150‑1 (MQ=255)

                                                                                                                                  aGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATg                    >  1:519685/1‑150 (MQ=255)

                                                                                                                                  aGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATg                    <  2:4856810/150‑1 (MQ=255)

                                                                                                                                    gcagcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGcc                  >  1:3777254/3‑150 (MQ=255)

                                                                                                                                    cGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGcc                  <  1:3111927/150‑1 (MQ=255)

                                                                                                                                     gAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGccc                 >  2:2061606/1‑150 (MQ=255)

                                                                                                                                     cagcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGccc                 <  1:2755689/149‑1 (MQ=255)

                                                                                                                                     cagcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGccc                 <  2:3023356/149‑1 (MQ=255)

                                                                                                                                     cagcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGccc                 >  2:579571/2‑150 (MQ=255)

                                                                                                                                      agcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCt                <  1:1297844/150‑1 (MQ=255)

                                                                                                                                      agcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCt                <  1:3260685/150‑1 (MQ=255)

                                                                                                                                      aGCAAGCCCCAGCCAGGGTCGCGGCCTGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAACGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCt                <  1:2685843/150‑1 (MQ=255)

                                                                                                                                      aGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGGAGGTTTCATAGGTGCCCTGGACGATGCCCt                >  1:76110/1‑150 (MQ=255)

                                                                                                                                      aGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGACGGTTTCATAGGTGCCCTGGACGATGCCCt                >  2:82573/1‑150 (MQ=255)

                                                                                                                                       gCAAGCACCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTg               <  1:4305922/150‑1 (MQ=255)

                                                                                                                                         aaGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGc             <  2:4462383/150‑1 (MQ=255)

                                                                                                                                         aaGCACCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTgtg             >  2:604542/1‑148 (MQ=255)

                                                                                                                                         aaGCACCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTgtg             >  2:604982/1‑148 (MQ=255)

                                                                                                                                            ccccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGc          <  1:3123604/150‑1 (MQ=255)

                                                                                                                                            ccccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGc          >  1:3312586/1‑150 (MQ=255)

                                                                                                                                            ccccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGc          >  2:933880/1‑150 (MQ=255)

                                                                                                                                             cccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGcc         <  1:93835/150‑1 (MQ=255)

                                                                                                                                                 gccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATg     >  1:3455139/1‑150 (MQ=255)

                                                                                                                                                    agGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAg  <  1:19551/150‑1 (MQ=255)

GGAGATGGCCTTGCCTTCGGCGGTGTACTTGGCGATGCGTGCGAGGGCGTCGTCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAG  >  NC_002947/5738713‑5739010

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