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.2%	810.6 / NA	251	M1M (ATG→ATC) †	ysgA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); major base G (117/132); minor base T (0/2); total (117/134)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 268 total aligned reads displayed.

TAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGC‑ATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑CCTTATAAAAATGTTGCGCAATGTTAACTATA  >  NC_000913/4013625‑4013802

  aGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                                        >  2:1504768/1‑91 (MQ=255)

  aGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                                        >  1:581259/1‑91 (MQ=255)

  aGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                                        >  1:4142530/1‑91 (MQ=255)

  aGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCAGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                                        >  2:4701075/1‑91 (MQ=255)

    cAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcgatg                                                                                      <  2:903812/93‑5 (MQ=255)

            ttATCCCCTTGTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGt                                                                              <  2:3821006/93‑1 (MQ=255)

             tATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTg                                                                             >  2:213860/1‑93 (MQ=255)

             tATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTg                                                                             >  2:4790108/1‑93 (MQ=255)

             tATCCCCTAGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTg                                                                             >  1:6572122/1‑93 (MQ=255)

               tCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGCAGCGAGACGCGATGCAGCAGGTGCa                                                                           >  1:3176805/1‑93 (MQ=255)

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

                ccccTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaa                                                                          >  2:2786222/1‑93 (MQ=255)

                 cctttGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaaa                                                                         <  2:760829/90‑1 (MQ=255)

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

                    ttGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATcc                                                                      >  2:681663/1‑93 (MQ=255)

                    ttGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATcc                                                                      >  2:3113285/1‑93 (MQ=255)

                      gtgAAGGGATGGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAg                                                                    <  1:6258139/93‑1 (MQ=255)

                       tgAAGGAATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGa                                                                   <  1:2974874/93‑1 (MQ=255)

                        gaggggATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAt                                                                  <  1:4523283/90‑1 (MQ=255)

                        gAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAt                                                                  >  2:3965472/1‑93 (MQ=255)

                         aaGGGATAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAtt                                                                 <  2:1432365/93‑1 (MQ=255)

                          aGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttg                                                                >  2:5635626/1‑93 (MQ=255)

                           gggATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgt                                                               >  1:1563254/1‑93 (MQ=255)

                           gggATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgt                                                               >  2:6260525/1‑93 (MQ=255)

                             gATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgttg                                                             >  1:1224666/1‑93 (MQ=255)

                              gaggagGTGAGCCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgt                                                            <  2:1157098/90‑1 (MQ=255)

                              aTAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgt                                                            >  1:2815848/1‑93 (MQ=255)

                              aTAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgt                                                            <  2:576587/93‑1 (MQ=255)

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

                                 gagGTGAAGCCCGCCACAATCGCGTCGT‑CGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgctgtgttg                                                        >  1:2493576/1‑93 (MQ=255)

                                  gggTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgttg                                                        <  2:1089821/92‑1 (MQ=255)

                                   ggTGACGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTgttgtgttgt                                                       <  2:2734151/93‑1 (MQ=255)

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

                                       aaGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGcc                                                   <  1:703451/93‑1 (MQ=255)

                                        accccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGCGTTGTTGCCa                                                  <  2:4195285/91‑1 (MQ=255)

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

                                         gCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCAt                                                 >  2:4265511/1‑93 (MQ=255)

                                         gCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCAt                                                 <  2:5366824/93‑1 (MQ=255)

                                         ccccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCAt                                                 <  2:4051237/92‑1 (MQ=255)

                                         ccccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCAt                                                 <  2:1037102/92‑1 (MQ=255)

                                          cccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                                                >  1:107099/1‑93 (MQ=255)

                                          cccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                                                >  1:356783/1‑93 (MQ=255)

                                          cccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                                                >  1:5855756/1‑93 (MQ=255)

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

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

                                            cGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGt                                              <  2:3674501/93‑1 (MQ=255)

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

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

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

                                                acaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc                                          <  2:5344646/93‑1 (MQ=255)

                                                acaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc                                          <  2:6524927/93‑1 (MQ=255)

                                                 caATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtct                                         >  2:5612858/1‑93 (MQ=255)

                                                  aaTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGGCATGGTATtctc                                        >  1:2554774/1‑93 (MQ=255)

                                                  aaTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtctc                                        <  1:4972125/93‑1 (MQ=255)

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

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

                                                     cgcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGt                                     >  2:1959747/1‑93 (MQ=255)

                                                       cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑c                                 >  1:6276923/1‑93 (MQ=255)

                                                       cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑c                                 <  2:3266372/93‑1 (MQ=255)

                                                        gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑cc                                <  2:5914123/93‑1 (MQ=255)

                                                        gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑cc                                >  1:6378173/1‑93 (MQ=255)

                                                        gtcgtcCGGGGTCTGAACGATGGGCGAAGCGAGAGGCGATGCAGCAGGGGCAAATCCAGATTGTTGTGTTGTTGCCATGGGATTTTCCGG‑A‑cc                                >  1:101316/1‑93 (MQ=255)

                                                           gtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑ccctt                             <  1:644660/93‑4 (MQ=255)

                                                           gtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑ccctt                             >  1:4087056/1‑90 (MQ=255)

                                                            tcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑ccctta                            <  1:1894735/93‑5 (MQ=255)

                                                            tcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑ccctta                            <  1:348145/93‑5 (MQ=255)

                                                              cGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTGCGT‑ACCCTtata                           >  2:4326285/1‑93 (MQ=255)

                                                              cGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTtata                           >  1:543048/1‑93 (MQ=255)

                                                               ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATaa                          >  2:2967791/1‑93 (MQ=255)

                                                               ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATaa                          <  2:5427320/93‑1 (MQ=255)

                                                                gggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATaaa                         >  1:2302538/1‑93 (MQ=255)

                                                                gggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATaaa                         >  2:1812380/1‑93 (MQ=255)

                                                                  gTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATaaaaa                       >  1:3869899/1‑93 (MQ=255)

                                                                  gTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATaaaaa                       >  2:2499070/1‑93 (MQ=255)

                                                                      gAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGtt                   >  2:2350480/1‑93 (MQ=255)

                                                                      gAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCCTATAAAAATGtt                   <  2:3588743/93‑1 (MQ=255)

                                                                        aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTgc                 <  2:3912716/93‑1 (MQ=255)

                                                                        aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTgc                 >  2:3438882/1‑93 (MQ=255)

                                                                           aTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAATATGTTGCGCa              >  2:5819634/1‑93 (MQ=255)

                                                                            tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCGCCGT‑ACCCTTATAAAAATGTTGCGCaa             <  2:3361699/93‑1 (MQ=255)

                                                                             ggTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAAt            <  1:6347980/93‑1 (MQ=255)

                                                                             ggTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAAt            <  1:4251600/93‑1 (MQ=255)

                                                                              gTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATg           <  2:1684062/93‑1 (MQ=255)

                                                                                 gAAGCGAGAGGCGATGCAGCAGGTGCACATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTCC‑CCTTATAAAAATGTTGCGCAATGTTa        <  1:3752282/93‑1 (MQ=255)

                                                                                 gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTa        >  1:5109749/1‑93 (MQ=255)

                                                                                  aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTaa       <  1:2729266/93‑1 (MQ=255)

                                                                                  aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTaa       >  1:6063114/1‑93 (MQ=255)

                                                                                  aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTaa       >  2:698852/1‑93 (MQ=255)

                                                                                   aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAAc      >  1:6698738/1‑93 (MQ=255)

                                                                                   aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAAc      >  1:133010/1‑93 (MQ=255)

                                                                                   aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAAc      <  2:5157815/93‑1 (MQ=255)

                                                                                    gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAACt     <  2:1253521/93‑1 (MQ=255)

                                                                                    gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAACt     >  1:4669633/1‑93 (MQ=255)

                                                                                     cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAACta    <  2:6010213/93‑1 (MQ=255)

                                                                                     cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAACta    <  1:1627089/93‑1 (MQ=255)

                                                                                     cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAACta    <  1:3215560/93‑1 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAACtata  <  2:6130238/93‑1 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAACtata  <  2:1368762/93‑1 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAACtata  <  2:6684629/93‑1 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑ACCCTTATAAAAATGTTGCGCAATGTTAACtata  <  2:3329705/93‑1 (MQ=255)

TAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGC‑ATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGT‑A‑CCTTATAAAAATGTTGCGCAATGTTAACTATA  >  NC_000913/4013625‑4013802

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