BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,345,576	Δ1 bp	5.2%	coding (244/1587 nt)	ybiT ←	ABC‑F family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,345,576	0	T	.	5.2%	1001.3 / 43.8	287	coding (244/1587 nt)	ybiT	ABC‑F family ATPase
Reads supporting (aligned to +/- strand): ref base T (137/135); new base . (3/12); total (140/147)
Fisher's exact test for biased strand distribution p-value = 3.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 291 total aligned reads displayed.

GGCGGTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGATCA  >  NC_002947/3345427‑3345725

   ggTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGGCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCAGGATCACCGTG‑TCGATCAc                                                                                                                                                    >  2:5219180/1‑150 (MQ=255)

    gtgtAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACg                                                                                                                                                   <  2:908896/150‑1 (MQ=255)

    gtgtAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACg                                                                                                                                                   <  2:3154182/150‑1 (MQ=255)

       tAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTg                                                                                                                                                <  2:1446850/150‑1 (MQ=255)

       gaGCCGTCCATGTCGGCGAATTCGGGTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTg                                                                                                                                                <  2:579152/149‑1 (MQ=255)

        aGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGA‑CACGGTGa                                                                                                                                               <  1:4906103/150‑1 (MQ=255)

        aGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGA‑CACGGTGa                                                                                                                                               <  1:2075475/150‑1 (MQ=255)

          ccGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGA‑CACGGTGatc                                                                                                                                             <  1:5420720/150‑3 (MQ=255)

           cGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGA‑CACGGTGatca                                                                                                                                            <  2:1342675/150‑4 (MQ=255)

            gTCCATTTCGGTGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTc                                                                                                                                           >  2:312810/1‑150 (MQ=255)

            gTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTc                                                                                                                                           >  1:5161933/1‑150 (MQ=255)

            gTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTc                                                                                                                                           >  2:3929655/1‑150 (MQ=255)

             tCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTcc                                                                                                                                          <  2:427554/150‑1 (MQ=255)

               cATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTc                                                                                                                                        <  1:1747925/150‑1 (MQ=255)

               cATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTc                                                                                                                                        <  1:1316719/150‑1 (MQ=255)

                 tttCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGt                                                                                                                                      >  2:3568433/1‑150 (MQ=255)

                  ttCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTa                                                                                                                                     <  1:2352114/150‑1 (MQ=255)

                  ttCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTa                                                                                                                                     <  1:2352150/150‑1 (MQ=255)

                  ttCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCCTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACCGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTa                                                                                                                                     <  2:3267591/150‑1 (MQ=255)

                   tCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAg                                                                                                                                    <  1:5070385/150‑1 (MQ=255)

                      gCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGACTTCCTCGTAGGCg                                                                                                                                 >  2:2582983/1‑150 (MQ=255)

                       cGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGa                                                                                                                                <  1:3878053/150‑1 (MQ=255)

                       cGAATTCGGTGTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGa                                                                                                                                <  1:3377233/150‑1 (MQ=255)

                          aTTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACt                                                                                                                             >  1:570918/1‑150 (MQ=255)

                           ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                            >  2:2925990/1‑150 (MQ=255)

                           ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                            <  2:5060330/150‑1 (MQ=255)

                           ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                            >  2:4559747/1‑150 (MQ=255)

                             cGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTa                                                                                                                                     >  1:632987/1‑139 (MQ=255)

                             cGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGt                                                                                                                          >  2:1240851/1‑150 (MQ=255)

                               gTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTCGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTcc                                                                                                                        >  2:1994807/1‑150 (MQ=255)

                               gTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTcc                                                                                                                        <  2:4940693/150‑1 (MQ=255)

                                   ccAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGgc                                                                                                                    <  1:2607735/150‑1 (MQ=255)

                                        tCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCAGGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGc                                                                                                               >  2:3712886/1‑150 (MQ=255)

                                            cGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGc                                                                                                           <  2:554599/150‑1 (MQ=255)

                                                  ccATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGc                                                                                                     <  2:5165516/150‑1 (MQ=255)

                                                  ccATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGAAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGc                                                                                                     >  2:4789075/1‑150 (MQ=255)

                                                   cATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCg                                                                                                    >  1:136635/1‑150 (MQ=255)

                                                    aTGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCgg                                                                                                   <  1:4943259/150‑1 (MQ=255)

                                                       ccGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTa                                                                                                >  2:1698533/1‑150 (MQ=255)

                                                           cttcttCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCg                                                                                            >  2:2305503/1‑150 (MQ=255)

                                                              cttcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCt                                                                                         >  2:3213476/1‑150 (MQ=255)

                                                               ttcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTc                                                                                        <  1:1002588/150‑1 (MQ=255)

                                                               ttcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTc                                                                                        >  2:4071288/1‑150 (MQ=255)

                                                               ttcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTc                                                                                        >  1:2744765/1‑150 (MQ=255)

                                                               ttcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTc                                                                                        >  1:3575102/1‑150 (MQ=255)

                                                               ctcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTc                                                                                        <  1:36295/149‑1 (MQ=255)

                                                                 cGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGa                                                                                      <  2:2241466/150‑1 (MQ=255)

                                                                 cGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGa                                                                                      <  1:2532689/150‑1 (MQ=255)

                                                                 cGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGa                                                                                      <  1:3071788/150‑1 (MQ=255)

                                                                 cGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGa                                                                                      >  1:4673573/1‑150 (MQ=255)

                                                                  ggTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAg                                                                                     >  2:93196/1‑150 (MQ=255)

                                                                  ggTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAg                                                                                     >  2:610338/1‑150 (MQ=255)

                                                                    tCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCa                                                                                   >  2:5271882/1‑150 (MQ=255)

                                                                      ctttCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATg                                                                                 <  1:380763/149‑1 (MQ=255)

                                                                      aTTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCAGCTCGAGCATg                                                                                 <  1:4353748/150‑1 (MQ=255)

                                                                        ttCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGAc                                                                               >  1:4264015/1‑150 (MQ=255)

                                                                             gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                          <  2:3901507/150‑1 (MQ=255)

                                                                             gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                          <  2:3901995/150‑1 (MQ=255)

                                                                             gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                          <  1:4090500/150‑1 (MQ=255)

                                                                                  gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                     >  2:266/1‑150 (MQ=255)

                                                                                  gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                     >  2:3035231/1‑150 (MQ=255)

                                                                                  gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                     <  1:4479291/150‑1 (MQ=255)

                                                                                  gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                     >  1:924794/1‑150 (MQ=255)

                                                                                    gTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGa                                                                   <  1:84222/150‑1 (MQ=255)

                                                                                     tAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGAc                                                                  <  2:3381630/150‑1 (MQ=255)

                                                                                       gATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACgg                                                                >  2:656181/1‑150 (MQ=255)

                                                                                            ggTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCGGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCa                                                           >  2:4850451/1‑150 (MQ=255)

                                                                                                aCGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTc                                                       >  1:3250464/1‑150 (MQ=255)

                                                                                                 cGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCa                                                      >  1:1121368/1‑150 (MQ=255)

                                                                                                  gCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCat                                                     <  1:4198731/150‑2 (MQ=255)

                                                                                                  gCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCAc                                                     >  1:1793947/1‑150 (MQ=255)

                                                                                                     cGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCAccgc                                                  <  2:2438316/150‑1 (MQ=255)

                                                                                                      ggCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCAccgcc                                                 >  2:1871328/1‑150 (MQ=255)

                                                                                                        cTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGa                                               <  2:2797311/150‑1 (MQ=255)

                                                                                                         tttGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAg                                              >  2:4634465/1‑150 (MQ=255)

                                                                                                                ttCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTc                                       <  2:4413682/150‑1 (MQ=255)

                                                                                                                 tCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCa                                      <  1:2166395/150‑1 (MQ=255)

                                                                                                                 tCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCa                                      <  2:5108384/150‑1 (MQ=255)

                                                                                                                   cacaGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATg                                    >  2:3376806/1‑150 (MQ=255)

                                                                                                                      aGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAg                                 <  1:1058550/150‑1 (MQ=255)

                                                                                                                       gCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAgg                                >  2:757193/1‑150 (MQ=255)

                                                                                                                          ggCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCg                             >  2:1118611/1‑150 (MQ=255)

                                                                                                                          ggCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCg                             >  1:3783670/1‑150 (MQ=255)

                                                                                                                                ggCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGc                       <  2:4141094/150‑1 (MQ=255)

                                                                                                                                ggCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACAGCCGAGGATTTTCATGAAGGTCGACTTGc                       <  2:825654/150‑1 (MQ=255)

                                                                                                                                ggCCCATGATCACCGTG‑TCGATCACGGTGAATTCCCCGTAGCCGAACGGGTCCGGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCAGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGc                       <  2:3655210/150‑1 (MQ=255)

                                                                                                                                   ccATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTc                                       <  1:3036360/131‑1 (MQ=255)

                                                                                                                                   ccATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGc                    >  2:3501614/1‑150 (MQ=255)

                                                                                                                                   ccATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGc                    >  1:1571749/1‑150 (MQ=255)

                                                                                                                                      tGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAAc                 <  2:241774/150‑1 (MQ=255)

                                                                                                                                        ttCACCGTGATCGA‑CACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCg               >  1:2392411/2‑150 (MQ=255)

                                                                                                                                        aTCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCg               <  1:312974/150‑1 (MQ=255)

                                                                                                                                         tCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGt              >  1:1096476/1‑150 (MQ=255)

                                                                                                                                          cACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGtt             <  2:4484/150‑1 (MQ=255)

                                                                                                                                           aCCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTg            <  2:1147369/150‑1 (MQ=255)

                                                                                                                                                g‑tCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGcc       <  1:3929655/150‑1 (MQ=255)

                                                                                                                                                   cGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGa     <  2:3829593/150‑1 (MQ=255)

                                                                                                                                                   cGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGa     <  2:2143069/150‑1 (MQ=255)

                                                                                                                                                    gATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGAt    >  2:4184590/1‑150 (MQ=255)

                                                                                                                                                    gATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGAt    <  1:474646/150‑1 (MQ=255)

GGCGGTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTG‑TCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGATCA  >  NC_002947/3345427‑3345725

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