BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,013,717	+G	coding (3/999 nt)	ysgA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,013,717	1	.	G	99.6%	801.6 / ‑3.8	224	coding (3/999 nt)	ysgA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (110/113); total (110/114)
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.31e-01

Only 100 of 235 total aligned reads displayed.

TAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGC‑ATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑CC‑TTATAAAAATGTTGCGCAATGTTAACTATA  >  NC_000913/4013625‑4013802

tAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGc                                                                                          >  1:3539591/1‑93 (MQ=255)

  aGCAGGCTTTTTATCCCCTTGGGAAGGGATAGAGGTAAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                                        <  1:1144476/93‑3 (MQ=255)

  aGCAGGCGTTTTATCCCTTTGTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                                        <  1:1873916/93‑3 (MQ=255)

      ggCATTTTATCCCTTGGGGAAGGGATAGAGGTGAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgca                                                                                    <  2:3693112/93‑1 (MQ=255)

      ggCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGGGAGAGGCGATgca                                                                                    >  2:1215734/1‑93 (MQ=255)

      ggCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgca                                                                                    >  1:4631928/1‑93 (MQ=255)

      ggCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgca                                                                                    >  2:3791954/1‑93 (MQ=255)

      ggCATGTTACCCCCTTGTGAAGGGATAGAGGTGAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgca                                                                                    <  1:2459472/93‑1 (MQ=255)

         ttGTTCTCCCTTGGAGAAGGGATAGACGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagca                                                                                 <  2:2113734/92‑1 (MQ=255)

          tGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGcagcag                                                                                >  1:5749522/1‑93 (MQ=255)

          tGTTACCCCTTTGTGAAGGGATAGAGGTGAGGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGcagcag                                                                                <  2:2970761/93‑1 (MQ=255)

           tttACCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAgg                                                                               <  1:4913878/92‑1 (MQ=255)

            ttATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAggg                                                                              >  1:4683587/1‑92 (MQ=255)

            ttATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAggg                                                                              >  1:2220674/1‑92 (MQ=255)

            ttATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAggg                                                                              >  2:4623930/1‑92 (MQ=255)

            ttATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGt                                                                              >  2:2882001/1‑93 (MQ=255)

              aTCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGc                                                                            >  1:5625371/1‑93 (MQ=255)

               tCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCa                                                                           >  1:4403030/1‑93 (MQ=255)

                ccccTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGCGGCGATGCAGCAGGTGCaa                                                                          >  2:1628863/1‑93 (MQ=255)

                  ccTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAAt                                                                        >  2:5246592/1‑93 (MQ=255)

                   tttGTAAAGGGAAAGGGGTGACGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATc                                                                       <  1:3634568/92‑1 (MQ=255)

                   cTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATc                                                                       >  2:3446388/1‑93 (MQ=255)

                   cTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATc                                                                       <  2:1843859/93‑1 (MQ=255)

                   cTTGTGAAGGGATAGAGGTAACGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATc                                                                       <  2:878001/93‑1 (MQ=255)

                     tgtgAAGGGAGAGGTGTGAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCa                                                                     <  1:6050650/93‑1 (MQ=255)

                       tgAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGa                                                                   >  1:5572591/1‑93 (MQ=255)

                        gAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCACATCCAGAt                                                                  >  2:3266891/1‑93 (MQ=255)

                        gAAGGGAGAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAt                                                                  <  1:5217518/93‑1 (MQ=255)

                         gaGGGAAAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAtt                                                                 <  2:4320769/92‑1 (MQ=255)

                             gATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgttg                                                             <  1:889529/93‑1 (MQ=255)

                              aTAGAGGTGACGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgt                                                            <  1:3864174/93‑1 (MQ=255)

                               tAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtg                                                           >  2:4530662/1‑93 (MQ=255)

                                 gagGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                                                         >  2:4248382/1‑93 (MQ=255)

                                 gagGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                                                         >  1:1331620/1‑93 (MQ=255)

                                  gggTGAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgttg                                                        <  1:176249/92‑1 (MQ=255)

                                  agGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgttg                                                        >  1:2245270/1‑93 (MQ=255)

                                    gTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTgttgtt                                                      <  1:292174/93‑1 (MQ=255)

                                     tGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGttgttg                                                     >  2:4065814/1‑93 (MQ=255)

                                      gAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGc                                                    >  2:2281885/1‑93 (MQ=255)

                                       aaGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGcc                                                   >  2:2520187/1‑93 (MQ=255)

                                        aGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCa                                                  >  1:5388244/1‑93 (MQ=255)

                                          cccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                                                <  1:4508411/93‑1 (MQ=255)

                                          cccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                                                <  1:2875780/93‑1 (MQ=255)

                                           ccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATgg                                               >  2:3634479/1‑93 (MQ=255)

                                           ccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATgg                                               <  1:4723416/93‑1 (MQ=255)

                                           ccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATgg                                               >  2:5663360/1‑93 (MQ=255)

                                            cGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGt                                              >  1:944145/1‑93 (MQ=255)

                                               cacaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAtt                                           >  2:4052053/1‑93 (MQ=255)

                                               cacaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAtt                                           >  1:2739907/1‑93 (MQ=255)

                                                acaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc                                          >  1:5089787/1‑93 (MQ=255)

                                                acaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc                                          >  1:1407728/1‑93 (MQ=255)

                                                  aaTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtctc                                        >  2:3228433/1‑93 (MQ=255)

                                                  aaTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtctc                                        >  1:3446248/1‑93 (MQ=255)

                                                    tCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCg                                      <  2:5093336/93‑1 (MQ=255)

                                                    tCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCg                                      <  2:96813/93‑1 (MQ=255)

                                                    tCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCg                                      <  2:6029825/93‑1 (MQ=255)

                                                      gcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                                    <  2:2344027/93‑1 (MQ=255)

                                                        gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑cc                                 >  1:1027680/1‑93 (MQ=255)

                                                        gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑cc                                 >  1:4118910/1‑93 (MQ=255)

                                                         tcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccc                                >  2:1558643/1‑92 (MQ=255)

                                                         tcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccc                                <  1:2608580/93‑2 (MQ=255)

                                                         tcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccc                                <  1:1565057/93‑2 (MQ=255)

                                                            tcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccctta                             >  2:1351798/1‑89 (MQ=255)

                                                               ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATaa                          >  2:4146837/1‑93 (MQ=255)

                                                               ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATaa                          >  2:966403/1‑93 (MQ=255)

                                                               ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATaa                          >  1:2283413/1‑93 (MQ=255)

                                                                 ggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATaaaa                        >  2:2430244/1‑93 (MQ=255)

                                                                 ggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATaaaa                        >  1:4218384/1‑93 (MQ=255)

                                                                 ggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATaaaa                        <  1:5533450/93‑1 (MQ=255)

                                                                  gTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATaaaaa                       >  2:4827911/1‑93 (MQ=255)

                                                                  gTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATaaaaa                       <  1:3624164/93‑1 (MQ=255)

                                                                  gTCGGAACGATGGTCGAAGCTAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATaaaaa                       <  1:296111/93‑1 (MQ=255)

                                                                   tCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAAt                      <  1:4341061/93‑1 (MQ=255)

                                                                   tCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAAt                      >  2:2615448/1‑93 (MQ=255)

                                                                     tGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGt                    <  2:678197/93‑1 (MQ=255)

                                                                     tGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATCCTCCGTACCC‑TTATAAAAATGt                    <  1:2049591/93‑1 (MQ=255)

                                                                      gAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGtt                   >  1:3773178/1‑93 (MQ=255)

                                                                        aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTgc                 >  1:5892590/1‑93 (MQ=255)

                                                                        aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTgc                 <  1:4967677/93‑1 (MQ=255)

                                                                          gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTgcgc               <  2:472333/93‑1 (MQ=255)

                                                                            tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCaa             <  2:189645/93‑1 (MQ=255)

                                                                               tCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTACTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGt          >  2:5946744/1‑93 (MQ=255)

                                                                                cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGtt         <  1:228936/93‑1 (MQ=255)

                                                                                  aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTaa       <  2:217755/93‑1 (MQ=255)

                                                                                    gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACt     >  2:2490678/1‑93 (MQ=255)

                                                                                    gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACt     >  2:2637487/1‑93 (MQ=255)

                                                                                    gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACt     >  1:4307962/1‑93 (MQ=255)

                                                                                    gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑CCTTTATAAAAATGTTGCGCAATGTTAACt     <  1:1492865/93‑1 (MQ=255)

                                                                                     cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACta    <  2:3060995/93‑1 (MQ=255)

                                                                                     cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACta    <  2:105036/93‑1 (MQ=255)

                                                                                     cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACta    >  2:865039/1‑93 (MQ=255)

                                                                                      gagaGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACtat   >  1:1844938/1‑93 (MQ=255)

                                                                                      gagaGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACtat   <  2:5011050/93‑1 (MQ=255)

                                                                                      gagaGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACtat   <  2:109252/93‑1 (MQ=255)

                                                                                      gagaGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACtat   <  1:615067/93‑1 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAGTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACtata  <  2:5491740/93‑1 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACtata  >  1:200849/1‑93 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACtata  >  2:6132155/1‑93 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACtata  <  1:3346419/93‑1 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCC‑TTATAAAAATGTTGCGCAATGTTAACtata  <  1:3791954/93‑1 (MQ=255)

TAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGC‑ATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑CC‑TTATAAAAATGTTGCGCAATGTTAACTATA  >  NC_000913/4013625‑4013802

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