BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,544,357	+G	coding (163/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,356	1	.	G	99.5%	652.5 / ‑3.8	203	coding (162/1314 nt)	gntT	gluconate permease
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (107/95); total (108/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

Only 100 of 226 total aligned reads displayed.

GGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑GC‑GGACG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACAA  >  NC_000913/3544264‑3544437

ggCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc                                                                                      <  1:5149697/93‑1 (MQ=255)

ggCTTCATCGCTCTCGTCCTCGTGGCACTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc                                                                                      <  2:5091981/93‑1 (MQ=255)

 gCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑g                                                                                    <  2:2349000/93‑1 (MQ=255)

 gCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑g                                                                                    >  2:2781427/1‑93 (MQ=255)

  cTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑gg                                                                                   <  1:5594216/93‑2 (MQ=255)

   ttCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑ggc                                                                                  >  1:4019901/1‑91 (MQ=255)

   ttCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑ggc                                                                                  <  1:6294048/93‑3 (MQ=255)

   ttCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGCGTc‑ggc                                                                                  >  1:2646322/1‑91 (MQ=255)

    tCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑ggcg                                                                                 <  2:629744/93‑4 (MQ=255)

    tCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑ggcg                                                                                 >  1:5778734/1‑90 (MQ=255)

                   tCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGcc                                                                   <  1:1076592/93‑1 (MQ=255)

                   tCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGcc                                                                   <  1:5981087/93‑1 (MQ=255)

                    cGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCt                                                                  <  1:1939479/93‑1 (MQ=255)

                     gTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCtt                                                                 <  2:1395636/93‑1 (MQ=255)

                     gTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCtt                                                                 <  2:2319203/93‑1 (MQ=255)

                      tGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTg                                                                >  1:4511343/1‑93 (MQ=255)

                       ggCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGc                                                               >  1:3613977/1‑93 (MQ=255)

                        gcgcTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGcc                                                              <  2:2712239/93‑1 (MQ=255)

                        gcgcTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGcc                                                              >  2:2639925/1‑93 (MQ=255)

                          gcTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCt                                                            <  1:3007706/93‑1 (MQ=255)

                           cTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTg                                                           >  2:3674251/1‑93 (MQ=255)

                           cTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTg                                                           <  2:1909449/93‑1 (MQ=255)

                           cTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTg                                                           >  1:2760924/1‑93 (MQ=255)

                            ttGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGa                                                          >  2:3793616/1‑93 (MQ=255)

                            ttGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGa                                                          >  1:4146632/1‑93 (MQ=255)

                            ttGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGa                                                          >  2:4733843/1‑93 (MQ=255)

                            ttGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGa                                                          <  1:464947/93‑1 (MQ=255)

                             tGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGATATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGAt                                                         >  2:22460/1‑93 (MQ=255)

                               cTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCa                                                       >  2:4118880/1‑93 (MQ=255)

                                tGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCAt                                                      >  2:2552077/1‑93 (MQ=255)

                                 gtgggATTAATGAAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATg                                                     <  2:560604/90‑1 (MQ=255)

                                 gTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATg                                                     >  2:3988264/1‑93 (MQ=255)

                                  ttGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATgg                                                    >  1:4696970/1‑93 (MQ=255)

                                  ttGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATgg                                                    >  1:5679139/1‑93 (MQ=255)

                                     gATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtt                                                 >  2:3084703/1‑93 (MQ=255)

                                     gATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtt                                                 <  2:1870636/93‑1 (MQ=255)

                                     gATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtt                                                 >  2:763949/1‑93 (MQ=255)

                                     gATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtt                                                 >  1:3294092/1‑93 (MQ=255)

                                      aTTANTNCAAGGAATGCCGCTGGATAAAGTNATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGttt                                                <  1:76558/93‑1 (MQ=255)

                                      aTTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGttt                                                >  1:248808/1‑93 (MQ=255)

                                      aTTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGttt                                                <  2:1600147/93‑1 (MQ=255)

                                       ttAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtttt                                               <  2:5759995/93‑1 (MQ=255)

                                       ttAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtttt                                               <  2:3527988/93‑1 (MQ=255)

                                       ttAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtttt                                               >  1:3131629/1‑93 (MQ=255)

                                       ttAATGCAAGGAATGCCCCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGtttt                                               <  1:4989191/93‑1 (MQ=255)

                                        tAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTg                                              >  1:1610063/1‑93 (MQ=255)

                                          aTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTggg                                            >  1:806758/1‑92 (MQ=255)

                                          aTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgc                                            >  2:3365553/1‑93 (MQ=255)

                                          aTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgc                                            >  1:942426/1‑93 (MQ=255)

                                            gCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgcgc                                          >  2:2688101/1‑93 (MQ=255)

                                            gCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgcgc                                          <  1:2039450/93‑1 (MQ=255)

                                              aaGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCaa                                        >  1:2406792/1‑93 (MQ=255)

                                              aaGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCaa                                        >  1:1826030/1‑93 (MQ=255)

                                              aaGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCaa                                        <  1:2792620/93‑1 (MQ=255)

                                               aGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAAt                                       <  2:2485219/93‑1 (MQ=255)

                                                 gAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGc                                     <  2:4972147/93‑1 (MQ=255)

                                                 gAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGc                                     <  2:4526932/93‑1 (MQ=255)

                                                  aaTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCt                                    >  1:2506440/1‑93 (MQ=255)

                                                  aaTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCt                                    <  2:3649302/93‑1 (MQ=255)

                                                   aTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                                   <  2:3293258/93‑1 (MQ=255)

                                                   aTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                                   <  1:5773647/93‑1 (MQ=255)

                                                   aTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                                   <  2:4153680/93‑1 (MQ=255)

                                                    tGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTgg                                  <  2:5722655/93‑1 (MQ=255)

                                                     gCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGTCCTGATCATGGGTTTTGGCGCAATGCTggg                                 >  1:143511/1‑93 (MQ=255)

                                                         cTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaa                             >  1:2289332/1‑93 (MQ=255)

                                                          tGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaaa                            <  1:1196779/93‑1 (MQ=255)

                                                           tgATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAAt                           <  1:6209459/92‑1 (MQ=255)

                                                           ggATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAAt                           <  2:4768755/93‑1 (MQ=255)

                                                           ggATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAAt                           >  1:5445833/1‑93 (MQ=255)

                                                             aTAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGc                         >  1:1335014/1‑93 (MQ=255)

                                                             aTAAAGTTATTGGCTCCATCAAAGCCGGGGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTCTTGGCGCAATGCTGGGCAAAATGc                         >  1:2908612/1‑93 (MQ=255)

                                                              tAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGGTTTGGCGCAATGCTGGGCAAAATGCt                        >  1:139838/1‑93 (MQ=255)

                                                               aaaGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTg                       <  2:2973931/93‑1 (MQ=255)

                                                               aaaGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTg                       >  1:5823525/1‑93 (MQ=255)

                                                               aaaGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTg                       >  2:1979545/1‑93 (MQ=255)

                                                                aaGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg                      <  2:1822689/93‑1 (MQ=255)

                                                                aaGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg                      >  2:1409198/1‑93 (MQ=255)

                                                                 aGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGc                     <  1:1000206/93‑1 (MQ=255)

                                                                  gTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCa                    >  1:113338/1‑93 (MQ=255)

                                                                  atgaTTGTCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATCCTGGCa                    <  2:3403654/90‑1 (MQ=255)

                                                                    tATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGa                  >  2:4413680/1‑93 (MQ=255)

                                                                     aTTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGAc                 >  2:3997967/1‑93 (MQ=255)

                                                                     aTTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGAc                 <  1:2567916/93‑1 (MQ=255)

                                                                      ttGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGGTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACt                >  2:603693/1‑93 (MQ=255)

                                                                       tGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTg               <  2:3423295/93‑1 (MQ=255)

                                                                        ggCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGc              <  2:3146112/93‑1 (MQ=255)

                                                                          cTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCgg            <  1:366119/93‑1 (MQ=255)

                                                                           tCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGt           >  2:4891162/1‑93 (MQ=255)

                                                                           tCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGt           >  2:1231956/1‑93 (MQ=255)

                                                                           tCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGt           <  2:5208980/93‑1 (MQ=255)

                                                                           tCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGt           <  1:5537422/93‑1 (MQ=255)

                                                                              aTCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTggg        >  2:5889886/1‑92 (MQ=255)

                                                                              aTCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGgc        >  2:4365805/1‑93 (MQ=255)

                                                                               tCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGgcg       >  1:6154799/1‑93 (MQ=255)

                                                                               tCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGgcg       >  1:4823524/1‑93 (MQ=255)

                                                                                 aaaGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGca     >  1:3550823/1‑93 (MQ=255)

                                                                                  aaGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATTGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGcac    >  2:5509030/1‑93 (MQ=255)

                                                                                   aGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGcaca   <  1:4347354/93‑1 (MQ=255)

                                                                                    gCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACaa  <  2:3846686/93‑1 (MQ=255)

                                                                                    gCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACaa  <  1:1775848/93‑1 (MQ=255)

GGCTTCATCGCTCTCGTCCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑GC‑GGACG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACAA  >  NC_000913/3544264‑3544437

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