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.6%	797.1 / ‑4.2	235	A55G (GCG→GGG)	gntT	gluconate permease
Reads supporting (aligned to +/- strand): ref base . (0/0); major base G (128/106); minor base T (1/0); total (129/106)
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.91e-01

Only 100 of 235 total aligned reads displayed.

TCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑GC‑GGACG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACAA  >  NC_000913/3544283‑3544437

tCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGcc                                                                   >  2:2000658/1‑93 (MQ=255)

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

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

  gTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCtt                                                                 >  2:2402863/1‑93 (MQ=255)

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

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

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

      cgcTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGccc                                                             >  2:5932981/1‑93 (MQ=255)

      cgcTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGccc                                                             >  1:5322454/1‑93 (MQ=255)

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

          tccTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGAt                                                         <  1:2429901/91‑1 (MQ=255)

          tGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGAt                                                         <  2:7385903/93‑1 (MQ=255)

          tGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGAt                                                         <  2:6530366/93‑1 (MQ=255)

            cTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCa                                                       <  1:6537662/93‑1 (MQ=255)

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

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

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

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

               ttGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATgg                                                    <  2:4095028/93‑1 (MQ=255)

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

                tGGATTAGTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATggg                                                   <  2:3799065/93‑1 (MQ=255)

                tGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATggg                                                   <  2:4510117/93‑1 (MQ=255)

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

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

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

                      aaTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTgg                                             >  1:1386699/1‑93 (MQ=255)

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

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

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

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

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

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

                            aGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAAt                                       >  1:6932687/1‑93 (MQ=255)

                            aGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAAt                                       >  1:3660834/1‑93 (MQ=255)

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

                              gAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGc                                     >  1:3556152/1‑93 (MQ=255)

                              gAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGGGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGGTTTGGCGCAATGc                                     >  1:7211585/1‑93 (MQ=255)

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

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

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

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

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

                                    cGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCa                               >  1:2484148/1‑93 (MQ=255)

                                    cGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCa                               >  2:1464693/1‑93 (MQ=255)

                                     gCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaa                              >  2:2600273/1‑93 (MQ=255)

                                     gCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaa                              >  2:5270797/1‑93 (MQ=255)

                                     gCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaa                              >  2:847447/1‑93 (MQ=255)

                                     gCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaa                              >  1:7490909/1‑93 (MQ=255)

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

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

                                       tGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaaa                            >  2:5810480/1‑93 (MQ=255)

                                       tGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaaa                            >  2:5163707/1‑93 (MQ=255)

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

                                         gATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATg                          >  2:397004/1‑93 (MQ=255)

                                         gATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATg                          >  2:7221785/1‑93 (MQ=255)

                                         gATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATg                          >  1:5241974/1‑93 (MQ=255)

                                         gATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATg                          <  2:3155060/93‑1 (MQ=255)

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

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

                                           tAAAGTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCt                        >  1:3493482/1‑93 (MQ=255)

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

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

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

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

                                               gTTATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCa                    <  2:3689817/93‑1 (MQ=255)

                                                ttATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAg                   >  1:3290429/1‑93 (MQ=255)

                                                 tATTGGCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGa                  <  1:71223/93‑1 (MQ=255)

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

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

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

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

                                                     ggCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGc              >  2:3889008/1‑93 (MQ=255)

                                                     ggCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGc              >  1:7296821/1‑93 (MQ=255)

                                                     ggCTCCATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGc              >  1:4705884/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                                                         ccATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTg          >  1:1485729/1‑93 (MQ=255)

                                                          cATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTgg         >  2:1142310/1‑93 (MQ=255)

                                                          cATCAAAGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTgg         <  1:3061413/93‑1 (MQ=255)

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

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

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

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

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

                                                              aaaGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGca     <  1:2221569/93‑1 (MQ=255)

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

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

                                                               aaGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGcac    >  2:5534417/1‑93 (MQ=255)

                                                               aaGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGcac    >  1:6021199/1‑93 (MQ=255)

                                                               aaGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGcac    >  2:2759627/1‑93 (MQ=255)

                                                               aaGCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGcac    >  1:906010/1‑93 (MQ=255)

                                                                 gCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACaa  >  1:3360039/1‑93 (MQ=255)

                                                                 gCCGGTGTCGGCGGGACGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACaa  >  1:1777295/1‑93 (MQ=255)

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

TCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑GC‑GGACG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCAGACTGCGGTGGCGCACAA  >  NC_000913/3544283‑3544437

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