BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	363,169	A→G	R787R (CGT→CGC)	lacZ ←	beta‑D‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	363,169	0	A	G	81.7%	550.4 / 127.5	268	R787R (CGT→CGC)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand): ref base A (24/25); new base G (108/111); total (132/136)
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.03e-01

Only 100 of 270 total aligned reads displayed.

GGTAATGGCCCGCCGCCTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGC  >  NC_000913/363078‑363259

   aaTGGCCCGCCGCCTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCggg                                                                                        >  2:7343775/1‑93 (MQ=255)

      ggCCCGCCGCCTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGa                                                                                     <  2:4405868/93‑1 (MQ=255)

       gCCCGCCGCCTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGaa                                                                                    <  2:5569154/93‑1 (MQ=255)

        cccGCCGCCTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAAc                                                                                   >  1:3061625/1‑93 (MQ=255)

         ccgccgCCTTCCAGCGTTCGACCCAGGCGTTACGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACt                                                                                  <  2:2356879/93‑1 (MQ=255)

           gccgccTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGa                                                                                >  2:6965587/1‑93 (MQ=255)

           gccgccTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATACGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGa                                                                                <  1:542731/93‑1 (MQ=255)

            ccgccTTCCAGCGTTCGGCCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGAt                                                                               <  2:3697880/93‑1 (MQ=255)

            ccgccTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGAt                                                                               >  2:4732694/1‑93 (MQ=255)

             cgccTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATc                                                                              <  1:1604005/93‑1 (MQ=255)

               ccTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATcgc                                                                            <  1:5038804/93‑1 (MQ=255)

               ccTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATcgc                                                                            <  1:3784922/93‑1 (MQ=255)

               ccTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATcgc                                                                            <  2:3984164/93‑1 (MQ=255)

               ccTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATcgc                                                                            >  2:1035442/1‑93 (MQ=255)

                 ttCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCgcgc                                                                          >  2:3321226/1‑93 (MQ=255)

                 ttCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCgcgc                                                                          >  1:7143482/1‑93 (MQ=255)

                   ccAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAg                                                                        <  1:2432970/93‑1 (MQ=255)

                   ccAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAg                                                                        >  2:2822595/1‑93 (MQ=255)

                    cAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAgc                                                                       >  1:1084552/1‑93 (MQ=255)

                    cAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAgc                                                                       <  1:6596948/93‑1 (MQ=255)

                    cAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAgc                                                                       <  1:7410498/93‑1 (MQ=255)

                     aGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAgcg                                                                      >  2:2482725/1‑93 (MQ=255)

                      gCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAgcgg                                                                     <  1:2778472/93‑1 (MQ=255)

                      gCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAgcgg                                                                     <  1:3440336/93‑1 (MQ=255)

                         ttCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGt                                                                  >  1:2365526/1‑93 (MQ=255)

                           cGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTca                                                                <  1:4048605/93‑1 (MQ=255)

                           cGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTca                                                                <  1:4397099/93‑1 (MQ=255)

                            gACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTcag                                                               >  1:7077384/1‑93 (MQ=255)

                            gACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTcag                                                               <  2:2057351/93‑1 (MQ=255)

                            gACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTcag                                                               <  1:1943229/93‑1 (MQ=255)

                            gACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTcag                                                               <  2:6496022/93‑1 (MQ=255)

                             aCCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTcagc                                                              >  1:3276377/1‑93 (MQ=255)

                              cccAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTcagca                                                             >  2:4674430/1‑93 (MQ=255)

                              cccAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTcagca                                                             <  1:2170517/93‑1 (MQ=255)

                                cAGGCGTTAGGGTCAATGCGGGTCGTTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAgt                                                           <  1:5599511/93‑1 (MQ=255)

                                cAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAgt                                                           <  1:4786346/93‑1 (MQ=255)

                                cAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAgt                                                           >  1:6858021/1‑93 (MQ=255)

                                  ggCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAgttg                                                         >  2:5725976/1‑93 (MQ=255)

                                     gTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGttt                                                      >  1:292315/1‑93 (MQ=255)

                                     gTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGttt                                                      <  1:5921093/93‑1 (MQ=255)

                                       tAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGttttt                                                    >  2:1188594/1‑93 (MQ=255)

                                        aGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGtttttt                                                   <  1:94633/93‑1 (MQ=255)

                                         gggTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTa                                                  >  1:7285257/1‑93 (MQ=255)

                                         gggTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTa                                                  >  2:1508554/1‑93 (MQ=255)

                                           gTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATc                                                >  2:7031867/1‑93 (MQ=255)

                                            tCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCg                                               >  2:3486645/1‑93 (MQ=255)

                                            tCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCg                                               >  2:6102960/1‑93 (MQ=255)

                                               aTGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCa                                            >  1:2086691/1‑93 (MQ=255)

                                                tGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCaa                                           <  2:4433220/93‑1 (MQ=255)

                                                  cGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATc                                         >  1:1004361/1‑93 (MQ=255)

                                                   gggTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATcc                                        >  2:7110678/1‑93 (MQ=255)

                                                   gggTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATcc                                        <  1:3205232/93‑1 (MQ=255)

                                                    ggTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCca                                       >  2:762962/1‑93 (MQ=255)

                                                     gTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCcac                                      <  2:1598198/93‑1 (MQ=255)

                                                     gTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCcac                                      >  1:3494044/1‑93 (MQ=255)

                                                      tCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCcaca                                     <  1:776639/93‑1 (MQ=255)

                                                      tCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCcaca                                     <  2:6359436/93‑1 (MQ=255)

                                                        gCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATc                                   >  1:3658042/1‑93 (MQ=255)

                                                        gCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATc                                   >  2:4269390/1‑93 (MQ=255)

                                                        gCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATc                                   >  1:1245557/1‑93 (MQ=255)

                                                         cTTCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCt                                  <  2:2639350/93‑1 (MQ=255)

                                                          ttCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCtg                                 <  1:3027260/93‑1 (MQ=255)

                                                           tCACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCtgt                                >  1:6904536/1‑93 (MQ=255)

                                                            cACTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCtgtg                               >  2:2238212/1‑93 (MQ=255)

                                                             aCTTACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTga                              <  2:2445101/93‑1 (MQ=255)

                                                                tACGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTgaaag                           >  2:5051894/1‑93 (MQ=255)

                                                                  cGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTgaaagaa                         <  1:6545933/93‑1 (MQ=255)

                                                                  cGCCAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTgaaagaa                         <  2:5669450/93‑1 (MQ=255)

                                                                    ccAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGaaagaaag                       >  1:6380363/1‑93 (MQ=255)

                                                                    ccAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGaaagaaag                       <  1:3473533/93‑1 (MQ=255)

                                                                     cAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGc                      <  2:6251824/93‑1 (MQ=255)

                                                                     cAATGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGc                      >  2:1324971/1‑93 (MQ=255)

                                                                     cAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGc                      <  1:5789887/93‑1 (MQ=255)

                                                                      aaTGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGcc                     >  2:3367472/1‑93 (MQ=255)

                                                                      aaTGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGcc                     >  1:4968947/1‑93 (MQ=255)

                                                                       aTGTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCt                    >  2:1858330/1‑93 (MQ=255)

                                                                         gTCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGa                  <  1:6306813/93‑1 (MQ=255)

                                                                          tCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGAc                 >  1:565391/1‑93 (MQ=255)

                                                                          tCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGAc                 >  2:5860045/1‑93 (MQ=255)

                                                                          tCGTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGAc                 >  1:4256898/1‑93 (MQ=255)

                                                                          tCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGAc                 <  2:6700240/93‑1 (MQ=255)

                                                                            gTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTg               <  1:4671382/93‑1 (MQ=255)

                                                                            gTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTg               <  1:141502/93‑1 (MQ=255)

                                                                            gTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTg               <  1:2615910/93‑1 (MQ=255)

                                                                            gTTATCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTg               >  1:6165365/1‑93 (MQ=255)

                                                                                tCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCgg           >  1:4022350/1‑93 (MQ=255)

                                                                                tCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCgg           >  2:3805787/1‑93 (MQ=255)

                                                                                tCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCgg           <  2:1573673/93‑1 (MQ=255)

                                                                                tCCAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCgg           <  1:1400225/93‑1 (MQ=255)

                                                                                 ccAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGt          >  1:2847716/1‑93 (MQ=255)

                                                                                  cAGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGtt         <  2:3277154/93‑1 (MQ=255)

                                                                                  cAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGtt         >  1:6702995/1‑93 (MQ=255)

                                                                                   aGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTa        >  1:4518581/1‑93 (MQ=255)

                                                                                   aGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTa        <  1:5987821/93‑1 (MQ=255)

                                                                                   aGCGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTa        >  2:1285320/1‑93 (MQ=255)

                                                                                     cGGTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTaaa      <  2:3123082/93‑1 (MQ=255)

                                                                                      ggTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAt     >  1:4146838/1‑93 (MQ=255)

                                                                                      ggTGCGCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAt     <  1:5304495/93‑1 (MQ=255)

                                                                                       gTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAtt    <  1:1740252/93‑1 (MQ=255)

                                                                                         gCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGc  >  1:3912527/1‑93 (MQ=255)

GGTAATGGCCCGCCGCCTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGC  >  NC_000913/363078‑363259

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