BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,107,744	C→G	5.3%	noncoding (35/76 nt)	PP_t21 →	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,107,744	0	C	G	5.3%	1112.6 / 38.3	339	noncoding (35/76 nt)	PP_t21	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base C (98/221); new base G (6/12); total (106/233)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 346 total aligned reads displayed.

GGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGT  >  NC_002947/2107595‑2107893

ggCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                                                                       <  1:7795734/150‑1 (MQ=255)

ggCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                                                                       >  1:8143490/1‑150 (MQ=255)

  cAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCa                                                                                                                                                     <  2:4705032/150‑1 (MQ=255)

    gggAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAg                                                                                                                                                   <  1:1648898/150‑1 (MQ=255)

      gAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGt                                                                                                                                                 <  2:1455665/150‑1 (MQ=255)

        gTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                                                                                                                               <  2:503888/150‑1 (MQ=255)

         tGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgt                                                                                                                                              <  2:704573/150‑2 (MQ=255)

         tGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                              <  1:3196920/150‑1 (MQ=255)

             aGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgctc                                                                                                                                          <  1:6941688/150‑6 (MQ=255)

             aGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTc                                                                                                                                          <  1:1689907/150‑1 (MQ=255)

             aGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTc                                                                                                                                          <  1:2847045/150‑1 (MQ=255)

             aGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTc                                                                                                                                          <  2:5801169/150‑1 (MQ=255)

               caggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                        >  2:4628802/1‑150 (MQ=255)

               caggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                        <  2:6360434/150‑1 (MQ=255)

               caggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGACCCCGACCTTGGCAAGGTCgtgctcta                                                                                                                                        <  1:3667050/150‑8 (MQ=17)

                 ggcatTAATGCCAGCAGCGTGTCCCCTTCGTCGAGTGGCCTAGGACACCGCCCTTTCTCGGCGGTAACAGGGGTTCGAGTCCCCTTGGGGACTCCAATGCGGGAATAGCTAAGTTGGTAGAGCACGACATTGCCAAGGTCGGGGTCGCGa                                                                                                                                      <  1:4759754/150‑1 (MQ=14)

                 ggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                                                      <  1:7538040/150‑1 (MQ=255)

                 ggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                                                      <  2:4124168/150‑1 (MQ=255)

                  gcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAg                                                                                                                                     >  1:3863796/1‑150 (MQ=255)

                     gTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                                                                                  <  1:4073258/150‑1 (MQ=255)

                       aaTGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                                <  1:5588902/150‑1 (MQ=255)

                       aaTGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                                <  1:1479563/150‑1 (MQ=255)

                         tGCCAGCAGCGTGTCCCCTTCGTCTCGGGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                                                              <  2:6472672/150‑1 (MQ=255)

                         tGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                                                              <  1:2945503/150‑1 (MQ=255)

                           ccAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtct                                                                                                                            >  1:4394098/1‑150 (MQ=40)

                           ccAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGTCCTTGCCAAGGTCGGGGTCGCGAGGTCGAGtct                                                                                                                            >  1:3953425/1‑150 (MQ=255)

                            cagcagCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                                                           <  2:3957105/150‑1 (MQ=39)

                                 gcGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                      <  1:5979454/150‑1 (MQ=35)

                                  cGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                     <  1:2000049/150‑1 (MQ=34)

                                        cccTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTcc                                                                                                               <  2:2105965/150‑1 (MQ=33)

                                           ttCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGt                                                                                                            <  1:1576924/150‑1 (MQ=33)

                                             cGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGttt                                                                                                          <  1:727970/150‑1 (MQ=33)

                                              gTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTtc                                                                                                         <  1:5185717/150‑1 (MQ=34)

                                                 tAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTcc                                                                                                      <  2:8533570/150‑1 (MQ=37)

                                                  aGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCg                                                                                                     <  1:4565117/150‑1 (MQ=38)

                                                  aGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCg                                                                                                     <  1:8007884/150‑1 (MQ=38)

                                                        cTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTg                                                                                               <  2:4329323/150‑1 (MQ=255)

                                                         aaGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGc                                                                                              <  2:6546919/149‑1 (MQ=255)

                                                           ggACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCtt                                                                                            >  2:746387/1‑150 (MQ=255)

                                                            gACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGTGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTa                                                                                           <  1:4157938/150‑1 (MQ=255)

                                                                  gCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGGATAGCGGAGGTGCGAGAGCACGAGCTTGCCGCGGTCGGGGTCGCGAGTTCGCGTCTCGTTGCCCGCTCCCGTGTCTCCGGCCCTGCGGGgcggcg                                                                                     >  2:7069836/1‑150 (MQ=14)

                                                                  gCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAgcggcg                                                                                     <  2:6250495/150‑1 (MQ=255)

                                                                   cccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                                                                    <  1:3111624/150‑1 (MQ=255)

                                                                    ccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCggg                                                                                   <  2:4217033/150‑1 (MQ=255)

                                                                      tttCACGGCGGTAACAGGGGTTCGAGTCCCCTGGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGc                                                                                 >  2:7627027/1‑150 (MQ=255)

                                                                      tttCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTGCGAGTCTCGATTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGc                                                                                 >  1:2354815/1‑150 (MQ=255)

                                                                       ttCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCa                                                                                <  1:6449543/150‑1 (MQ=255)

                                                                        tCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTAGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAg                                                                               <  2:6091854/150‑1 (MQ=255)

                                                                         cACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAgg                                                                              <  2:2976089/150‑1 (MQ=255)

                                                                          aCGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGga                                                                             >  2:2256855/1‑150 (MQ=255)

                                                                            ggcggTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaa                                                                           <  1:5279056/150‑1 (MQ=255)

                                                                            ggcggTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaa                                                                           >  2:3584470/1‑150 (MQ=255)

                                                                                 tAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaagaaa                                                                      <  2:7344601/150‑1 (MQ=255)

                                                                                 tAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaagaaa                                                                      <  2:5677858/150‑1 (MQ=255)

                                                                                   aCAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgc                                                                                                                                            <  2:5628732/78‑4 (MQ=255)

                                                                                   aCAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                                                                    <  1:3813246/134‑1 (MQ=255)

                                                                                   aCAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                                                                    >  2:3813246/1‑134 (MQ=255)

                                                                                     aGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCCCGACCTTGGCAAGGTCgtgc                                                                                                                                            >  2:527445/1‑73 (MQ=11)

                                                                                     aGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAAccc                                                                  <  1:3443540/150‑1 (MQ=255)

                                                                                      ggggTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCa                                                                 <  1:5823404/150‑1 (MQ=255)

                                                                                      ggggTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGGGCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTCGCGGCGGGGCAGGAAAAGAAAAccct                                                                 >  1:2625872/1‑149 (MQ=255)

                                                                                       gggTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                                                           >  1:1277015/1‑91 (MQ=35)

                                                                                       gggTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAg                                                                <  2:3261955/150‑1 (MQ=255)

                                                                                                tCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGGTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCGGGAAGAGAAAACCCAGGCTGACTCg                                                       >  1:1750658/1‑150 (MQ=255)

                                                                                                  cccTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTg                                                     <  2:5776453/150‑1 (MQ=255)

                                                                                                  cccTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTg                                                     <  1:6819848/150‑1 (MQ=255)

                                                                                                  cccTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCGCCAGTTACTCCGGCACTGCTGAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTg                                                     >  1:509423/1‑150 (MQ=255)

                                                                                                   ccgaGGGGACGCCTCTGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGa                                                    <  2:7695851/147‑1 (MQ=255)

                                                                                                   ccTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGa                                                    <  1:3501283/150‑1 (MQ=255)

                                                                                                    cTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAg                                                   <  1:7116767/150‑1 (MQ=255)

                                                                                                        gggACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGTGGtgc                                               >  2:5583497/1‑148 (MQ=255)

                                                                                                          gACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCAt                                             <  1:7441205/150‑1 (MQ=255)

                                                                                                             gCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                                >  1:8631781/1‑64 (MQ=35)

                                                                                                                aaTGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTg                                       <  1:973049/150‑1 (MQ=255)

                                                                                                                aaTGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTg                                       <  1:7618384/150‑1 (MQ=255)

                                                                                                                 aTGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTgg                                      <  2:2464722/150‑1 (MQ=255)

                                                                                                                  tGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGt                                     <  1:5005246/150‑1 (MQ=255)

                                                                                                                          tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGaa                             <  2:8143490/150‑1 (MQ=255)

                                                                                                                                aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAgtgt                       <  1:2109324/150‑1 (MQ=255)

                                                                                                                                aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAgtgt                       >  1:4649799/1‑150 (MQ=255)

                                                                                                                                 gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgtt                      <  1:1052982/150‑1 (MQ=255)

                                                                                                                                 gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgtt                      <  1:6264151/150‑1 (MQ=255)

                                                                                                                                  ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgttg                     <  2:6754312/150‑1 (MQ=255)

                                                                                                                                  ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgttg                     <  2:7315725/150‑1 (MQ=255)

                                                                                                                                   tGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgttgt                    <  2:8616221/150‑1 (MQ=255)

                                                                                                                                   tGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgttgt                    <  2:7828437/150‑1 (MQ=255)

                                                                                                                                   tGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgttgt                    <  1:2763691/150‑1 (MQ=255)

                                                                                                                                   tGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgttgt                    <  1:4440426/150‑1 (MQ=255)

                                                                                                                                   tGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgttgt                    <  1:8670976/150‑1 (MQ=255)

                                                                                                                                      tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTtgtgtg                 >  2:8341409/1‑150 (MQ=255)

                                                                                                                                          gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTccc             >  2:6811673/1‑150 (MQ=255)

                                                                                                                                          ccccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCa                                                                                <  2:6623538/80‑1 (MQ=255)

                                                                                                                                           cccGACCTTGGCAAGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTcccc            >  1:6524441/3‑150 (MQ=255)

                                                                                                                                           cccGACCTTGGCAAGGTCAGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTATCCGGCGCTGCTTCGCGGCGGTGCAGGAAACGAAAACCCCGGCAGACGGGTGAGGGTCATGTGCTGGTGCACTGGAAAGTGTGGTGTGTcacc            >  2:182203/3‑150 (MQ=21)

                                                                                                                                           cccGACCTTGGCAAGGTAGGGCGCGCGAGTTCGAGTCTCGTTACCCGCTCCAGTGTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGGGAGGTTCAGGTTCTGGGGCACTGAGGGGTGGTGTGTGTcccc            >  2:2170728/3‑150 (MQ=255)

                                                                                                                                              gACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTc         <  2:2337659/150‑1 (MQ=255)

                                                                                                                                              gACCGTGCCACGGTCGGGGTCGCGAGTTCGAGTCTGGTGTCCCGCTCCAGTTGCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTc         <  2:122660/150‑1 (MQ=39)

                                                                                                                                                  ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGTGAAGTTAACCAAGGGTGAATCGTGTGGTTCTGGTTGTGGGGGACGGAGAAGTGTTGTGTGTCCCCGTGGTCt     >  2:2883559/1‑150 (MQ=21)

                                                                                                                                                   tGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTa    <  2:7209950/150‑1 (MQ=255)

                                                                                                                                                    gCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGG‑AAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGt  <  1:8525202/150‑1 (MQ=255)

GGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGT  >  NC_002947/2107595‑2107893

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