BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,163,459	T→G	26.5%	L159R (CTG→CGG)	PP_3661 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,163,459	0	T	G	26.5%	300.5 / 35.9	200	L159R (CTG→CGG)	PP_3661	membrane protein
Reads supporting (aligned to +/- strand): ref base T (68/73); new base G (41/16); total (111/89)
Fisher's exact test for biased strand distribution p-value = 2.69e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTC  >  NC_002947/4163315‑4163607

cGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCggcgcggcgcc                                                                                                                                                 >  1:2635484/1‑149 (MQ=255)

cGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCa                                                                                                                                                 <  1:3891152/150‑1 (MQ=255)

cGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCa                                                                                                                                                 <  1:42878/150‑1 (MQ=255)

  gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                               >  2:2769853/1‑150 (MQ=255)

     cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcc                                                                                                                                            <  2:395774/150‑1 (MQ=255)

     cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcc                                                                                                                                            >  2:2665025/1‑150 (MQ=255)

       ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccga                                                                                                                                          <  2:4677818/150‑3 (MQ=255)

       ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccc                                                                                                                                          <  1:886665/150‑1 (MQ=255)

       ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccc                                                                                                                                          >  2:2989966/1‑150 (MQ=255)

       ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGcccc                                                                                                                                          >  1:4516961/1‑150 (MQ=255)

       ggTAGTCACCAGTACCCGGCTGAGCGGCCGGTGGTGTGGCATGAAAGTGTTCACGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccc                                                                                                                                          <  2:169218/150‑1 (MQ=255)

        gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccgat                                                                                                                                         <  2:3991147/150‑4 (MQ=255)

         tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                        >  2:4273992/1‑150 (MQ=255)

          tgTCACCAGTACCCTGCTGTTCGGCCTGTGGTGTGGCATGATTGTGTTCAAGCTCGGCCGCGATACCGCCGTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGTCATTCGAATCGGCGCTgcgccgcgccgattc                                                                                                                                       <  1:4746042/149‑6 (MQ=255)

          aGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCAGCGCCCCCCa                                                                                                                                       <  2:4143925/150‑1 (MQ=255)

            tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGGCGCGGCCCCCaga                                                                                                                                     >  1:2073079/1‑148 (MQ=255)

            tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccccac                                                                                                                                     >  1:3659519/1‑147 (MQ=255)

               cctgTACCCTGCTGATCGGCCTGTGGTGTGGCTTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGGTGACCGCCGCCGGCTGCGCATTCTGCGGCGCCGCGGCCTTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                                  <  2:70394/147‑1 (MQ=255)

               ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                                  <  2:2094784/150‑1 (MQ=255)

                cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGTCGCAGCGCCCCCCagcaaac                                                                                                                                 >  1:4482228/1‑144 (MQ=255)

                 aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAgc                                                                                                                                <  2:2001075/150‑1 (MQ=255)

                   tACCCTGCTGATCGGCCTGGGGTGTGGCATGAACGGGTTCACGCTGGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCGGGCTTTGGAATCGGGGCTGCGCAGCGCCCCCCATAAAAgcgc                                                                                                                              <  2:4056283/150‑1 (MQ=255)

                      ccTGCTGTTGGGCGTGTGTTGTGGCATTAATGTGTTCATGCTCGACGGCGATACCGCCTTGCTGACCGCGGCGGGCTGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATGGGCGCTGCGATGCGCCCCCCATAAAAGCGCCAt                                                                                                                           <  2:5340179/150‑1 (MQ=255)

                      ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAt                                                                                                                           <  2:4614799/150‑1 (MQ=255)

                      ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCCCCCCCCCAAAAAACCCCaa                                                                                                                           >  1:201077/1‑149 (MQ=255)

                        ttgtgATCGGCCTGTGGTGTGGCATTAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGAAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATgg                                                                                                                         <  2:5497350/147‑1 (MQ=255)

                        tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAGCGcccagg                                                                                                                         >  1:3121832/1‑146 (MQ=255)

                         gctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCCAAAAAGCGCCCTGGc                                                                                                                        >  1:779226/1‑150 (MQ=255)

                           tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCACCGCCCCCCCCAAAAACCCCAACGCAg                                                                                                                      >  1:2458629/1‑150 (MQ=255)

                            gATCGGCCTGTGGTGTGGCTTGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTATGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGt                                                                                                                     <  2:4205296/150‑1 (MQ=255)

                              tgggCCTGTGGTGTGGCATGAATGTGTTCACGCTCGACCGGGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCGGGCATTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCg                                                                                                                   <  2:3528310/148‑1 (MQ=255)

                                      tggtgtggCTTGAAAGTGTTCACGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                           <  2:5030859/150‑1 (MQ=255)

                                      tggtgtggCATGAATGTTTTCAAGCTCGACCGCGATACCGCCTTGCTGCCCGCCGCCGGCCGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCATTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                           <  2:150088/150‑1 (MQ=255)

                                       ggtgtggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAg                                                                                                          <  2:632627/150‑1 (MQ=255)

                                        gtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGt                                                                                                         >  2:4633069/1‑150 (MQ=255)

                                           tggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCTTGTTGGCATTCGAATCGGCGCTGCGAAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                      <  2:439563/150‑1 (MQ=255)

                                            ggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCACCGCCCCCCCAAAAAACGCCCAGGCCGTCGGAAGCGTAGttcgg                                                                                                     >  1:3361668/1‑146 (MQ=255)

                                               aTGAACGTGTTCAAGCTCGACCGCGATACCGCCTGGCTGCCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                  <  2:2850422/150‑1 (MQ=255)

                                               aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCGGTTc                                                                                                  >  2:2377334/1‑150 (MQ=255)

                                                 taatttgtTCACGCTCGACCGCGATCCCGCGTTGCTGACCGCCGCCGGCAGCGCCTTGTGCGGCGCCGCGGCGGTGTTGGATTTCGATTCGGGGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCgg                                                                                                <  1:3589673/145‑1 (MQ=255)

                                                 gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCCAGGCAGTCGGACGCGTAGTGGTGTTCgg                                                                                                >  1:4559664/1‑150 (MQ=255)

                                                  aaaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                               <  2:1611338/150‑1 (MQ=255)

                                                   aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCACCGCCCCCCAAAAAAACGCCAAGGCACTCGGAACCGGAGTGCTGTTAGGTa                                                                                              >  1:4486714/1‑150 (MQ=255)

                                                   aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCCAAAACCGCCCAGGCCGGCGGCAGCGGAGTGCGGTTCgggc                                                                                              >  1:2741748/1‑148 (MQ=255)

                                                     gtgtTCAAGCTCGACCGCGATACCGCCTTCCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCACCCCTAAAAGCGCCAAGGCACTCGGAAGCGAAGTTCTGTGAGGTAcc                                                                                            >  2:2940242/1‑150 (MQ=255)

                                                       gtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCGAAAAGCGCCCAGGCAGTCGGGAGCGGAGTGCTGTTGGGttgcat                                                                                          >  1:2886380/1‑145 (MQ=255)

                                                        ttCAAGCTCGACTGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCTAAAAGCGCCATGGCCGGCGGGAGCGTGGTGGTGGTCGGTACCtgg                                                                                         >  2:5485320/1‑148 (MQ=255)

                                                        ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCCCCCCCCCAAAAATCGCCCAGGCAGGAGGAATCCTAGTGCTGTTggggaacctg                                                                                         >  1:2949787/1‑143 (MQ=255)

                                                        ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCCCGCCCCCCAGAAAATCGCCCAGGCCGACGGCAGCCGAGTGGTGTTGGGtgacctg                                                                                         >  1:3402169/1‑144 (MQ=255)

                                                          cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGATTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                       <  1:3745685/150‑1 (MQ=255)

                                                          cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCTTGTc                                                                                       >  2:1883817/1‑150 (MQ=255)

                                                           aaGCTCGACCGCTTTACCGCCTTGCTGCCCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCACCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                                      <  1:1258546/150‑1 (MQ=255)

                                                           aaGCTCGACCGCTATACCGCCTTGCTGCCCGCCGCCGGCAGCGCATTCTGGGGCGCCGCCGCGTTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                                      <  1:4112995/150‑1 (MQ=255)

                                                           aaGCACGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTAGGCAGCGTAGTGCGGTTCGGTACCTTGTCg                                                                                      >  2:2395696/1‑150 (MQ=255)

                                                               tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCATAAAATCGCCATGGCAGGCGGAAGCGTAGTGGTGTTCGGTACCGTGTCGATGt                                                                                  >  1:2800122/1‑150 (MQ=255)

                                                                cGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGtt                                                                                 <  2:3133372/150‑1 (MQ=255)

                                                                     gcgATACCGCCTTGTTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCttt                                                                            >  2:2223784/1‑150 (MQ=255)

                                                                      cgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCTTTCGACGGGGCGCTGCCCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                           <  1:1287356/150‑1 (MQ=255)

                                                                      cgACACCGCCTGGGTGCCCGCCGCCGCCAGCGCCAGCTGCGGCGCCGCCGCGGTGGTGGATTTCGACGGGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                           <  1:1539502/150‑1 (MQ=255)

                                                                            cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCAGAAAAGCGCCAAGGCAGGCGGAAGCGTAGTGCTGGTCGGTACCTTGTCGGGGTTGCGTTAGCCACt                                                                     >  1:3176612/1‑150 (MQ=255)

                                                                              ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAGAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGGTGTTCGGTACCATGTCGAGGTTCCGTTAACCACAgg                                                                   >  1:1365106/1‑150 (MQ=255)

                                                                              ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCTGCGCCGCGCCCCCCCTACAAGCGCCATGGCAGGCGGAAGCGGAGTGCGGGTCGGTCCCGGGGCGCTGTGCCGTTACCCACAgg                                                                   >  2:3705439/1‑150 (MQ=255)

                                                                                ttGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACAGGCg                                                                 >  2:492966/1‑150 (MQ=255)

                                                                                 tggtGACGGCGGCCGGCTGCGCGATGTGCGGCGCGGCCGCGGTGCGGGATTTCGATTGGGGGCTGCGCGGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGa                                                                <  1:5495443/147‑1 (MQ=255)

                                                                                      aCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                           >  2:2113545/1‑150 (MQ=255)

                                                                                        cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGc                                                         >  2:2441790/1‑150 (MQ=255)

                                                                                         gccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGcc                                                        <  2:5207607/150‑1 (MQ=255)

                                                                                                gCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggc                                                 <  1:3703801/150‑1 (MQ=255)

                                                                                                 cAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggct                                                <  2:4496017/150‑1 (MQ=255)

                                                                                                  aGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCGGTTCGGTACCGTGTCGAAGTTCCTTTACCCACAGGCGATCAATGCCGgctggctg                                               >  2:2297956/1‑150 (MQ=255)

                                                                                                  aGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCCCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTAACCACAGGCGATCAATGCCGgctggctg                                               >  2:1145224/1‑150 (MQ=255)

                                                                                                    cgcCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCAGAAAAGCGCCCAGGCCGTCGGAAGCGTAGTGCTGGTCGGTACCTTGTCGGTGGTCGTTTAACCACCGGCGATCAATTCCGGCTGGCTGcc                                             >  1:5107924/1‑149 (MQ=255)

                                                                                                          gggcggcgCCGCCGCGGTGCGGGCATTCTTCGCGGCGCGGCGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                       <  1:1895245/148‑1 (MQ=255)

                                                                                                          cTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAACCCCCATGGCAGTCGGAAGCCTAGTGCTGTTCGGTACCTTGTCCATGTTACTTTACCCACAGGCGATCAATGCCCGCTGGCGGCGTCTGGa                                       >  2:3473635/1‑150 (MQ=255)

                                                                                                          cTGCGGCGCCGCCGCCGTGCTGGAATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                       <  2:1400573/150‑1 (MQ=255)

                                                                                                            gcggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATa                                     >  2:1767234/1‑150 (MQ=255)

                                                                                                             cggcgCCGCCGCGGTGGGGGCATGCCACGCGGCGCGGCGCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATAc                                    <  1:3856750/150‑1 (MQ=255)

                                                                                                              ggcgCCGCGGCCGTGCGGGCTTTCTAATCGGCGGGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATAcc                                   <  2:1616232/150‑1 (MQ=255)

                                                                                                                         cgTGCTGGCTTTCGAATCGGGGCTGCGCGGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                        <  2:2601290/150‑1 (MQ=255)

                                                                                                                          gTGCTGGCATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                       <  2:1861174/150‑1 (MQ=255)

                                                                                                                            gCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATACAAGCGCCATGGCAGTCGGACGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACGGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                     <  2:4119311/150‑1 (MQ=255)

                                                                                                                            gCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                     <  2:1484620/150‑1 (MQ=255)

                                                                                                                              tGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgc                   >  2:4591021/1‑150 (MQ=255)

                                                                                                                              tGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGggcggc                   >  2:455988/1‑150 (MQ=255)

                                                                                                                                   ttCTAAGCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg              <  1:5337593/150‑1 (MQ=255)

                                                                                                                                   ttCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg              >  2:2262632/1‑150 (MQ=255)

                                                                                                                                   tacgccgcGGCGCTGCGCCCCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg              <  1:2514484/143‑1 (MQ=255)

                                                                                                                                    tCGACTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg             <  1:5627443/150‑1 (MQ=255)

                                                                                                                                    tCGACTCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg             <  2:3189551/150‑1 (MQ=255)

                                                                                                                                    tCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg             >  2:3218234/1‑150 (MQ=255)

                                                                                                                                      gattgggcgcggcgcTGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc           <  2:702102/145‑1 (MQ=255)

                                                                                                                                      gAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTAGGTACCATTTCGATGTTACTTTACCCACAGGCGATCAAATCCGGCTGGCTGCAAATTGATACCAAGGGCGCCGGGCTGCTGCTGggc           >  1:827567/1‑150 (MQ=255)

                                                                                                                                      gAATCGGCGCGGCGCCGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc           <  2:4516961/150‑1 (MQ=255)

                                                                                                                                          cGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa       <  2:2545298/150‑1 (MQ=255)

                                                                                                                                           ggCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAc      >  2:1205125/1‑150 (MQ=255)

                                                                                                                                           ggCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAc      >  1:4620335/1‑150 (MQ=255)

                                                                                                                                            gcgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt     >  2:1509978/1‑150 (MQ=255)

                                                                                                                                              gcTGCGCAGCGCCCCCCATAACAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt   <  2:3567429/150‑1 (MQ=255)

                                                                                                                                              gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt   >  1:4720274/1‑150 (MQ=255)

                                                                                                                                               cTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc  >  2:2498781/1‑150 (MQ=255)

CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTC  >  NC_002947/4163315‑4163607

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