BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,738,914	T→A	9.9%	V170V (GTA→GTT) ‡	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,738,914	0	T	A	9.9%	896.8 / 71.4	304	V170V (GTA→GTT) ‡	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base T (112/162); new base A (7/23); total (119/185)
Fisher's exact test for biased strand distribution p-value = 7.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01

Only 100 of 307 total aligned reads displayed.

TCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAG  >  NC_002947/5738765‑5739063

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

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

  gAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTAc                                                                                                                                                     <  2:8272379/150‑1 (MQ=255)

        aGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACtt                                                                                                                                               <  2:4622937/150‑1 (MQ=255)

        aGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGAGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACtt                                                                                                                                               >  2:277660/1‑150 (MQ=255)

            gCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCt                                                                                                                                           >  2:4810287/1‑150 (MQ=255)

            gCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCt                                                                                                                                           <  2:5774919/150‑1 (MQ=255)

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

               tGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGATGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTc                                                                                                                                        >  1:6988153/1‑150 (MQ=255)

                gCTCGTCGACGTAGCGGGTTTCCAGGCGGATGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCa                                                                                                                                       <  2:7774421/150‑1 (MQ=255)

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

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

                    gtcgACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCt                                                                                                                                   <  2:5424792/150‑1 (MQ=255)

                        aCGTAGCGGGTTTCCAGGCGGAAGTCGATATGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTgccg                                                                                                                               >  2:679895/1‑150 (MQ=255)

                             gCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTa                                                                                                                          >  1:4241286/1‑150 (MQ=255)

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

                             gCGGGGGTCCAGGCGGCCGTCGATCCGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTa                                                                                                                          <  1:4581667/150‑1 (MQ=255)

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

                                ggTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTg                                                                                                                       >  1:4853409/1‑150 (MQ=255)

                                 gTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGc                                                                                                                      <  2:6622102/150‑1 (MQ=255)

                                  tttCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCt                                                                                                                     <  2:8637410/150‑1 (MQ=255)

                                  tttCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCt                                                                                                                     >  1:5479964/1‑150 (MQ=255)

                                   ttCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGCTGGTCCGCGAGCAAGCCCCAGCCCGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTg                                                                                                                    <  1:142863/150‑1 (MQ=255)

                                   ttCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGGTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCCGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTg                                                                                                                    <  2:6207272/150‑1 (MQ=255)

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

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

                                       aGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGa                                                                                                                >  2:9086176/1‑150 (MQ=255)

                                         gCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGAcc                                                                                                              <  1:7213880/150‑1 (MQ=255)

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

                                                tCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcg                                                                                                       <  2:6541413/150‑1 (MQ=255)

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

                                                    tACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcgacga                                                                                                   >  2:1312383/1‑150 (MQ=255)

                                                    tACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcgacga                                                                                                   <  2:7619967/150‑1 (MQ=255)

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

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

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

                                                         ctctGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGtt                                                                                              <  1:1774641/150‑1 (MQ=255)

                                                            tgctgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTc                                                                                                                                        <  1:6621659/105‑1 (MQ=255)

                                                              ctgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATa                                                                                         >  2:874644/1‑150 (MQ=255)

                                                               tgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                                                                        >  1:9147825/1‑150 (MQ=255)

                                                               tgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                                                                        <  1:9034990/150‑1 (MQ=255)

                                                               tgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                                                                        >  2:766251/1‑150 (MQ=255)

                                                               tgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                                                                        <  1:414682/150‑1 (MQ=255)

                                                                  cATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTg                                                                                     <  1:868725/150‑1 (MQ=255)

                                                                    ttCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGcc                                                                                   <  1:51821/150‑1 (MQ=255)

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

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

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

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

                                                                         tGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGa                                                                              <  1:1607652/150‑1 (MQ=255)

                                                                          gTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGAc                                                                             <  1:3156766/150‑1 (MQ=255)

                                                                            tCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGa                                                                           <  1:8393492/150‑1 (MQ=255)

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

                                                                               gCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGc                                                                        >  1:8466007/1‑150 (MQ=255)

                                                                                    cAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTgg                                                                   <  1:7988581/150‑1 (MQ=255)

                                                                                    cAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTgg                                                                   >  1:7078174/1‑150 (MQ=255)

                                                                                    cAAGCACCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTgt                                                                   <  2:4755473/150‑2 (MQ=255)

                                                                                     aaGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGc                                                                  >  1:5041108/1‑150 (MQ=255)

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

                                                                                         cccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGcc                                                              >  1:4521663/1‑150 (MQ=255)

                                                                                           cagccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGa                                                            <  1:6521379/150‑1 (MQ=255)

                                                                                           cagccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGa                                                            <  2:1910915/150‑1 (MQ=255)

                                                                                              ccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGt                                                         <  1:5892602/150‑1 (MQ=255)

                                                                                               cagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTa                                                        >  2:380243/1‑150 (MQ=255)

                                                                                                 gggTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGa                                                      <  1:3673061/150‑1 (MQ=255)

                                                                                                 gggTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGa                                                      >  2:7847212/1‑150 (MQ=255)

                                                                                                  ggTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGAt                                                     >  1:8717010/1‑150 (MQ=255)

                                                                                                     cgcgGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCa                                                  <  1:2783469/150‑1 (MQ=255)

                                                                                                      gcgGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAg                                                 <  2:806790/150‑1 (MQ=255)

                                                                                                      gcgGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAg                                                 >  2:6425485/1‑150 (MQ=255)

                                                                                                           cAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGcc                                            >  2:1029035/1‑150 (MQ=255)

                                                                                                                gCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCcggcgg                                       <  1:5626829/150‑1 (MQ=255)

                                                                                                                 cTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGt                                      <  1:6163966/150‑1 (MQ=255)

                                                                                                                  tgtgctccgccCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTc                                     <  2:8543352/144‑1 (MQ=255)

                                                                                                                  tGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGt                                                         >  2:5904598/1‑130 (MQ=255)

                                                                                                                   gggCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCa                                    <  1:1910029/150‑1 (MQ=255)

                                                                                                                   gggCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCa                                    <  2:2956359/150‑1 (MQ=255)

                                                                                                                    ggCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCAt                                   <  2:4490113/150‑1 (MQ=255)

                                                                                                                    ggCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCAt                                   <  2:409584/150‑1 (MQ=255)

                                                                                                                    ggCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCAt                                   >  2:7389546/1‑150 (MQ=255)

                                                                                                                       tccgccCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTg                                >  2:6058904/2‑150 (MQ=255)

                                                                                                                          gccCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGcc                             >  2:2458512/1‑150 (MQ=255)

                                                                                                                            ccATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCAt                           <  1:3851752/150‑1 (MQ=255)

                                                                                                                             cATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATa                          <  1:8083458/150‑1 (MQ=255)

                                                                                                                              aTGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATAc                         <  1:8965906/150‑1 (MQ=255)

                                                                                                                                gccgccAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACAt                       <  2:4977158/150‑1 (MQ=255)

                                                                                                                                gccgccAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACAt                       <  2:2701140/150‑1 (MQ=255)

                                                                                                                                   gccCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTc                                     >  2:1468995/1‑133 (MQ=255)

                                                                                                                                     ccAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCa                  <  1:873914/150‑1 (MQ=255)

                                                                                                                                     caaGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGcga                  <  1:2384275/148‑3 (MQ=255)

                                                                                                                                     caaGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGcga                  <  1:339731/148‑3 (MQ=255)

                                                                                                                                     caaGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGcga                  <  1:6415481/148‑3 (MQ=255)

                                                                                                                                      cAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAg                 <  1:4023757/150‑1 (MQ=255)

                                                                                                                                       aGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAga                >  1:6812368/1‑149 (MQ=255)

                                                                                                                                       aGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAga                <  1:1113172/150‑2 (MQ=255)

                                                                                                                                           cAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCt            <  1:1621628/150‑1 (MQ=255)

                                                                                                                                            aGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCtt           <  1:766251/150‑1 (MQ=255)

                                                                                                                                                gTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCg       <  1:7802555/150‑1 (MQ=255)

                                                                                                                                                    gTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCa   <  2:2999831/150‑1 (MQ=255)

                                                                                                                                                    gTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCa   <  1:6342318/150‑1 (MQ=255)

TCGAGGTCAGTGGCCTGCTCGTCGACGTAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAG  >  NC_002947/5738765‑5739063

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