BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,613,300	G→C	5.1%	intergenic (‑162/+26)	cyaB ← / ← waaA	cyclolysin secretion/processing ATP‑binding protein CyaB/3‑deoxy‑D‑manno‑octulosonic acid transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,613,300	0	G	C	5.1%	1001.6 / 30.9	296	intergenic (‑162/+26)	cyaB/waaA	cyclolysin secretion/processing ATP‑binding protein CyaB/3‑deoxy‑D‑manno‑octulosonic acid transferase
Reads supporting (aligned to +/- strand): ref base G (161/119); new base C (9/6); total (170/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 302 total aligned reads displayed.

CCCACAAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAGGCGGCGC  >  NC_002947/5613151‑5613449

cccACAAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCg                                                                                                                                                       <  2:1999202/150‑1 (MQ=255)

 ccACAAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCgg                                                                                                                                                      <  2:3833987/150‑1 (MQ=255)

  cacaAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCggg                                                                                                                                                     >  1:6611630/1‑150 (MQ=255)

   acaAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGt                                                                                                                                                    <  1:4374911/150‑1 (MQ=255)

   acaAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGt                                                                                                                                                    >  2:3100396/1‑150 (MQ=255)

   acaAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGt                                                                                                                                                    >  1:3004704/1‑150 (MQ=255)

       gCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCttt                                                                                                                                                <  1:342865/150‑1 (MQ=255)

         ttCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGa                                                                                                                                              >  2:1214038/1‑150 (MQ=255)

           ccGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGt                                                                                                                                            >  2:465572/1‑150 (MQ=255)

           ccGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGt                                                                                                                                            <  1:6834396/150‑1 (MQ=255)

           ccGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGt                                                                                                                                            <  2:6772958/150‑1 (MQ=255)

               cGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGat                                                                                                                                        >  1:2634938/1‑150 (MQ=255)

                 cAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGatat                                                                                                                                      <  1:3899192/150‑1 (MQ=255)

                     cGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCAA‑AGACCCGGTC‑TAGAGTGGATATGTct                                                                                                                                  <  2:3729817/150‑2 (MQ=38)

                     cGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCg                                                                                                                                  <  1:468146/150‑1 (MQ=255)

                        gtgACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGtt                                                                                                                               >  2:605999/1‑150 (MQ=255)

                          gACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTAt                                                                                                                             >  2:1781085/1‑150 (MQ=255)

                           acacTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATc                                                                                                                            <  1:5617982/150‑1 (MQ=255)

                            cacTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCa                                                                                                                           <  2:1897158/150‑1 (MQ=255)

                              cTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGc                                                                                                                         <  2:6303121/150‑1 (MQ=255)

                               ttCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCg                                                                                                                        <  2:6413650/150‑1 (MQ=255)

                                tCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCgg                                                                                                                       >  2:4084751/1‑150 (MQ=255)

                                 cTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGa                                                                                                                      <  1:3751176/150‑1 (MQ=255)

                                 cTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGa                                                                                                                      <  1:1022631/150‑1 (MQ=255)

                                   gACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATc                                                                                                                    <  2:2147294/150‑1 (MQ=255)

                                      ttttAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAg                                                                                                                 >  1:2333126/1‑150 (MQ=255)

                                      ttttAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAg                                                                                                                 <  1:6377501/150‑1 (MQ=255)

                                      ttttAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAg                                                                                                                 <  2:50607/150‑1 (MQ=255)

                                      ttttAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAg                                                                                                                 <  2:6673231/150‑1 (MQ=255)

                                         tAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGcc                                                                                                              >  1:3277773/1‑150 (MQ=255)

                                           aaGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGcccc                                                                                                            <  1:4727057/150‑1 (MQ=255)

                                             gTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCa                                                                                                          <  2:1485927/150‑1 (MQ=255)

                                              tCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAg                                                                                                         >  1:1146601/1‑150 (MQ=255)

                                              tCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAg                                                                                                         >  2:7391/1‑150 (MQ=255)

                                                 aTTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGAcc                                                                                                      <  1:5414541/150‑1 (MQ=255)

                                                  ttCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGCTCAAGGCCCCCAGAccc                                                                                                     >  2:6191559/1‑150 (MQ=255)

                                                  ttCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGAccc                                                                                                     >  1:1350143/1‑150 (MQ=255)

                                                    cGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTc                                                                                                   <  1:2009570/150‑1 (MQ=255)

                                                    cGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTc                                                                                                   <  1:1809633/150‑1 (MQ=255)

                                                      tAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGa                                                                                                 >  1:5780712/1‑150 (MQ=255)

                                                               tCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCt                                                                                        >  2:2331772/1‑150 (MQ=255)

                                                                 cTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGc                                                                                      >  2:4804555/1‑150 (MQ=255)

                                                                  tAAAGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCACTCGCTGCa                                                                                     >  2:5245693/1‑150 (MQ=255)

                                                                      ggCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGc                                                                                 >  2:2153741/1‑150 (MQ=255)

                                                                          agaAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGccccc                                                                             >  1:728004/1‑150 (MQ=255)

                                                                                 cAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTgg                                                                      <  1:5324134/150‑1 (MQ=255)

                                                                                    tgtAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGccc                                                                   <  1:5041624/150‑1 (MQ=255)

                                                                                       aGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCa                                                                <  2:6560557/150‑1 (MQ=255)

                                                                                         gAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTCTCAGCGGATCAAGGCCCCCAGACCCTCGAGCCATCGCTGCAAAGCCCCCTGCGTGGCCCGCagc                                                              >  2:5595281/1‑148 (MQ=255)

                                                                                             ttACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCg                                                          >  1:5998403/1‑150 (MQ=255)

                                                                                             ttACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCg                                                          >  2:2752622/1‑150 (MQ=255)

                                                                                               aCAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGcc                                                        >  1:4164438/1‑150 (MQ=255)

                                                                                                cAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGccc                                                       >  2:870955/1‑150 (MQ=255)

                                                                                                 aaCTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTGTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCt                                                      >  2:3006219/1‑150 (MQ=255)

                                                                                                       cGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGc                                                <  2:5779187/150‑1 (MQ=255)

                                                                                                              ggAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCAccc                                         <  1:944836/150‑1 (MQ=255)

                                                                                                                    aataaCCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATgcgc                                   >  2:1379899/1‑150 (MQ=255)

                                                                                                                       aaCCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTgg                                <  1:905044/150‑1 (MQ=255)

                                                                                                                         ccAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGc                              <  2:54795/150‑1 (MQ=255)

                                                                                                                         ccAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGc                              >  2:3139617/1‑150 (MQ=255)

                                                                                                                          cAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCg                             >  2:4021651/1‑150 (MQ=255)

                                                                                                                          cAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCg                             >  1:2743757/1‑150 (MQ=255)

                                                                                                                          cAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCg                             >  1:2897177/1‑150 (MQ=255)

                                                                                                                          cAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCg                             >  2:4254978/1‑150 (MQ=255)

                                                                                                                          cAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCg                             >  1:6758823/1‑150 (MQ=255)

                                                                                                                          cAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCg                             >  1:824909/1‑150 (MQ=255)

                                                                                                                          cAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGAGCAAGGCCCCCAGACCCTCGAGCAATCGCTGCCAAGCCCCCTGATTGGCCCGCATCACCGCCCGGCCCGCCTCACCCCTGCGCTGGGCg                             >  2:15633/1‑150 (MQ=255)

                                                                                                                             aaTCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTcc                          <  1:6676780/150‑1 (MQ=255)

                                                                                                                               tCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTg                        <  2:2488853/150‑1 (MQ=255)

                                                                                                                                 aGACATATCCACTC‑T‑AGACCGGTTCTGTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCAGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCg                      >  2:2631701/1‑150 (MQ=255)

                                                                                                                                 aGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCg                      >  1:4572259/1‑150 (MQ=255)

                                                                                                                                  gACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCgg                     <  1:605999/150‑1 (MQ=255)

                                                                                                                                   aCATATCCACTCAA‑AGACCCGGTC‑TAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGc                    >  1:6119255/1‑150 (MQ=38)

                                                                                                                                   aCATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTCATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGc                    >  1:3076288/1‑150 (MQ=255)

                                                                                                                                    cATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCa                   <  1:2153741/150‑1 (MQ=255)

                                                                                                                                    cATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCa                   <  2:5329053/150‑1 (MQ=255)

                                                                                                                                      tatCCACTC‑T‑AGACCGGGTCATTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCACTCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGt                 <  2:6855458/150‑1 (MQ=255)

                                                                                                                                       atCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATg                                      <  1:3583186/128‑1 (MQ=255)

                                                                                                                                       atCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGtt                >  2:6805440/1‑150 (MQ=255)

                                                                                                                                        tCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTc               >  2:565182/1‑150 (MQ=255)

                                                                                                                                        tCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTc               >  2:1432200/1‑150 (MQ=255)

                                                                                                                                        tCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTc               <  2:4667406/150‑1 (MQ=255)

                                                                                                                                        tCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTc               >  2:348368/1‑150 (MQ=255)

                                                                                                                                        tCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCACATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTc               <  1:5631965/150‑1 (MQ=255)

                                                                                                                                        tCCACTC‑AAAGACCCGGTC‑TAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTc               >  2:10929/1‑150 (MQ=255)

                                                                                                                                          cACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTcca             >  2:4883922/1‑150 (MQ=255)

                                                                                                                                           cactcaaaGACCCGGTC‑TAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTcca             >  1:3604731/8‑150 (MQ=255)

                                                                                                                                           aCTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTccac            <  1:4792322/150‑1 (MQ=255)

                                                                                                                                           aCTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTccac            >  1:2489438/1‑150 (MQ=255)

                                                                                                                                           aCTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTccac            >  2:3572515/1‑150 (MQ=255)

                                                                                                                                           aCTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTccac            >  2:3288335/1‑150 (MQ=255)

                                                                                                                                             tc‑t‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTccacca          >  1:6244903/1‑150 (MQ=255)

                                                                                                                                             tc‑t‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCAGCGGCAGTTccacca          >  1:6591949/1‑150 (MQ=255)

                                                                                                                                              c‑t‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAg         >  1:1627320/1‑150 (MQ=255)

                                                                                                                                                 aaGACCCGGTC‑TAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAggc       >  2:2936593/2‑150 (MQ=255)

                                                                                                                                                  aGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCCAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGGTCCACCAggc       <  2:5704180/150‑1 (MQ=255)

                                                                                                                                                      cGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCa                                        <  2:3630031/113‑1 (MQ=255)

                                                                                                                                                      cGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAGgcggcg   >  1:2596106/1‑150 (MQ=255)

                                                                                                                                                      cGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAGgcggcg   >  1:5499206/1‑150 (MQ=255)

                                                                                                                                                      cGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAGgcggcg   >  2:4466309/1‑150 (MQ=255)

CCCACAAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTC‑T‑AGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAGGCGGCGC  >  NC_002947/5613151‑5613449

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