BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_008463	2,787,877	C→G	pseudogene (955/3090 nt)	PA14_RS13085 →	CusA/CzcA family heavy metal efflux RND transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_008463	2,787,877	0	C	G	99.7%	1303.4 / NA	330	pseudogene (955/3090 nt)	PA14_RS13085	CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand): ref base C (0/1); major base G (145/183); minor base T (1/0); total (146/184)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

CATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCC‑A‑GGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑G‑TGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAAC  >  NC_008463/2787741‑2788023

cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑Tgaagaa                                                                                                                                        >  2:7085524/1‑150 (MQ=255)

cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑Tgaagaa                                                                                                                                        >  1:115673/1‑150 (MQ=255)

 aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑Tgaagaag                                                                                                                                       <  2:353052/150‑1 (MQ=255)

 aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑Tgaagaag                                                                                                                                       >  2:5091084/1‑150 (MQ=255)

    cTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCTTGGCCACGG‑TGAAGAAGAAc                                                                                                                                    >  2:6875958/1‑150 (MQ=255)

     tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAAcc                                                                                                                                   >  2:5924295/1‑150 (MQ=255)

     tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAAcc                                                                                                                                   <  1:7527251/150‑1 (MQ=255)

     tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAAcc                                                                                                                                   <  2:6074452/150‑1 (MQ=255)

      gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAAcct                                                                                                                                  <  2:10253475/150‑2 (MQ=255)

       ctCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAAcctg                                                                                                                                 <  1:7332361/149‑3 (MQ=255)

        tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAAcctgg                                                                                                                                <  1:9680769/150‑4 (MQ=255)

        tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAAcctgg                                                                                                                                <  2:158879/150‑4 (MQ=255)

        tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAAcctgg                                                                                                                                <  1:5103316/150‑4 (MQ=255)

         cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAAcctggt                                                                                                                               >  1:10409157/1‑146 (MQ=255)

             gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGaa                                                                                                                          >  1:8701677/1‑150 (MQ=255)

              agaACAGCCGGGACGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAg                                                                                                                         <  2:4377660/150‑1 (MQ=255)

               gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAgg                                                                                                                        <  2:2881623/150‑1 (MQ=255)

               gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAgg                                                                                                                        <  1:2584817/150‑1 (MQ=255)

               gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAgg                                                                                                                        >  1:2936645/1‑150 (MQ=255)

                aaCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggc                                                                                                                       <  1:2724145/150‑1 (MQ=255)

                aaCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggc                                                                                                                       >  2:7563017/1‑150 (MQ=255)

                   aGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACCAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggcgcg                                                                                                                    >  2:7231009/1‑150 (MQ=255)

                   aGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggcgcg                                                                                                                    >  2:4757088/1‑150 (MQ=255)

                   aGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggcgcg                                                                                                                    >  2:2159371/1‑150 (MQ=255)

                    gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggcgcgg                                                                                                                   >  1:9216421/1‑150 (MQ=255)

                    gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggcgcgg                                                                                                                   >  2:8139063/1‑150 (MQ=255)

                    gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggcgcgg                                                                                                                   <  1:6359744/150‑1 (MQ=255)

                      cGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggcgcggcg                                                                                                                 <  1:8627936/150‑1 (MQ=255)

                       gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggcgcggcgc                                                                                                                <  1:1371503/150‑1 (MQ=255)

                       gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAggcgcggcgc                                                                                                                <  2:6744925/150‑1 (MQ=255)

                         gAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTg                                                                                                              <  1:2827080/150‑1 (MQ=255)

                          aaGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCTCCCTGGTGGACAAGGGC‑GTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTgg                                                                                                             <  2:2631520/150‑1 (MQ=255)

                           aGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGt                                                                                                            >  2:4157741/1‑150 (MQ=255)

                            gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTg                                                                                                           >  2:8390878/1‑150 (MQ=255)

                              ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTAGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGAt                                                                                                         >  2:6513290/1‑150 (MQ=255)

                                  caggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcgcg                                                                                                     >  2:3298359/4‑150 (MQ=255)

                                  caggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGAGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcgcg                                                                                                     >  1:5621181/4‑150 (MQ=255)

                                   aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCgg                                                                                                    >  2:2755897/3‑150 (MQ=255)

                                   aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCgg                                                                                                    <  2:4255833/148‑1 (MQ=255)

                                   aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCgg                                                                                                    >  1:9461001/3‑150 (MQ=255)

                                    ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGt                                                                                                   >  1:6164472/2‑150 (MQ=255)

                                    ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGt                                                                                                   <  1:2347966/149‑1 (MQ=255)

                                        gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctc                                                                                               >  2:3108235/1‑150 (MQ=255)

                                         tcggCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctct                                                                                              <  2:1047435/147‑1 (MQ=255)

                                          cggCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCtt                                                                                             >  2:887433/3‑150 (MQ=255)

                                            gCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                                                                           <  1:10155696/150‑1 (MQ=255)

                                             ccAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                                                                          <  2:1163803/150‑1 (MQ=255)

                                              cAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTg                                                                                         >  2:5253300/1‑150 (MQ=255)

                                              cAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTg                                                                                         <  2:4194234/150‑1 (MQ=255)

                                              cAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTg                                                                                         >  2:3670681/1‑150 (MQ=255)

                                               aGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                                                                        >  1:7459318/1‑150 (MQ=255)

                                                 cgcCTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTc                                                                                      <  1:839161/150‑1 (MQ=255)

                                                   ccTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCt                                                                                    >  2:6223411/1‑150 (MQ=255)

                                                   ccTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGAC‑ATGGCCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCt                                                                                    >  1:4779166/1‑150 (MQ=255)

                                                    cTGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTc                                                                                   <  1:8103450/150‑1 (MQ=255)

                                                     tGGAGGAGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                                                                  <  2:1835008/150‑1 (MQ=255)

                                                      ggaggaGATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                                                                 <  2:7771620/150‑1 (MQ=255)

                                                       gaggagATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGc                                                                                >  1:9855859/1‑150 (MQ=255)

                                                       gaggagATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGc                                                                                <  1:7622995/150‑1 (MQ=255)

                                                       gaggagATCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGc                                                                                <  1:5477092/150‑1 (MQ=255)

                                                            gaTCAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACAt                                                                           <  2:768909/150‑1 (MQ=255)

                                                              tCACCCGGACCCTG‑CCC‑ACGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCAC‑GCTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGAGCTCTTCCTGTTCCTCGGCAACATcc                                                                         <  2:1431904/150‑1 (MQ=255)

                                                               cAACCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                                                                        <  2:5522030/150‑1 (MQ=255)

                                                                aaCCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgc                                                                       >  2:9409267/1‑150 (MQ=255)

                                                                aaCCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCGGTTCATCGGCAACATCcgc                                                                       >  2:2474490/1‑150 (MQ=255)

                                                                 aCCGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg                                                                      >  1:3752950/1‑150 (MQ=255)

                                                                  ccGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc                                                                     <  1:4463935/150‑1 (MQ=255)

                                                                   cGGACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg                                                                    <  2:9182735/150‑1 (MQ=255)

                                                                    ggACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCACCATCCGCgcgg                                                                   >  2:4630603/1‑150 (MQ=255)

                                                                     gACCCTGCCCC‑A‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                                                                  <  2:3876354/150‑1 (MQ=255)

                                                                     gACCCTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                                                                  <  1:3165117/150‑1 (MQ=255)

                                                                        ccTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCt                                                               <  1:1696574/150‑1 (MQ=255)

                                                                         cTG‑CCCAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCGACATCCGCGCGGCGCTg                                                              <  1:9553264/150‑1 (MQ=255)

                                                                            gcccaa‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGAt                                                            <  1:6594604/145‑1 (MQ=255)

                                                                            gcccaa‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGAt                                                            <  1:4230130/145‑1 (MQ=255)

                                                                                 a‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg                                                       <  1:893615/150‑1 (MQ=255)

                                                                                 a‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg                                                       <  1:7127666/150‑1 (MQ=255)

                                                                                       gtcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGTCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacc                                                  >  2:10340869/7‑150 (MQ=255)

                                                                                        tcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacca                                                 >  1:10188869/6‑150 (MQ=255)

                                                                                        tcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacca                                                 <  1:9806935/145‑1 (MQ=255)

                                                                                         caaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcat                                                <  1:6781261/146‑1 (MQ=255)

                                                                                         caaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcat                                                >  1:2109000/5‑150 (MQ=255)

                                                                                           aggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatca                                              >  1:7768887/3‑150 (MQ=255)

                                                                                           aggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatca                                              <  1:6414742/148‑1 (MQ=255)

                                                                                                aTCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCCTCATCCCg                                         <  1:5996934/150‑1 (MQ=255)

                                                                                                    ccGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGt                                     <  2:6443035/150‑1 (MQ=255)

                                                                                                    ccGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGt                                     <  1:10449640/150‑1 (MQ=255)

                                                                                                     cGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                                    <  2:8225448/150‑1 (MQ=255)

                                                                                                      gTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCg                                   <  1:329079/150‑1 (MQ=255)

                                                                                                        cTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGAt                                 <  1:303220/150‑1 (MQ=255)

                                                                                                         tACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATg                                >  2:3660436/1‑150 (MQ=255)

                                                                                                         tACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATg                                >  2:7496537/1‑150 (MQ=255)

                                                                                                          aCGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGc                               >  1:5210580/1‑150 (MQ=255)

                                                                                                          aCGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGc                               >  1:356515/1‑150 (MQ=255)

                                                                                                             aCCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt                            <  2:86730/150‑1 (MQ=255)

                                                                                                                 caccacCCTGGTGGACCAGG‑CCGTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttcac                        >  2:8014657/1‑150 (MQ=255)

                                                                                                                                  aaggccgTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTc        <  2:5606678/144‑1 (MQ=255)

                                                                                                                                  aaggccgTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTc        >  2:9509558/7‑150 (MQ=255)

                                                                                                                                      ccgTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCa    >  1:1271575/3‑150 (MQ=255)

                                                                                                                                        gTGGCCACGG‑TGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAAc  <  1:7939766/150‑1 (MQ=255)

CATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCC‑A‑GGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑G‑TGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAAC  >  NC_008463/2787741‑2788023

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