BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,738,944	G→A	8.1%	H160H (CAC→CAT)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,738,944	0	G	A	8.1%	541.5 / 27.0	184	H160H (CAC→CAT)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base G (73/96); new base A (7/8); total (80/104)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 185 total aligned reads displayed.

CGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCA  >  NC_002947/5738795‑5739090

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

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

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

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

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

       ccAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGc                                                                                                                                             <  1:2330632/150‑1 (MQ=255)

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

          ggCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCAACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGAc                                                                                                                                          <  1:3907582/150‑1 (MQ=255)

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

             ggAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCgg                                                                                                                                       >  1:2502212/1‑150 (MQ=255)

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

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

                  tCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCggacatg                                                                                                                                  <  2:2234658/150‑6 (MQ=255)

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

                   cGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcga                                                                                                                                 <  1:2214603/150‑1 (MQ=255)

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

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

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

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

                         ggCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGTCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAgg                                                                                                                           >  2:1466514/1‑150 (MQ=255)

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

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

                            tctGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGttt                                                                                                                        <  1:4603089/150‑1 (MQ=255)

                               gctgGCATTCGATGTTCAGCGAGCAAGCCCCGGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCAt                                                                                                                     >  1:2090542/1‑150 (MQ=255)

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

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

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

                                        cgttGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCt                                                                                                            <  1:129942/147‑1 (MQ=255)

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

                                               cAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGAt                                                                                                     >  2:1296976/1‑150 (MQ=255)

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

                                                 ggcagcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCATGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGc                                                                                                   >  2:292782/4‑150 (MQ=255)

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

                                                   gAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGccc                                                                                                 <  1:3705089/150‑1 (MQ=255)

                                                   gAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGccc                                                                                                 <  1:962995/150‑1 (MQ=255)

                                                   gAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGccc                                                                                                 <  1:1162070/150‑1 (MQ=255)

                                                    cgCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCt                                                                                                <  1:1057631/149‑1 (MQ=255)

                                                    aGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCt                                                                                                <  2:2570120/150‑1 (MQ=255)

                                                    aGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCt                                                                                                <  1:1057684/150‑1 (MQ=255)

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

                                                      cAAGCACCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTGGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTgt                                                                                              >  2:3408698/1‑149 (MQ=255)

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

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

                                                           cccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCCCTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTGCGACGACGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGcc                                                                                         >  2:3865902/1‑150 (MQ=255)

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

                                                            ccagccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTg                                                                                                                                                  >  1:5260669/1‑92 (MQ=255)

                                                             cagccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGTGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGa                                                                                       <  1:2106354/150‑1 (MQ=255)

                                                             cagccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGa                                                                                       >  2:4078441/1‑150 (MQ=255)

                                                             cagccagGGTCGCGGCCAGGGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGGACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGa                                                                                       <  2:5275251/150‑1 (MQ=255)

                                                              agccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGAt                                                                                      <  1:4373523/150‑1 (MQ=255)

                                                                ccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGt                                                                                    >  2:1718948/1‑150 (MQ=255)

                                                                 cagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTa                                                                                   <  2:967373/150‑1 (MQ=255)

                                                                  agGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAg                                                                                  >  1:2421720/1‑150 (MQ=255)

                                                                      tCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATcc                                                                              <  2:1455628/150‑1 (MQ=255)

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

                                                                                 ggCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCcggcg                                                                   <  1:83682/150‑1 (MQ=255)

                                                                                  gCTGGGCTCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCcggcgg                                                                  <  2:1788904/150‑1 (MQ=255)

                                                                                     gggCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGtt                                                                                                                         <  1:2603214/92‑1 (MQ=255)

                                                                                           cgccCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGc                                                         >  2:3016421/1‑150 (MQ=255)

                                                                                            gccCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGcc                                                        <  1:3703838/150‑1 (MQ=255)

                                                                                              ccATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCAt                                                      <  2:1162885/150‑1 (MQ=255)

                                                                                              ccATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCAt                                                      >  1:4772162/1‑150 (MQ=255)

                                                                                                   gcgccAAGGCCACCGGTCAGAATCCATTTGCCCTTGTGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACAtg                                                 >  2:4664198/2‑149 (MQ=255)

                                                                                                   ccgccCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATc                                                 <  2:3309833/150‑1 (MQ=255)

                                                                                                   ccgccCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATc                                                 <  2:65678/150‑1 (MQ=255)

                                                                                                    cgccCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCg                                                >  1:2476077/1‑150 (MQ=255)

                                                                                                       ccAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCcgg                                                                     >  1:4102421/1‑126 (MQ=255)

                                                                                                         aGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAgg                                           <  2:2476077/150‑1 (MQ=255)

                                                                                                           gCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCGATGTATGGCCAgatg                                         <  2:295785/150‑4 (MQ=255)

                                                                                                               ccGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCttc                                     <  2:4180041/149‑2 (MQ=255)

                                                                                                                 ggTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGc                                   <  2:4015103/150‑1 (MQ=255)

                                                                                                                  gTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGTCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCg                                  <  2:3596721/150‑1 (MQ=255)

                                                                                                                    gAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTc                                >  2:4552244/1‑150 (MQ=255)

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

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

                                                                                                                        aCCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGt                            >  1:4855892/1‑150 (MQ=255)

                                                                                                                         cccACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGtt                           >  2:3932345/1‑150 (MQ=255)

                                                                                                                          ccACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTc                          >  1:309770/1‑150 (MQ=255)

                                                                                                                           cACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCTAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCg                         >  2:1942991/1‑150 (MQ=255)

                                                                                                                           cACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCg                         >  2:1010360/1‑150 (MQ=255)

                                                                                                                            aCTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGt                        <  2:5383023/150‑1 (MQ=255)

                                                                                                                              ttGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTc                          >  2:2708243/1‑146 (MQ=255)

                                                                                                                              ttGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTc                          <  1:2708243/146‑1 (MQ=255)

                                                                                                                               tGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGa                     <  2:2964232/150‑1 (MQ=255)

                                                                                                                                gCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGaa                    <  1:253467/150‑1 (MQ=255)

                                                                                                                                 cccttgaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGaaa                   <  2:5430788/144‑2 (MQ=255)

                                                                                                                                 cccttgaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGaaa                   <  1:3992904/144‑2 (MQ=255)

                                                                                                                                 ccTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAg                   <  2:80388/150‑1 (MQ=255)

                                                                                                                                   tttCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTTGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCTGAGGTCATCTGGCCATACATCGCAAGGCCCATGGCGACCAGTGCGTTGAAttg                 >  1:5394533/1‑150 (MQ=255)

                                                                                                                                   tttCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAgtg                 >  2:140532/1‑150 (MQ=255)

                                                                                                                                   tttCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCAAACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAgtg                 >  1:283650/1‑150 (MQ=255)

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

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

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

                                                                                                                                      cAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTc              <  2:5071676/150‑1 (MQ=255)

                                                                                                                                        ggCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTccc            <  1:4995542/150‑1 (MQ=255)

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

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

                                                                                                                                                  gTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCtt                                      <  1:1293611/114‑1 (MQ=255)

                                                                                                                                                  gTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCa  <  1:205447/150‑1 (MQ=255)

CGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCA  >  NC_002947/5738795‑5739090

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