BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,544,360	+G	coding (166/1314 nt)	gntT →	gluconate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,544,358	1	.	G	99.5%	703.2 / ‑3.7	194	coding (164/1314 nt)	gntT	gluconate permease
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (90/103); total (90/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 196 total aligned reads displayed.

TCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑G‑‑C‑GGACG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACAA  >  NC_000913/3544283‑3544437

tCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGcc                                                                   <  2:317085/93‑1 (MQ=255)

tCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGcc                                                                   >  1:1757879/1‑93 (MQ=255)

 cGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCt                                                                  <  2:990068/93‑1 (MQ=255)

 cGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCt                                                                  <  2:1176655/93‑1 (MQ=255)

  gTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCtt                                                                 >  2:1234695/1‑93 (MQ=255)

   tGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGATGCTCGGTAGCCTTg                                                                >  2:2631687/1‑93 (MQ=255)

    ggCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGc                                                               <  2:3146517/93‑1 (MQ=255)

    ggCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGc                                                               >  1:2568519/1‑93 (MQ=255)

       gcTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCt                                                            >  1:6145588/1‑93 (MQ=255)

       gcTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCt                                                            <  1:1982322/93‑1 (MQ=255)

        cTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTg                                                           >  1:6050876/1‑93 (MQ=255)

        cTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTg                                                           >  1:6391831/1‑93 (MQ=255)

        cTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTg                                                           <  2:2722308/93‑1 (MQ=255)

         ttGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGa                                                          >  1:4500436/1‑93 (MQ=255)

           gCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATc                                                        <  2:429641/93‑1 (MQ=255)

           gCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATc                                                        >  1:1422522/1‑93 (MQ=255)

           gCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATc                                                        >  2:1568719/1‑93 (MQ=255)

           gCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATc                                                        >  2:2963480/1‑93 (MQ=255)

            cTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATct                                                       >  2:3692683/1‑92 (MQ=255)

            cTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCa                                                       <  2:1364412/93‑1 (MQ=255)

            cTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCa                                                       >  2:1327937/1‑93 (MQ=255)

            cTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCa                                                       <  1:4133975/93‑1 (MQ=255)

             tGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCAt                                                      >  2:22570/1‑93 (MQ=255)

              gTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATg                                                     <  1:3960345/93‑1 (MQ=255)

              gTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATg                                                     <  2:4155502/93‑1 (MQ=255)

                tGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATggg                                                   <  2:2458899/93‑1 (MQ=255)

                tGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATggg                                                   <  1:4954845/93‑1 (MQ=255)

                 ggATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGt                                                  <  1:4276255/93‑1 (MQ=255)

                 ggATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGt                                                  >  2:5524799/1‑93 (MQ=255)

                 ggATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGt                                                  <  1:5958219/93‑1 (MQ=255)

                  gATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtt                                                 >  2:1143453/1‑93 (MQ=255)

                  gATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtt                                                 >  2:5061425/1‑93 (MQ=255)

                  gATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtt                                                 <  2:3405919/93‑1 (MQ=255)

                   aTTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGttt                                                <  1:1549352/93‑1 (MQ=255)

                    ttAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtttt                                               <  1:4902173/93‑1 (MQ=255)

                    ttAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGgttt                                               >  1:245279/1‑93 (MQ=255)

                     tAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGGTTTg                                              >  1:2643464/1‑93 (MQ=255)

                       aTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgc                                            <  1:5196773/93‑1 (MQ=255)

                        tGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgcg                                           >  1:1134039/1‑93 (MQ=255)

                           aaGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCaa                                        <  2:446941/93‑1 (MQ=255)

                            aGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAAt                                       <  2:4019118/93‑1 (MQ=255)

                            aGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAAt                                       <  1:2051785/93‑1 (MQ=255)

                            aGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAAt                                       >  2:3328427/1‑93 (MQ=255)

                            aGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAAt                                       <  2:5571123/93‑1 (MQ=255)

                            aGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAAt                                       <  1:67665/93‑1 (MQ=255)

                              gAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGc                                     >  2:5942645/1‑93 (MQ=255)

                                aTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                                   >  2:5301526/1‑93 (MQ=255)

                                aTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                                   <  1:4404572/93‑1 (MQ=255)

                                aTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                                   <  1:6258505/93‑1 (MQ=255)

                                aTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                                   <  2:5417538/93‑1 (MQ=255)

                                aTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGCGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                                   <  1:5639240/93‑1 (MQ=255)

                                 tGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTgg                                  <  1:6452548/93‑1 (MQ=255)

                                  gCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTggg                                 >  1:3533903/1‑93 (MQ=255)

                                   ccGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGc                                >  1:5310863/1‑93 (MQ=255)

                                   ccGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGc                                <  1:5564505/93‑1 (MQ=255)

                                    cGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCa                               <  1:5820960/93‑1 (MQ=255)

                                     gCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaa                              >  2:3255430/1‑93 (MQ=255)

                                      cTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaa                             <  1:3309386/93‑1 (MQ=255)

                                      cTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaa                             >  2:3471556/1‑93 (MQ=255)

                                       tGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaaa                            >  2:4090987/1‑93 (MQ=255)

                                       tGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaaa                            >  1:1900557/1‑93 (MQ=255)

                                        ggATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGGCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAAt                           >  1:631741/1‑93 (MQ=255)

                                        ggATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTTGGCAAAAt                           <  1:6502106/93‑1 (MQ=255)

                                        ggATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAAt                           <  2:6471741/93‑1 (MQ=255)

                                        ggATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAAt                           >  1:1158355/1‑93 (MQ=255)

                                          aTAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGc                         <  2:5715382/93‑1 (MQ=255)

                                           tAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCt                        >  1:5894980/1‑93 (MQ=255)

                                           tAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCt                        <  1:4595026/93‑1 (MQ=255)

                                            aaaGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTg                       >  2:1384112/1‑93 (MQ=255)

                                             aaGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg                      <  2:2414102/93‑1 (MQ=255)

                                             aaGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg                      >  2:193902/1‑93 (MQ=255)

                                              aGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGc                     <  1:761220/93‑1 (MQ=255)

                                              aGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGc                     <  2:3853344/93‑1 (MQ=255)

                                               gTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGcg                    >  2:6380640/1‑92 (MQ=255)

                                               gTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCa                    >  2:1773733/1‑93 (MQ=255)

                                                 tATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTTATCATGGGTTTGGGCGCAATGCTGGGCAAAATGCTGGCAGa                  <  1:3026918/93‑1 (MQ=255)

                                                 tATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGCTTTGGCGCAATGCTGGGCAAAATGCTGGCAGa                  >  2:832099/1‑93 (MQ=255)

                                                  aTTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGAc                 <  2:5323470/93‑1 (MQ=255)

                                                  aTTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGAc                 <  1:2490911/93‑1 (MQ=255)

                                                  aTTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGAc                 <  2:1666485/93‑1 (MQ=255)

                                                   ttGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACt                <  1:5092165/93‑1 (MQ=255)

                                                   ttGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACt                <  2:1017110/93‑1 (MQ=255)

                                                    tGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTg               >  2:1897430/1‑93 (MQ=255)

                                                      gCTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCg             >  2:5006057/1‑93 (MQ=255)

                                                       cTCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCgg            <  1:167480/93‑1 (MQ=255)

                                                        tCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGt           >  1:5508448/1‑93 (MQ=255)

                                                        tCCATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGt           >  1:2573785/1‑93 (MQ=255)

                                                         ccATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTg          <  1:3127504/93‑1 (MQ=255)

                                                          cATCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTgg         <  1:4385345/93‑1 (MQ=255)

                                                          cATCAAAGCCGGTGTC‑GGAG‑GGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTgg         <  2:5179049/93‑1 (MQ=255)

                                                           aTCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGgc        <  2:1636767/93‑1 (MQ=255)

                                                           aTCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGgc        <  2:6640734/93‑1 (MQ=255)

                                                            tCAAAGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGgcg       <  1:2906419/93‑1 (MQ=255)

                                                              aaaGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGca     >  1:1469318/1‑93 (MQ=255)

                                                              aaaGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGca     >  1:1332448/1‑93 (MQ=255)

                                                               aaGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGcac    >  1:2680876/1‑93 (MQ=255)

                                                               aaGCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTAGGCAAAATGCTGGCAGACTGCGGCGgcgccc    >  2:6296775/1‑91 (MQ=255)

                                                                 gCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACaa  <  2:4732363/93‑1 (MQ=255)

                                                                 gCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACaa  <  2:6019347/93‑1 (MQ=255)

                                                                 gCCGGTGTCGG‑‑CGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACaa  <  2:1744600/93‑1 (MQ=255)

TCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑G‑‑C‑GGACG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACAA  >  NC_000913/3544283‑3544437

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