BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,530,410	C→T	5.8%	R67R (CGC→CGT)	PP_4858 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,530,410	0	C	T	5.8%	745.9 / 24.2	268	R67R (CGC→CGT)	PP_4858	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (81/160); new base T (11/4); total (102/166)
Fisher's exact test for biased strand distribution p-value = 3.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 277 total aligned reads displayed.

TCGCCCTGCCGGCGTCCGAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTC  >  NC_002947/5530262‑5530559

tCGCCCTGCCGGCGTCCGAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGcc                                                                                                                                                      <  1:3197319/150‑1 (MQ=255)

  gCCCTGCCGGCGTCCGAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTg                                                                                                                                                    >  1:7140770/1‑150 (MQ=255)

        ccGGCGTCCGAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGaa                                                                                                                                              <  2:7386142/150‑1 (MQ=255)

              tCCGAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCt                                                                                                                                        <  2:4896521/150‑1 (MQ=255)

              tCCGAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCt                                                                                                                                        >  1:5457758/1‑150 (MQ=255)

               ccGAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCGTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTg                                                                                                                                       <  2:4443358/150‑1 (MQ=255)

               ccGAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTg                                                                                                                                       <  1:1934455/150‑1 (MQ=255)

                 gAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGc                                                                                                                                     >  2:8010905/1‑150 (MQ=255)

                     aaCACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCcgcg                                                                                                                                 >  1:6514179/1‑150 (MQ=255)

                     aaCACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCcgcg                                                                                                                                 >  1:4385983/1‑150 (MQ=255)

                     aaCACCTGACGTCAGAAGCGCGCTCCAGAGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCcgcg                                                                                                                                 >  1:6632457/1‑150 (MQ=255)

                         ccTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGc                                                                                                                             <  1:1674452/150‑1 (MQ=255)

                             ccGTCAGAAGCGCTCTGCAGCTCTGACGTAATCGCCCAGCCGCTTACCACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGaa                                                                                                                         <  2:3216023/149‑1 (MQ=17)

                               gTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAgg                                                                                                                       >  1:2926571/1‑150 (MQ=255)

                                tCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGaa                                                                                                                         <  1:7309368/147‑1 (MQ=255)

                                tCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGt                                                                                                                      <  2:2366101/150‑1 (MQ=255)

                                  aGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTgg                                                                                                                    <  1:3031084/150‑1 (MQ=255)

                                  aGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTgg                                                                                                                    <  2:5530641/150‑1 (MQ=255)

                                  aGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTgg                                                                                                                    <  1:4771033/150‑1 (MQ=255)

                                  aGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGAGTTGCTGAAGGTgg                                                                                                                    >  1:4148362/1‑150 (MQ=255)

                                    aaGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCt                                                                                                                  >  2:5458264/1‑150 (MQ=255)

                                        gcgcTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGAGTTGCAGAAGGTGGCTCCGa                                                                                                              >  1:6835737/1‑150 (MQ=255)

                                         cgcTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGAt                                                                                                             <  1:4742393/150‑1 (MQ=255)

                                          gcTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                                                                                                            <  1:4596474/150‑1 (MQ=255)

                                          gcTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                                                                                                            <  2:8449000/150‑1 (MQ=255)

                                            tCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGc                                                                                                          <  2:7087397/150‑1 (MQ=255)

                                             ccAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGcc                                                                                                         <  1:250078/150‑1 (MQ=255)

                                              cAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCt                                                                                                        <  1:6116580/150‑1 (MQ=255)

                                                gcgcTGCCGAAATAGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGa                                                                                                      <  2:6302485/150‑1 (MQ=255)

                                                 cgcTGCCGAAATCGACCAGCCGCTGAACACAGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGAt                                                                                                     <  1:8372353/150‑1 (MQ=255)

                                                       cGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAGGCGGCTCCGATCGTCTGa                                                                                                      >  2:1647529/1‑143 (MQ=11)

                                                       cGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAGGCGGCTCCGATCGTCTGa                                                                                                      <  1:1647529/143‑1 (MQ=11)

                                                        gAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGt                                                                                              >  1:7037085/1‑150 (MQ=255)

                                                         aaaTCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGtt                                                                                             >  1:1414767/1‑150 (MQ=255)

                                                         aaaTCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGtt                                                                                             <  1:2403881/150‑1 (MQ=255)

                                                              gCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGc                                                                                        <  1:3575574/150‑1 (MQ=255)

                                                              gCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGCGTTGCAGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGc                                                                                        >  1:2894187/1‑150 (MQ=255)

                                                                ccAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGAt                                                                                                     >  2:4988178/1‑135 (MQ=255)

                                                                 cAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCgt                                                                                     <  1:2328760/150‑1 (MQ=255)

                                                                 cAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCgt                                                                                     <  1:3941375/150‑1 (MQ=255)

                                                                 cAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGGCTGATCGAACGTTCTGGCCgg                                                                                     >  1:857071/1‑149 (MQ=255)

                                                                      gCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTAGCAGAAGGTGGCTCCGATCGCCTGATCGAACGTTATGGCCGAGTTGCGGAAGGAGGCTCGGATCGGCTGATCGAACTATCTGGCCGGGTTGc                                                                                >  2:5871967/1‑150 (MQ=14)

                                                                        tGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTg                                                                              <  1:5693591/150‑1 (MQ=255)

                                                                          aaCACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGaa                                                                            <  2:4151702/150‑1 (MQ=255)

                                                                             acCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGt                                                                         <  1:799325/150‑1 (MQ=255)

                                                                               cGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTATGGCCGCGTTGCTGAAGGAGGCTCCGATCGACTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGGGTTGCTGTAGgggg                                                                       >  2:7708856/1‑150 (MQ=14)

                                                                                  aGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTc                                                                    <  1:1331274/150‑1 (MQ=40)

                                                                                       caGGCGGCTCCGAGCGCCTGATCGATCGTTCGGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                               <  2:1831447/149‑1 (MQ=255)

                                                                                       caGGCGGCTCCGAGCGCCAGATCGAACGTTCGGGCCGCGTTGCTGAAGGGGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                               <  2:5562636/149‑1 (MQ=255)

                                                                                       aaGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                               <  1:5246828/150‑1 (MQ=39)

                                                                                       aaGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTAGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                               >  1:5834451/1‑150 (MQ=255)

                                                                                         ggcggcTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGt                                                             <  1:4956669/150‑1 (MQ=39)

                                                                                          gcggcTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTc                                                            <  2:2692327/150‑1 (MQ=38)

                                                                                          gcggcTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTc                                                            <  1:5093797/150‑1 (MQ=18)

                                                                                              cTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCGGAAGGTGGCTCCGATCGCCTGCTCGAACGCGCTGGCCGGCTTGCCGAACGTGGCTCCGAGCGTcggcg                                                        >  1:7429947/1‑146 (MQ=11)

                                                                                               tCCGATCGCCTGTTCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGTCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATc                                                       <  2:8080416/150‑1 (MQ=255)

                                                                                               tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATc                                                       >  1:7533824/1‑150 (MQ=37)

                                                                                                  gatgcccTGCTCGATCGTGCTTTCAGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGaa                                                    <  2:7014326/145‑1 (MQ=11)

                                                                                                  gATCGCCTGTCCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGaa                                                    <  2:2108055/150‑1 (MQ=255)

                                                                                                      gCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGATCTGGCCGCGTTGCTGAAGGTGCCTCCGATCGACTGATCGAACgat                                                >  2:6517904/1‑148 (MQ=255)

                                                                                                        cTGATCGAACGTTCTCGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCt                                              <  2:403922/150‑1 (MQ=37)

                                                                                                        cTGATCGAACGTTCGGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCt                                              <  2:5191894/150‑1 (MQ=37)

                                                                                                         tGATCGACCGGTCTGGCCGCGTTGCTGACGGGGGCTCCGCTCGGCTGATCGAACGTTCTGGCCGCGTTGCTGCAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTg                                             <  2:7980891/150‑1 (MQ=11)

                                                                                                           aTCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGc                                           <  2:680861/150‑1 (MQ=37)

                                                                                                              gAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCcgc                                        <  2:674626/150‑1 (MQ=38)

                                                                                                              gAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGAGTTGCTGAAGGTGGCTCCGATCGGCTGATCGAACTTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCcgg                                        >  2:3841020/1‑149 (MQ=12)

                                                                                                               aaCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCcgcg                                       <  2:8562102/150‑1 (MQ=38)

                                                                                                                   ttCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCGGATCGCCTGGTTGAActca                                                                                             <  2:1353889/92‑4 (MQ=11)

                                                                                                                        gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCGGATCGCCt                                                                                                        <  2:8703826/76‑1 (MQ=255)

                                                                                                                        gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCt                                                                                                        >  1:8703826/1‑76 (MQ=39)

                                                                                                                        gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGCAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCGGacg                              >  1:1057311/1‑148 (MQ=255)

                                                                                                                        gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAg                              <  2:5329136/150‑1 (MQ=255)

                                                                                                                          cgcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGt                            <  1:7796927/150‑1 (MQ=255)

                                                                                                                          cgcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGATCTGGCCGCGTTGCTGAAGGt                            >  2:7600181/1‑150 (MQ=255)

                                                                                                                            cgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTgg                          <  1:5540513/150‑1 (MQ=255)

                                                                                                                               tGCTGAATGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAgggacagg                       <  2:4529925/150‑7 (MQ=255)

                                                                                                                               tGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGCTCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTc                       <  2:2926825/150‑1 (MQ=255)

                                                                                                                                 cTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCCGAAGGTGGCTCCg                     >  1:950347/1‑150 (MQ=38)

                                                                                                                                  tGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGGGTTGCTGAAGGTGGCTCCGa                    >  2:3353350/1‑150 (MQ=255)

                                                                                                                                  tGAAGGCGGCTCCGATCGCCTGATCGAACGGTCTGGCCGCGTTGCTGAAGGTGGCTCCGCTCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGa                    <  2:6313225/150‑1 (MQ=255)

                                                                                                                                   gAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACTTTCTGGCCGTGTTGCGGAAGGTGGCTCCGATCGTCTGATCGAACTATCTGGCCGAGTAGCTGAAGGTGGCTCCGAt                   >  2:5966639/1‑150 (MQ=11)

                                                                                                                                    aaGGTGGCTCCGATCGCCTGATCGACCGGTCTCGCCGCGTTGCTGAAGGGGGCTCCGATCGCCTGATCGACCGGTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                  <  2:7431899/150‑1 (MQ=11)

                                                                                                                                    aaGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                  <  2:5697589/150‑1 (MQ=39)

                                                                                                                                    aaGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                  <  1:7909809/150‑1 (MQ=39)

                                                                                                                                      ggTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGt                <  1:2217329/150‑1 (MQ=39)

                                                                                                                                       gTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTc               >  1:19824/1‑150 (MQ=38)

                                                                                                                                       gTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTc               <  1:8252727/150‑1 (MQ=38)

                                                                                                                                       gTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTc               >  2:3985557/1‑150 (MQ=38)

                                                                                                                                       gTGGCGCGGATCGGCTGTTCGAACGTTCTGGCCGCGTTGCTGTAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTc               <  2:3576660/150‑1 (MQ=11)

                                                                                                                                        tGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGGGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCGGAAGGTGGCTCCGATCGTCt              >  2:1730137/1‑150 (MQ=255)

                                                                                                                                        cggCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCt              <  2:1153418/149‑1 (MQ=25)

                                                                                                                                          gCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGa            <  2:4257265/150‑1 (MQ=37)

                                                                                                                                          gCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGa            <  1:5183597/150‑1 (MQ=37)

                                                                                                                                           cTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGAt           >  1:5392481/1‑150 (MQ=37)

                                                                                                                                            tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGGGTTGCTGAAGGTGGCTCCGATCGTCTGATc          >  2:4594714/1‑150 (MQ=17)

                                                                                                                                            tCCGATCGCCTGAGCGAACGTTCTGGCCGCGTTGCTGAAGGGGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATc          <  2:2273430/150‑1 (MQ=255)

                                                                                                                                               gATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCCGAAg                                                                           >  2:7243177/1‑82 (MQ=18)

                                                                                                                                                aTCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGGGTGGCTGAAGGTGGCTCCGATCGTCTGATCGAAc      >  2:2613188/1‑150 (MQ=21)

                                                                                                                                                   gCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGtt   <  1:1567591/150‑1 (MQ=37)

                                                                                                                                                    ccTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTc  >  1:7351489/1‑150 (MQ=37)

TCGCCCTGCCGGCGTCCGAGCAACACCTGACGTCAGAAGCGCGCTCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTC  >  NC_002947/5530262‑5530559

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`