BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,904,345	T→C	intergenic (‑150/+136)	bglG ← / ← phoU	BglG/PhoU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,904,345	0	T	C	98.3%	733.2 / 0.1	238	intergenic (‑150/+136)	bglG/phoU	BglG/PhoU
Reads supporting (aligned to +/- strand): ref base T (1/3); new base C (55/179); total (56/182)
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.76e-01

Only 100 of 243 total aligned reads displayed.

TCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTAT  >  NC_000913/3904254‑3904434

tCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTg                                                                                          <  1:6223692/93‑1 (MQ=255)

    gTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCaga                                                                                      <  2:904513/93‑1 (MQ=255)

        tGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGaaa                                                                                  <  1:4460594/93‑1 (MQ=255)

         gCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGaaaa                                                                                 <  2:2739688/93‑1 (MQ=255)

          ccTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAAt                                                                                <  1:5729130/93‑1 (MQ=255)

          ccTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAAt                                                                                <  2:6166725/93‑1 (MQ=255)

           cTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATa                                                                               <  1:4907750/93‑1 (MQ=255)

            tGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATaa                                                                              >  2:6078093/1‑92 (MQ=255)

             gCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGc                                                                             <  1:2770807/93‑1 (MQ=255)

               gAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCaa                                                                           <  1:2943986/93‑1 (MQ=255)

                aaTGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAAt                                                                          >  2:6690115/1‑93 (MQ=255)

                 aTGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATg                                                                         <  2:1240249/93‑1 (MQ=255)

                  tGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATgc                                                                        <  2:3898571/93‑1 (MQ=255)

                  tGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATgc                                                                        <  2:3217824/93‑1 (MQ=255)

                  tGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGAAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATgc                                                                        <  1:53110/93‑1 (MQ=255)

                    cAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATgcgc                                                                      <  2:536717/93‑1 (MQ=255)

                    cAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATgcgc                                                                      <  2:5101553/93‑1 (MQ=255)

                       tAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTAt                                                                   <  2:7186035/93‑1 (MQ=255)

                        aaCAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTAtt                                                                  >  2:6621471/1‑93 (MQ=255)

                         aCAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTg                                                                 <  1:6416418/93‑1 (MQ=255)

                          cAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGa                                                                <  2:6942465/93‑1 (MQ=255)

                          cAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGa                                                                >  1:840374/1‑93 (MQ=255)

                           aaTCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGAt                                                               >  2:2917374/1‑93 (MQ=255)

                            aTCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATa                                                              >  1:4700026/1‑93 (MQ=255)

                               cAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATaaaa                                                           >  1:3562158/1‑93 (MQ=255)

                                aGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATaaaaa                                                          <  2:4417312/93‑1 (MQ=255)

                                  tCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAata                                                        <  1:6580324/93‑1 (MQ=255)

                                  tCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAata                                                        >  2:6427983/1‑93 (MQ=255)

                                  tCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAata                                                        >  1:6763270/1‑93 (MQ=255)

                                  tCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAata                                                        >  2:6933986/1‑93 (MQ=255)

                                  tCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAata                                                        <  2:3324019/93‑1 (MQ=255)

                                  tCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAata                                                        >  1:142978/1‑93 (MQ=255)

                                    aTTTATTAATGGTTTTTATAACGAATATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATg                                                      <  2:1472837/93‑1 (MQ=255)

                                    aTTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATg                                                      >  2:6831558/1‑93 (MQ=255)

                                     tttATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGa                                                     <  1:2498867/93‑1 (MQ=255)

                                     tttATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGa                                                     >  2:1108453/1‑93 (MQ=255)

                                     tttATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGa                                                     >  1:365133/1‑93 (MQ=255)

                                      ttattaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGAc                                                    >  2:6422302/1‑93 (MQ=255)

                                      ttattaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGAc                                                    <  1:2540939/93‑1 (MQ=255)

                                      ttattaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATT‑ATAAAAATATGaac                                                   >  2:2537075/1‑91 (MQ=255)

                                        attaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCa                                                  <  2:7546109/93‑1 (MQ=255)

                                        attaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCa                                                  <  2:7118536/93‑1 (MQ=255)

                                        attaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCa                                                  <  1:1511583/93‑1 (MQ=255)

                                         ttaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCAt                                                 <  1:747031/93‑1 (MQ=255)

                                         ttaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCAt                                                 <  1:519378/93‑1 (MQ=255)

                                          taATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATg                                                <  2:1259033/93‑1 (MQ=255)

                                            aTGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCt                                              >  2:3870110/1‑93 (MQ=255)

                                             tGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTc                                             >  1:7206963/1‑93 (MQ=255)

                                             tGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTc                                             <  1:7809315/93‑1 (MQ=255)

                                                tttttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCa                                          <  1:4569544/93‑1 (MQ=255)

                                                tttttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCa                                          <  1:1649918/93‑1 (MQ=255)

                                                 ttttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAg                                         <  1:7719574/93‑1 (MQ=255)

                                                 ttttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAg                                         <  1:7115101/93‑1 (MQ=255)

                                                 ttttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAg                                         >  1:911548/1‑93 (MQ=255)

                                                  tttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGt                                        <  2:432585/93‑1 (MQ=255)

                                                  tttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGt                                        <  1:3104034/93‑1 (MQ=255)

                                                  tttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGt                                        <  2:1707286/93‑1 (MQ=255)

                                                   ttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGtt                                       <  1:4236899/93‑1 (MQ=255)

                                                   ttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGtt                                       <  2:7241136/93‑1 (MQ=255)

                                                    tataACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGtta                                      >  1:5146790/1‑93 (MQ=255)

                                                      taACGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGttatt                                    >  1:3469567/1‑93 (MQ=255)

                                                       aaCGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGttatta                                   <  2:492312/93‑1 (MQ=255)

                                                        aCGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTaa                                  <  2:1910366/93‑1 (MQ=255)

                                                         cGAACATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAAc                                 <  1:2094521/93‑1 (MQ=255)

                                                           aaCATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACtt                               <  2:3496938/93‑1 (MQ=255)

                                                           aaCATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACtt                               <  2:3793658/93‑1 (MQ=255)

                                                            aCATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACttt                              <  1:6792947/93‑1 (MQ=255)

                                                            aCATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACttt                              <  2:7444146/93‑1 (MQ=255)

                                                             cATCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTtg                             <  1:1845780/93‑1 (MQ=255)

                                                               tCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTtgtg                           <  2:403346/93‑1 (MQ=255)

                                                               tCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTtgtg                           <  1:3866619/93‑1 (MQ=255)

                                                               tCCAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTtgtg                           <  1:7398037/93‑1 (MQ=255)

                                                                 cAGGTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTa                         <  1:807524/93‑1 (MQ=255)

                                                                   ggTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAAt                       <  2:5739957/93‑1 (MQ=255)

                                                                    gTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAAtt                      <  2:5262008/93‑1 (MQ=255)

                                                                    gTTCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAAtt                      <  1:4827882/93‑1 (MQ=255)

                                                                     ttCGGAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAAttt                     <  2:1636466/93‑1 (MQ=255)

                                                                        ggAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAg                  >  2:1951789/1‑93 (MQ=255)

                                                                         gAAATTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAgg                 >  2:3983499/1‑93 (MQ=255)

                                                                          aaaTTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGa                <  1:1228853/93‑1 (MQ=255)

                                                                          aaaTTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGa                <  1:1248070/93‑1 (MQ=255)

                                                                          aaaTTAATTTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGa                <  1:1548507/93‑1 (MQ=255)

                                                                           aattaattTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGaa               <  2:3590806/93‑1 (MQ=255)

                                                                            attaattTAATTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAAt              >  2:6265574/1‑93 (MQ=255)

                                                                             ttaatttaatTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAAtt             <  1:5863992/93‑1 (MQ=255)

                                                                             ttaatttaatTGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAAtt             <  1:5669032/93‑1 (MQ=255)

                                                                              taatttaattGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAAttt            <  1:7285584/93‑1 (MQ=255)

                                                                                  ttaattGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATaa        <  2:256446/93‑1 (MQ=255)

                                                                                   taattGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATaaa       <  1:7295172/93‑1 (MQ=255)

                                                                                   taattGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATaaa       <  2:2349548/93‑1 (MQ=255)

                                                                                   taattGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATaaa       <  1:2896899/93‑1 (MQ=255)

                                                                                    aattGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAg      >  1:6190740/1‑93 (MQ=255)

                                                                                    aattGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAg      >  2:5908009/1‑93 (MQ=255)

                                                                                    aattGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAg      <  1:196169/93‑1 (MQ=255)

                                                                                      ttGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGtt    <  1:5053803/93‑1 (MQ=255)

                                                                                      ttGCGCGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGtt    <  1:3294794/93‑1 (MQ=255)

                                                                                        gcgcgcAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTtat  <  1:5571123/93‑1 (MQ=255)

                                                                                        gcgcgcAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTtat  >  2:4804456/1‑93 (MQ=255)

                                                                                        gcgcgcAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTtat  <  1:1024001/93‑1 (MQ=255)

                                                                                        gcgcgcAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTtat  <  1:761279/93‑1 (MQ=255)

TCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTAT  >  NC_000913/3904254‑3904434

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