BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,104,717	C→G	5.1%	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,104,717	0	C	G	5.1%	501.6 / 13.4	157	noncoding (35/76 nt)	PP_t21	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base C (62/86); new base G (2/6); total (65/92)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 160 total aligned reads displayed.

CAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGT  >  NC_002947/2104570‑2104866

cAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCa                                                                                                                                                     <  2:1396993/150‑1 (MQ=255)

 aGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCaa                                                                                                                                                    >  1:578484/1‑150 (MQ=255)

 aGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCaa                                                                                                                                                    <  2:2731206/150‑1 (MQ=255)

       tGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgt                                                                                                                                              <  1:608313/150‑2 (MQ=255)

       tGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgt                                                                                                                                              <  1:1287473/150‑2 (MQ=255)

       tGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                              <  1:5188062/150‑1 (MQ=255)

       tGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                              <  2:2778147/150‑1 (MQ=255)

        gAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtg                                                                                                                                             <  1:3670806/150‑3 (MQ=255)

        gAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                             <  2:4132830/150‑1 (MQ=255)

        gAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                             <  1:968585/150‑1 (MQ=255)

         aaaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                                                                            >  1:5051997/1‑150 (MQ=255)

          aaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGt                                                                                                                                           >  2:4993649/1‑150 (MQ=255)

          aaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGt                                                                                                                                           >  1:2901661/1‑150 (MQ=255)

          aaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTAGAGTCCCCTAGTGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGt                                                                                                                                           >  1:2902130/1‑150 (MQ=255)

             caggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                        <  1:5138237/150‑1 (MQ=255)

             caggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                        <  1:3117769/150‑1 (MQ=255)

             caggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                        <  1:4545355/150‑1 (MQ=255)

             caggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                        <  2:4854684/150‑1 (MQ=255)

               ggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                                                      <  1:1038532/150‑1 (MQ=255)

               ggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                                                      <  1:2590737/150‑1 (MQ=255)

                gcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAg                                                                                                                                     >  1:1644494/1‑150 (MQ=255)

                gcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAg                                                                                                                                     >  1:3860680/1‑150 (MQ=255)

                   gTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                                                                                  <  2:1832322/150‑1 (MQ=255)

                   gTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                                                                                  <  2:4616098/150‑1 (MQ=255)

                    tAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTGGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                                                                                 >  1:3865104/1‑150 (MQ=255)

                     aaTGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                                <  1:4174101/150‑1 (MQ=255)

                      aTGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                                                                               <  2:1608846/150‑1 (MQ=255)

                       tGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                                                              >  2:1324814/1‑150 (MQ=255)

                        gCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTGCGAGTCCCCTCGGGGCCGCACAAGGC‑CGAATAGCTCAGTTGGTATAGCACGAACTTGCCAAGGTCGGGGTCGCGGGGTCGAGtg                                                                                                                             >  1:1999506/1‑149 (MQ=11)

                          cagcagCGGGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTGCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                                                           <  2:3780492/150‑1 (MQ=255)

                               gcGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                      <  2:813240/150‑1 (MQ=35)

                                 gtgtCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCCCGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTccc                                                                                                                    <  2:2027489/150‑1 (MQ=18)

                                          tCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCAAGTCCCCTAGGGGACGC‑CACTGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGtt                                                                                                           <  2:1722991/150‑1 (MQ=18)

                                                     ccTAGGACACCGCCCTTTCACGGCGGTAAAAGGTGTTCGAGTGCCCTAGGGGACTC‑CAATGCGGGAATAGCTGAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACt                                                                                                >  2:2170062/1‑150 (MQ=14)

                                                      cTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTg                                                                                               <  1:2272208/150‑1 (MQ=255)

                                                      cTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTg                                                                                               >  1:4873239/1‑150 (MQ=255)

                                                          gACACAGCCCTTTCAAGGCGGTAAAAGGTGTTCGAGTCACCTAGGGGACGC‑CAATGAGGGAAAAGCTCAGGTGGTAGAGCACGACCTTGCCAAGGTCGGTGTCGCGAGTTCGAGTCTCGTTTCCCGCTACAGTTTCTCCTGCACAGCTTa                                                                                           >  2:102007/1‑150 (MQ=14)

                                                           acacCGCCCGTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CACTGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAg                                                                                          <  2:2825336/150‑1 (MQ=255)

                                                                 cccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                                                                    <  2:3445020/150‑1 (MQ=255)

                                                                 cccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                                                                    <  2:2490282/150‑1 (MQ=255)

                                                                 cccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTGGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                                                                    <  2:1672114/150‑1 (MQ=255)

                                                                  ccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                                <  2:4038279/105‑1 (MQ=35)

                                                                  ccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                                >  1:4038279/1‑105 (MQ=35)

                                                                   cTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                     <  2:4374584/115‑1 (MQ=34)

                                                                   cTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                     >  1:4374584/1‑115 (MQ=34)

                                                                   cTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCgggg                                                                                  <  1:1324814/150‑1 (MQ=255)

                                                                     ttCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCa                                                                                <  1:1607582/150‑1 (MQ=255)

                                                                       cACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAgg                                                                              <  1:3059884/150‑1 (MQ=255)

                                                                          ggcggTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaa                                                                           <  2:3647627/150‑1 (MQ=255)

                                                                          ggcggTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTCCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaa                                                                           <  2:5256273/150‑1 (MQ=255)

                                                                           gcggTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaa                                                                          <  2:3860680/150‑1 (MQ=255)

                                                                                 aCAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAAc                                                                    >  2:3050187/1‑150 (MQ=255)

                                                                                   aGGGGTTTGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAAccc                                                                  >  2:3619874/1‑150 (MQ=255)

                                                                                   aGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAAccc                                                                  >  2:4184712/1‑150 (MQ=255)

                                                                                   aGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGATCGAGTCTCGTTTCAAGCTCCAGAATCTCCCGCACTGCGTAGCGGCGGGGCAGGAAAAGAAAAccc                                                                  >  1:1783546/1‑150 (MQ=18)

                                                                                     gggTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAg                                                                <  1:1038931/150‑1 (MQ=255)

                                                                                        ttCGAGTCCCCTAGGGGACGC‑CAATTCGGGAATAGCTCAGTTGGTAGAGCAAGACCTTCACAAGCTCGCCGTCGAGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGGGGGAAAAGACAGCCCggggg                                                             >  1:3244830/1‑148 (MQ=17)

                                                                                           gAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGaa                                                          <  1:5175816/150‑1 (MQ=255)

                                                                                               ccccTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGt                                                      <  1:4848773/150‑1 (MQ=255)

                                                                                                cccTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTg                                                     >  2:3817190/1‑150 (MQ=255)

                                                                                                 ccTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTTGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGa                                                    <  2:2790031/150‑1 (MQ=255)

                                                                                                 ccTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCTGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGa                                                    <  2:1643151/150‑1 (MQ=255)

                                                                                                 ccTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCt                                                                                             <  1:3244225/109‑1 (MQ=255)

                                                                                                 ccTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGa                                                    <  2:2701161/150‑1 (MQ=255)

                                                                                                      gggACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTc                                               <  1:1048808/150‑1 (MQ=255)

                                                                                                      gggACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTc                                               <  2:4389675/150‑1 (MQ=255)

                                                                                                           gC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                                                      >  1:871731/1‑58 (MQ=35)

                                                                                                           gC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                                                      <  2:871731/58‑1 (MQ=35)

                                                                                                           gC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGtt                                          >  1:158359/1‑150 (MQ=255)

                                                                                                              cAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                                                                                                        >  2:4585098/1‑70 (MQ=31)

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

                                                                                                                        aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGa                              >  1:2212206/1‑150 (MQ=255)

                                                                                                                          aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGaaa                            <  1:2759531/150‑1 (MQ=255)

                                                                                                                            cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaa                                                                          <  2:3498668/102‑1 (MQ=255)

                                                                                                                            cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaa                                                                          >  1:3498668/1‑102 (MQ=255)

                                                                                                                             tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAgt                         >  1:1146353/1‑150 (MQ=255)

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

                                                                                                                                gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGTGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgtt                      <  1:102007/150‑1 (MQ=255)

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

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

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

                                                                                                                                   ggTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTtcaa                                                                                                       >  2:3877978/1‑64 (MQ=14)

                                                                                                                                    gTAGAGCACGACCTTGCCATGGTCGGGGTCGCGAGTTGGAGTCTCGTTTCCCGGTCCTGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGATAACCCAGGCTTATTCGTGAGGTTGATGTTCTTGTGCTCTGTAAAGTGTtgtgt                  <  2:745767/150‑1 (MQ=255)

                                                                                                                                    gTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTGGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGgtgtgt                  >  1:4247282/1‑150 (MQ=255)

                                                                                                                                      agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTgtgtgt                >  2:3097648/1‑150 (MQ=255)

                                                                                                                                       gaccccGACCTTGGCAAGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTgg                                      >  2:5154032/6‑127 (MQ=255)

                                                                                                                                          cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGaa                             >  2:2925757/3‑133 (MQ=255)

                                                                                                                                          cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGaa                             <  1:2925757/131‑1 (MQ=255)

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

                                                                                                                                             gACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGATTCCCGATCCAGTATCTCCTCCACTGCTAAGCGGCGGGGCAGGAAAATAAAACCAAGGCAGAATAGTGGGGGTCATGTTAGGGTGCACTGACAAGTGATGTGTGTCCCCTTc         >  2:495898/1‑150 (MQ=25)

                                                                                                                                               ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGt       <  2:578484/150‑1 (MQ=255)

                                                                                                                                               ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTACAGTTTCTCAGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTGTGGTGCACTGAATTGTGTTGTGTGTCCCCTTCGt       >  2:5428219/1‑150 (MQ=38)

                                                                                                                                               ccTTGCCAAGGTCGGGGTCGCGAGTTAGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTTTGTGTCCCCTTCGt       <  2:2902130/150‑1 (MQ=255)

                                                                                                                                                cTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTc      >  1:4935940/1‑150 (MQ=255)

                                                                                                                                                 ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCt     >  1:5201432/1‑150 (MQ=255)

                                                                                                                                                 ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCt     >  1:5151677/1‑150 (MQ=255)

                                                                                                                                                   gCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTc                                               <  2:3938041/106‑1 (MQ=255)

                                                                                                                                                    ccAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGt  <  2:3043429/150‑1 (MQ=255)

                                                                                                                                                    ccAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGt  <  2:1644494/150‑1 (MQ=255)

                                                                                                                                                    ccAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGt  <  1:3817190/150‑1 (MQ=255)

CAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGC‑CAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGT  >  NC_002947/2104570‑2104866

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