BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,763,637	A→C	5.5%	K338T (AAG→ACG) ‡	PP_5057 →	M23/M37 family peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,763,637	0	A	C	5.5%	345.8 / 9.3	109	K338T (AAG→ACG) ‡	PP_5057	M23/M37 family peptidase
Reads supporting (aligned to +/- strand): ref base A (50/52); new base C (2/4); total (52/57)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

Only 100 of 109 total aligned reads displayed.

CCGTACCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGC  >  NC_002947/5763491‑5763782

ccGTACCACCCTCGGGCCTATGGTTTCCAGCGATGGCGCAAACTACGGCGGCGAATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGTTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGt                                                                                                                                                 <  3:224310/149‑1 (MQ=255)

ccGTACCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGt                                                                                                                                                 <  5:59239/149‑1 (MQ=255)

ccGTACCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGt                                                                                                                                                 <  5:88449/149‑1 (MQ=255)

     ccaccCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGa                                                                                                                                             <  2:252516/148‑1 (MQ=255)

      caccCTCGGGCCGATGGTTTCCAGCGCTGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGGCAGTCTATGGTCGATTGCTGGCACGTGTCGGTGATGCCCGCGGTGGCGACGCGCGTGCCAAGGGGGACg                                                                                                                                           >  8:330250/1‑149 (MQ=255)

      caccCTCGGGCCGATGGTTTCCAGCCATGGCGCCAAATACGGCGGCGCATTTTCTGATGCGCGTGGAAAACTTCCATGGCAAGTCAATGGTCGATTGCTGCCACGTTTCGGTGATGCCCGCTGTAGCGACGCCCGTGCCAAGTGGCACg                                                                                                                                           >  2:374021/1‑149 (MQ=255)

        cccTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACggg                                                                                                                                         >  6:314876/1‑149 (MQ=255)

        cccTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACggg                                                                                                                                         <  5:430324/149‑1 (MQ=255)

        cccTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACggg                                                                                                                                         <  3:98222/149‑1 (MQ=255)

         ccTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACgggg                                                                                                                                        <  4:252197/149‑1 (MQ=255)

         ccTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACgggc                                                                                                                                        <  6:78241/149‑2 (MQ=255)

           tCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGt                                                                                                                                       >  6:252248/1‑148 (MQ=255)

           tCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGt                                                                                                                                       >  5:353980/1‑148 (MQ=255)

              ggCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATg                                                                                                                                   <  2:380998/149‑1 (MQ=255)

              ggCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTCCCACTTGGCACgggcgtcgc                                                                                                                                   <  3:139357/149‑7 (MQ=255)

               gCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGa                                                                                                                                  >  2:266412/1‑149 (MQ=255)

                 cGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGAGGCCCGTGCCAAGTGGGACGGGGTTATGATc                                                                                                                                <  1:38934/149‑1 (MQ=255)

                 cGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGAt                                                                                                                                 >  8:225442/1‑148 (MQ=255)

                 cGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATc                                                                                                                                <  7:88884/149‑1 (MQ=255)

                      gTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAgcgc                                                                                                                           >  1:281422/1‑149 (MQ=255)

                       tttCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGTCGGGGTTATGATCAgcgc                                                                                                                           >  3:97771/1‑148 (MQ=255)

                         tCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCaa                                                                                                                        >  8:373218/1‑149 (MQ=255)

                          ccAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAc                                                                                                                       >  3:27714/1‑149 (MQ=255)

                            aGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAccc                                                                                                                     >  7:218412/1‑149 (MQ=255)

                            aGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAccc                                                                                                                     >  7:212307/1‑149 (MQ=255)

                            aGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAcc                                                                                                                      >  5:169945/1‑148 (MQ=255)

                             gCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAAcc                                                                                                                      >  8:378734/1‑147 (MQ=255)

                             gCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCa                                                                                                                    <  5:271487/149‑1 (MQ=255)

                                aTGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCAAgg                                                                                                                  >  3:47922/1‑148 (MQ=255)

                                         aCTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTa                                                                                                        <  8:212307/149‑1 (MQ=255)

                                         aCTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTa                                                                                                        <  8:218412/149‑1 (MQ=255)

                                          cTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCCGGCACCCAGGTAc                                                                                                       <  4:349565/149‑1 (MQ=255)

                                           tACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcg                                                                                                      >  5:57509/1‑149 (MQ=255)

                                             cggcggCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCAGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgt                                                                                                     <  4:47922/148‑1 (MQ=255)

                                              ggcggcGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgc                                                                                                   >  8:119588/1‑149 (MQ=255)

                                                cggcgCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgccg                                                                                                 >  8:242165/1‑149 (MQ=255)

                                                cggcgCATTTTCTGCTGCGCGGGGCAAACTTCAATGTCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgccg                                                                                                 <  7:302020/149‑1 (MQ=255)

                                                 ggcgCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgccgt                                                                                                <  5:252248/149‑1 (MQ=255)

                                                 ggcgCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTAcgtgccgt                                                                                                >  5:124970/1‑149 (MQ=255)

                                                    gcATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCa                                                                                             >  1:57822/1‑149 (MQ=255)

                                                    gcATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGCTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCa                                                                                             <  7:395239/149‑1 (MQ=255)

                                                       ttttCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCAcgg                                                                                          <  5:314876/149‑1 (MQ=255)

                                                                 cgcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCgtgt                                                                                 <  1:266412/148‑1 (MQ=255)

                                                                 cgcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGtgtg                                                                                <  5:364189/149‑1 (MQ=255)

                                                                  gcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTgg                                                                               <  8:64008/149‑1 (MQ=255)

                                                                  gcgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTgg                                                                               <  4:247862/149‑1 (MQ=255)

                                                                   cgGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTgg                                                                               >  6:216378/1‑148 (MQ=255)

                                                                    ggggCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTa                                                                             >  7:350841/1‑149 (MQ=255)

                                                                     gggCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTAt                                                                            >  1:413051/1‑149 (MQ=255)

                                                                           aCTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCgg                                                                      <  5:216378/149‑1 (MQ=255)

                                                                             ttCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGAc                                                                    <  2:281422/149‑1 (MQ=255)

                                                                               ccATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGAc                                                                    <  7:242165/147‑1 (MQ=255)

                                                                                cATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTgg                                                                 <  6:124970/149‑1 (MQ=255)

                                                                                cATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTg                                                                  <  6:432657/148‑1 (MQ=255)

                                                                                 aTGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCTGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGt                                                                >  3:360967/1‑149 (MQ=255)

                                                                                    gCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCTGTAGCGACGCCCGTGCCATGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGATTGGTTgc                                                             <  1:374021/149‑1 (MQ=255)

                                                                                    gCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGtgtg                                                                                >  7:316462/1‑130 (MQ=255)

                                                                                    gCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGtgtg                                                                                <  8:316462/130‑1 (MQ=255)

                                                                                       aGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGcgcg                                                          <  2:64118/149‑1 (MQ=255)

                                                                                       aGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGcgcg                                                          <  2:57822/149‑1 (MQ=255)

                                                                                       aGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGcgcg                                                          <  7:373218/149‑1 (MQ=255)

                                                                                          cAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTg                                                       >  6:439980/1‑149 (MQ=255)

                                                                                           aaTGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTg                                                       <  2:413051/148‑1 (MQ=255)

                                                                                           aaTGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGc                                                      >  8:404786/1‑149 (MQ=255)

                                                                                              ggTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCgg                                                   >  5:243729/1‑149 (MQ=255)

                                                                                                tCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGc                                                 <  7:119588/149‑1 (MQ=255)

                                                                                                  gATTGCTGGCACGTTTCGGTGATGCCCGCGGGAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCtt                                               >  1:267650/1‑149 (MQ=255)

                                                                                                   aTTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                              >  7:144915/1‑149 (MQ=255)

                                                                                                      gCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACATGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTgg                                           >  4:88741/1‑149 (MQ=255)

                                                                                                        tGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                              <  2:267178/144‑1 (MQ=255)

                                                                                                        tGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                              >  1:267178/1‑144 (MQ=255)

                                                                                                        tGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGt                                          >  5:387549/1‑148 (MQ=255)

                                                                                                          gCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCa                                        <  1:57757/148‑1 (MQ=255)

                                                                                                           cACGTTTCGGTGATGCCCGCGGTAGCGACGGCCGTGCCCAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCAtt                                      <  3:314808/149‑1 (MQ=255)

                                                                                                            aCGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGATATCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATtc                                     >  1:363399/1‑149 (MQ=255)

                                                                                                              gTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                              <  2:105647/138‑1 (MQ=255)

                                                                                                              gTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTc                                              >  1:105647/1‑138 (MQ=255)

                                                                                                               tttCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATtctc                                   >  4:410995/1‑148 (MQ=255)

                                                                                                                 tCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGAc                                >  3:169257/1‑149 (MQ=255)

                                                                                                                     tGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATg                            <  6:169945/149‑1 (MQ=255)

                                                                                                                      gATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                           <  6:243729/149‑1 (MQ=255)

                                                                                                                      gATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAAGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                           >  2:341199/1‑149 (MQ=255)

                                                                                                                       aTGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCATGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                           <  1:341199/148‑1 (MQ=255)

                                                                                                                         gCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTc                                         >  7:80000/1‑132 (MQ=255)

                                                                                                                         gCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTc                                         <  8:80000/132‑1 (MQ=255)

                                                                                                                         gCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCaa                        >  4:275182/1‑149 (MQ=255)

                                                                                                                          cccGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                           <  1:70136/145‑1 (MQ=255)

                                                                                                                          cccGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATgg                           >  2:70136/1‑145 (MQ=255)

                                                                                                                            cgcgGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACgg                     >  7:87568/1‑149 (MQ=255)

                                                                                                                             gcgGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGt                    <  7:378734/149‑1 (MQ=255)

                                                                                                                              cgGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTTTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGtt                   <  4:27714/149‑1 (MQ=255)

                                                                                                                                gTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGTGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTAc                 <  7:143321/149‑1 (MQ=255)

                                                                                                                                    cGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCt               >  3:101348/1‑147 (MQ=255)

                                                                                                                                    ataaccccGTCCCACTTGGCACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCt               >  6:393121/6‑147 (MQ=255)

                                                                                                                                     taaccccGTCCCACTTGGCACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAg            >  6:222303/5‑149 (MQ=255)

                                                                                                                                     gACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAg            >  2:318966/1‑149 (MQ=255)

                                                                                                                                      aCGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGc           <  6:57509/149‑1 (MQ=255)

                                                                                                                                      aCGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGc           >  6:36973/1‑149 (MQ=255)

                                                                                                                                        gCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCt         <  6:124792/149‑1 (MQ=255)

                                                                                                                                               ccAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGc  <  4:169257/149‑1 (MQ=255)

CCGTACCACCCTCGGGCCGATGGTTTCCAGCGATGGCGCCAACTACGGCGGCGCATTTTCTGCTGCGCGGGGCAAACTTCCATGGCCAGTCAATGGTCGATTGCTGGCACGTTTCGGTGATGCCCGCGGTAGCGACGCCCGTGCCAAGTGGGACGGGGTTATGATCAGCGCCAACCCAGGCACCCAGGTACGTGCCGTGCACGGCGGGCGTGTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGC  >  NC_002947/5763491‑5763782

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