BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,139,669	+C	coding (233/330 nt) coding (132/159 nt)	yijI → b4558 ←	hypothetical protein hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,139,667	1	.	C	99.6%	931.3 / ‑5.0	253	L77F (TTG→TTC) A45G (GCA→GGA)	yijI b4558	hypothetical protein hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); major base C (105/146); minor base A (1/0); total (106/147)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 289 total aligned reads displayed.

GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑CCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAT  >  NC_000913/4139575‑4139754

gcTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTg                                                                                          <  2:2911366/93‑1 (MQ=255)

  ttCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑cc                                                                                       >  1:5469836/1‑93 (MQ=255)

  ttCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑cc                                                                                       >  1:4192213/1‑93 (MQ=255)

   tCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑ccc                                                                                      >  2:3549519/1‑92 (MQ=255)

   tCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑ccc                                                                                      <  2:5027600/93‑2 (MQ=255)

    ccTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑ccct                                                                                     <  2:4477208/93‑3 (MQ=255)

     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑ccctt                                                                                    <  2:3240543/93‑4 (MQ=255)

     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑ccctt                                                                                    >  1:4773826/1‑90 (MQ=255)

     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑ccctt                                                                                    <  2:4651587/93‑4 (MQ=255)

      tGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑cccttg                                                                                   >  2:1307899/1‑89 (MQ=255)

      tGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑cccttg                                                                                   <  2:3054914/93‑5 (MQ=255)

      tGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATAAGCGCAAGGGCGCTTTTGGCGTTGACCTTg                                                                                    >  1:1228933/1‑93 (MQ=255)

      tGTTCGCAGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGGGCAAGGGCCCTTTTGGCGTTG‑cccttg                                                                                   <  1:4473640/93‑5 (MQ=255)

         tcgtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCgg                                                                                 >  2:4965994/1‑93 (MQ=255)

          cgtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGt                                                                                <  1:4322442/93‑1 (MQ=255)

           gtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTc                                                                               >  2:872850/1‑93 (MQ=255)

            tcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCa                                                                              <  1:4102999/93‑1 (MQ=255)

             cgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAg                                                                             <  1:3592662/93‑1 (MQ=255)

             cgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAg                                                                             <  1:1005526/93‑1 (MQ=255)

             cgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAg                                                                             <  1:1993685/93‑1 (MQ=255)

              ggCCCCGGAAATCAGTAACTGGCAGTCATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGt                                                                            <  1:5358554/93‑1 (MQ=255)

               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                           >  1:3190247/1‑93 (MQ=255)

               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                           <  1:3071338/93‑1 (MQ=255)

               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                           <  1:5886155/93‑1 (MQ=255)

               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                           <  1:514119/93‑1 (MQ=255)

                 cccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgcg                                                                         >  1:2270202/1‑93 (MQ=255)

                    ggAAATCAGTAACTGGCAGTTATCGCCCGTCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAg                                                                      >  2:1378495/1‑93 (MQ=255)

                    ggAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAg                                                                      <  2:3094942/93‑1 (MQ=255)

                    ggAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAg                                                                      >  2:2129238/1‑93 (MQ=255)

                     gAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAgg                                                                     >  1:4669774/1‑93 (MQ=255)

                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                  >  1:5280200/1‑93 (MQ=255)

                         tCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTAt                                                                 <  2:4138190/93‑1 (MQ=255)

                           aGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGc                                                               <  1:6320250/93‑1 (MQ=255)

                           aGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGc                                                               >  1:1630978/1‑93 (MQ=255)

                            gTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGcc                                                              >  1:5148665/1‑93 (MQ=255)

                               aCTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACt                                                           <  1:6450874/93‑1 (MQ=255)

                                cTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTc                                                          <  2:5656359/93‑1 (MQ=255)

                                cTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTc                                                          >  2:2707591/1‑93 (MQ=255)

                                 tGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCa                                                         <  2:1887354/93‑1 (MQ=255)

                                  ggCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCaa                                                        <  2:125123/93‑1 (MQ=255)

                                  ggCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCaa                                                        <  2:6518385/93‑1 (MQ=255)

                                   gCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAt                                                       >  1:1494048/1‑93 (MQ=255)

                                    cAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCCCTCAAtt                                                      <  1:680050/93‑1 (MQ=255)

                                      gTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtg                                                    >  2:4576/1‑93 (MQ=255)

                                      gTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtg                                                    <  1:6158183/93‑1 (MQ=255)

                                        tATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgtg                                                  >  2:4029116/1‑93 (MQ=255)

                                         aTCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtga                                                 <  1:1272825/93‑1 (MQ=255)

                                            gcccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgatga                                              >  1:1119423/1‑93 (MQ=255)

                                            gcccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgatga                                              <  1:2966890/93‑1 (MQ=255)

                                             cccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGTAGGTTATGCCACTCAATTTGTGATGaa                                             <  1:895653/93‑1 (MQ=255)

                                             cccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaa                                             <  1:6551576/93‑1 (MQ=255)

                                               cgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaa                                           >  2:5057866/1‑93 (MQ=255)

                                               cgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaa                                           <  1:2508558/93‑1 (MQ=255)

                                                gccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaaa                                          <  1:6061882/93‑1 (MQ=255)

                                                gccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaaa                                          <  1:758892/93‑1 (MQ=255)

                                                gccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaaa                                          <  1:2861899/93‑1 (MQ=255)

                                                  cAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGTCACTCAATTTGTGATGAAAAAGt                                        >  1:4387734/1‑93 (MQ=255)

                                                  cAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGt                                        >  2:2405758/1‑93 (MQ=255)

                                                  cAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGt                                        >  2:3435722/1‑93 (MQ=255)

                                                  cAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGt                                        >  2:1690557/1‑93 (MQ=255)

                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                      >  1:6206387/1‑93 (MQ=255)

                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                      >  1:3747757/1‑93 (MQ=255)

                                                     cGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAtt                                     <  2:557512/93‑1 (MQ=255)

                                                      gTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTa                                    <  1:4962996/93‑1 (MQ=255)

                                                      gTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTa                                    <  1:5852411/93‑1 (MQ=255)

                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                   >  1:3618059/1‑93 (MQ=255)

                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                   <  1:3546328/93‑1 (MQ=255)

                                                       tATCGGCGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCACTTTGTGATGAAAAAGTATTac                                   >  1:2216999/1‑93 (MQ=255)

                                                        aTCGGTGCCAAGCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTaca                                  <  2:4212006/93‑1 (MQ=255)

                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 <  2:544308/93‑1 (MQ=255)

                                                          cGGTGTCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTAcaca                                <  2:2567542/93‑1 (MQ=255)

                                                          cGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTAcaca                                <  1:3808264/93‑1 (MQ=255)

                                                           ggTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAg                               <  1:4224027/93‑1 (MQ=255)

                                                           ggTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAg                               <  2:6233900/93‑1 (MQ=255)

                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGc                              <  1:3359228/93‑1 (MQ=255)

                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGc                              <  2:1133791/93‑1 (MQ=255)

                                                              gCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGt                            <  2:2270740/93‑1 (MQ=255)

                                                              gCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGt                            <  2:6131874/93‑1 (MQ=255)

                                                               ccAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGtt                           >  2:5537384/1‑93 (MQ=255)

                                                               ccAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGtt                           >  2:2999442/1‑93 (MQ=255)

                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                         <  1:6414146/93‑1 (MQ=255)

                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                         <  1:3898445/93‑1 (MQ=255)

                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                         >  1:5762033/1‑93 (MQ=255)

                                                                  aTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCa                        >  2:5721530/1‑93 (MQ=255)

                                                                  aTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCa                        >  2:220180/1‑93 (MQ=255)

                                                                   tCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCaa                       >  2:2362762/1‑93 (MQ=255)

                                                                    cAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAAc                      >  1:2774444/1‑93 (MQ=255)

                                                                      gcgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGt                    <  1:2549069/93‑1 (MQ=255)

                                                                       cgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTg                   >  2:5181161/1‑93 (MQ=255)

                                                                         cAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGc                 >  2:4780420/1‑93 (MQ=255)

                                                                          aaGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCt                <  2:590639/93‑1 (MQ=255)

                                                                          aaGGGCGCTTTTGGCGGTGCCCTTGCGGGCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACCCAGCGGTTCAACGTGGCt                <  1:2641237/93‑1 (MQ=255)

                                                                           aGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCtt               <  2:1943028/93‑1 (MQ=255)

                                                                           aGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCtt               >  2:1991381/1‑93 (MQ=255)

                                                                             ggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGc             <  2:2775905/93‑1 (MQ=255)

                                                                             ggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGc             <  1:1795033/93‑1 (MQ=255)

                                                                                gcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCgg          >  2:712181/1‑93 (MQ=255)

                                                                                 cTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGa         <  1:3488606/93‑1 (MQ=255)

                                                                                 cTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGa         <  1:1842588/93‑1 (MQ=255)

                                                                                        cGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAt  <  2:5341642/93‑1 (MQ=255)

GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTG‑CCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAT  >  NC_000913/4139575‑4139754

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