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.8%	1809.8 / ‑5.3	526	M1M (ATG→ATC) †	ysgA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); major base G (281/244); minor base T (0/1); total (281/245)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 559 total aligned reads displayed.

TAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGC‑ATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑CCTTATAAAAATGTTGCGCAATGTTAACTATA  >  NC_000913/4013625‑4013802

tAAGCAGGCATTTTACCCCTTTGTAAGGGGATAGAGGTAAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGc                                                                                         <  2:4939465/93‑1 (MQ=255)

gaaGCAGGGATTTTTCCCCCTTGGAAAGGAAAAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGc                                                                                         <  1:5698097/92‑1 (MQ=255)

  aGCAGGCATGTTATCCCCTTGTGAAGGAATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                                       <  1:7097009/93‑3 (MQ=255)

  aGCAGGCATGTTATCCCCTGGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                                       <  1:6410836/93‑3 (MQ=255)

   gCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcgat                                                                                      <  1:5451431/93‑4 (MQ=255)

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

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

          tGTTATCACCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGcagcag                                                                               >  2:3128491/1‑93 (MQ=255)

              aTCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGc                                                                           <  2:6760091/93‑1 (MQ=255)

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

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

                 cccTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaaa                                                                        >  2:5895907/1‑93 (MQ=255)

                 cccTGGTGAGGGGATAGAGGTGAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaaa                                                                        <  2:6654004/93‑1 (MQ=255)

                  ctttGTGAGGGGAAAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAAt                                                                       <  2:1586288/91‑1 (MQ=255)

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

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

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

                    ttGTGAGGTGAGGGAGGTGAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATcc                                                                     <  1:4721074/93‑1 (MQ=255)

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

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

                     tgtgAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCa                                                                    >  1:5771632/1‑93 (MQ=255)

                       gaaaGGGAGAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGa                                                                  <  2:7267423/91‑1 (MQ=255)

                         aaGGGATGAGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAtt                                                                <  1:8064710/93‑1 (MQ=255)

                          aGGGATAGAGGTGACCCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttg                                                               <  2:35942/93‑1 (MQ=255)

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

                             gATAGAGTGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgttg                                                            <  2:8501684/93‑1 (MQ=255)

                              agagagGTGGAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgt                                                           <  2:7644799/91‑1 (MQ=255)

                              aTAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCACCAGGGGCAAATCCAGATtgttgt                                                           >  1:2450499/1‑93 (MQ=255)

                              aTAGAGGGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgt                                                           <  2:895577/93‑1 (MQ=255)

                               tAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtt                                                          >  2:1895893/1‑92 (MQ=255)

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

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

                                aggggTGACGCCCGCCACAAGCGCGTCGTCCGGGGTCTAAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATGCAGATtgttgtgt                                                         <  2:2081062/90‑1 (MQ=255)

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

                                 gagGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                                                        <  1:910577/93‑1 (MQ=255)

                                  agGTGAACACCTCCATAATCGCGTCGTCCTGTCTTTGAACGATGGTCGAAGCTAGAGTCTATGCAGCATGTGCAAATCCAGATtgttgtgttg                                                       <  1:729432/93‑1 (MQ=255)

                                  agGTAAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgttg                                                       <  2:6708713/93‑1 (MQ=255)

                                     tGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGGAAATCCAGATTGTTGTGttgttg                                                    >  2:5996609/1‑93 (MQ=255)

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

                                      gAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGc                                                   <  2:6421650/93‑1 (MQ=255)

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

                                        aGCCCCCCCCCATCTCGTCGTCCGGGGTCTGGACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCa                                                 <  1:8684829/93‑1 (MQ=255)

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

                                         gCCCGCCACAATCGCGTCGGCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCAt                                                <  2:7484229/93‑1 (MQ=255)

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

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

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

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

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

                                             gCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTa                                            <  1:8487975/93‑1 (MQ=255)

                                             gCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGAAAATCCAGATTGTTGTGTTGTTGCCATGGTa                                            <  2:8137444/93‑1 (MQ=255)

                                              ccACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAt                                           >  2:3365981/1‑93 (MQ=255)

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

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

                                                 caACCGCGTCGTCCGGGGTCTGAACGATGGTCTAAGCGAGCGGCGGTGCAGCAGGCGCAAACACAGATTGTTGTGTTGTTGCCATGGTAttgt                                        >  2:7947005/1‑91 (MQ=255)

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

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

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

                                                   aTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTcc                                      >  1:6599687/1‑93 (MQ=255)

                                                   aTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTcc                                      >  1:6814508/1‑93 (MQ=255)

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

                                                    tCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCg                                     >  1:6969251/1‑93 (MQ=255)

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

                                                       cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑c                                 >  2:4219179/1‑93 (MQ=255)

                                                       cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑c                                 <  2:698281/93‑1 (MQ=255)

                                                       cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑c                                 <  2:4801141/93‑1 (MQ=255)

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

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

                                                          cgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccct                              >  2:4023072/1‑91 (MQ=255)

                                                          cgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccct                              <  1:8267657/93‑3 (MQ=255)

                                                           gtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccctt                             >  2:566245/1‑90 (MQ=255)

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

                                                              cGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTtata                           >  1:7190781/1‑93 (MQ=255)

                                                               ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa                          >  1:101893/1‑93 (MQ=255)

                                                                   tCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt                      >  1:3263887/1‑93 (MQ=255)

                                                                    cTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATg                     <  1:1516726/93‑1 (MQ=255)

                                                                     tGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGt                    >  2:5201538/1‑93 (MQ=255)

                                                                      gAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGtt                   <  2:8739458/93‑1 (MQ=255)

                                                                      gAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGtt                   >  1:4920749/1‑93 (MQ=255)

                                                                       aaCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTg                  >  1:6046355/1‑93 (MQ=255)

                                                                        aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgc                 >  1:4266788/1‑93 (MQ=255)

                                                                        aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgc                 >  1:3927023/1‑93 (MQ=255)

                                                                         cGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgcg                <  2:5110206/93‑1 (MQ=255)

                                                                           aTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCa              <  2:258765/93‑1 (MQ=255)

                                                                           aTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCa              <  1:5125149/93‑1 (MQ=255)

                                                                             ggTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAAt            <  1:2838277/93‑1 (MQ=255)

                                                                                cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGtt         <  1:2250844/93‑1 (MQ=255)

                                                                                 gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTa        >  2:3014256/1‑93 (MQ=255)

                                                                                 gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTa        >  1:2186888/1‑93 (MQ=255)

                                                                                 gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTa        <  1:1868303/93‑1 (MQ=255)

                                                                                  aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTaa       >  1:2973160/1‑93 (MQ=255)

                                                                                  aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTaa       >  1:4333116/1‑93 (MQ=255)

                                                                                  aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTaa       <  2:8304626/93‑1 (MQ=255)

                                                                                   aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAAc      <  1:6133029/93‑1 (MQ=255)

                                                                                   aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAAc      <  1:4520949/93‑1 (MQ=255)

                                                                                    gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACt     >  1:2339304/1‑93 (MQ=255)

                                                                                    gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACt     <  1:2021820/93‑1 (MQ=255)

                                                                                     cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACta    <  1:5561050/93‑1 (MQ=255)

                                                                                      gagaGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACtat   <  1:3739025/93‑1 (MQ=255)

                                                                                       agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACtata  <  1:5204171/93‑1 (MQ=255)

TAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGC‑ATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑CCTTATAAAAATGTTGCGCAATGTTAACTATA  >  NC_000913/4013625‑4013802

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