BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,954,264	A→T	D881E (GAT→GAA)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,954,264	0	A	T	99.7%	1148.9 / ‑4.1	319	D881E (GAT→GAA)	ptr	protease III
Reads supporting (aligned to +/- strand): ref base A (1/0); new base T (152/166); total (153/166)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 326 total aligned reads displayed.

CTGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGATCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTGCGC  >  NC_000913/2954174‑2954353

 tGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCg                                                                                        >  1:6857025/1‑93 (MQ=255)

  gATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGa                                                                                       >  2:2654875/1‑93 (MQ=255)

  gATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGa                                                                                       >  2:5448501/1‑93 (MQ=255)

  gATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGa                                                                                       >  1:4887751/1‑93 (MQ=255)

   gtGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGaa                                                                                      <  2:5997581/92‑1 (MQ=255)

     ggAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAAt                                                                                    <  2:3661061/93‑1 (MQ=255)

      gaagaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATc                                                                                   <  1:1083377/93‑1 (MQ=255)

      gaagaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATc                                                                                   >  2:5674859/1‑93 (MQ=255)

      gaagaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATc                                                                                   <  2:7662463/93‑1 (MQ=255)

        agaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCtt                                                                                 <  2:4738421/93‑1 (MQ=255)

        agaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCtt                                                                                 <  2:6136309/93‑1 (MQ=255)

         gaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCttt                                                                                <  2:3277082/93‑1 (MQ=255)

          aaaTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTa                                                                               <  1:6496793/93‑1 (MQ=255)

           aaTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTAc                                                                              <  2:4803424/93‑1 (MQ=255)

            aTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACt                                                                             <  2:5631142/93‑1 (MQ=255)

             tCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGGCGCGATTTTATCACGCGAATCGAAGCGCATATGGCCTCGTTCGAAATGTTTACtt                                                                            >  2:320749/1‑93 (MQ=255)

              cAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTa                                                                           <  1:162060/93‑1 (MQ=255)

               aGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCTAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTaa                                                                          <  1:3909591/93‑1 (MQ=255)

               aGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTac                                                                          >  2:3398621/1‑92 (MQ=255)

               aGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTaa                                                                          <  2:5696743/93‑1 (MQ=255)

                gCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAAc                                                                         <  2:5279786/93‑1 (MQ=255)

                gCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAAc                                                                         <  1:5061711/93‑1 (MQ=255)

                 cAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACt                                                                        <  2:1341608/93‑1 (MQ=255)

                   aGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTc                                                                      >  2:3833043/1‑93 (MQ=255)

                   aGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTc                                                                      <  1:2389128/93‑1 (MQ=255)

                   aGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTc                                                                      >  2:426502/1‑93 (MQ=255)

                   aGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTc                                                                      >  1:7798364/1‑93 (MQ=255)

                     tttttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGa                                                                    >  2:5077225/1‑93 (MQ=255)

                     tttttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGa                                                                    <  2:1534513/93‑1 (MQ=255)

                        ttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGc                                                                 <  1:72519/93‑1 (MQ=255)

                        ttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGc                                                                 <  1:3406580/93‑1 (MQ=255)

                            ggcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttct                                                             <  2:3005279/93‑1 (MQ=255)

                              cgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCttcttc                                                           <  1:2657462/93‑1 (MQ=255)

                              cgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCttcttc                                                           <  2:208739/93‑1 (MQ=255)

                              cgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCttcttc                                                           <  1:3084506/93‑1 (MQ=255)

                                tCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGc                                                         >  1:3747299/1‑93 (MQ=255)

                                  agcagTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCg                                                       >  2:183413/1‑93 (MQ=255)

                                   gcagTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGa                                                      <  2:1436384/93‑1 (MQ=255)

                                   gcagTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGa                                                      <  2:6522448/93‑1 (MQ=255)

                                   gcagTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGa                                                      >  1:2878696/1‑93 (MQ=255)

                                   gcagTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGa                                                      >  2:2926225/1‑93 (MQ=255)

                                    cagTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAg                                                     <  1:7416254/93‑1 (MQ=255)

                                     agTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGc                                                    <  1:1813225/93‑1 (MQ=255)

                                     agTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGc                                                    >  1:7385950/1‑93 (MQ=255)

                                     agTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGc                                                    >  2:1038110/1‑93 (MQ=255)

                                       ttttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGCAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGt                                                  >  2:4940387/1‑93 (MQ=255)

                                       ttttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGt                                                  >  1:526142/1‑93 (MQ=255)

                                       ttttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGt                                                  >  2:7348654/1‑93 (MQ=255)

                                        tttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGtt                                                 <  2:1724098/93‑1 (MQ=255)

                                        tttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGtt                                                 <  1:3474578/93‑1 (MQ=255)

                                        tttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGtt                                                 <  2:6652045/93‑1 (MQ=255)

                                        tttATCTGGGCCACGATCTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGtt                                                 <  1:6959882/93‑1 (MQ=255)

                                          tATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTg                                               <  2:6558891/93‑1 (MQ=255)

                                            tCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCg                                             >  2:4514903/1‑93 (MQ=255)

                                            tCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCg                                             >  2:7122548/1‑93 (MQ=255)

                                            tCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCg                                             >  1:4540391/1‑93 (MQ=255)

                                            tCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCg                                             >  2:6923047/1‑93 (MQ=255)

                                             cTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCgg                                            <  2:7528703/93‑1 (MQ=255)

                                                 gCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGcc                                        <  2:7810808/93‑1 (MQ=255)

                                                   cACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTg                                      >  2:7909671/1‑93 (MQ=255)

                                                    aCGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgc                                     >  1:1742680/1‑93 (MQ=255)

                                                     cGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgca                                    >  1:3279135/1‑93 (MQ=255)

                                                      gATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcag                                   >  1:2207603/1‑93 (MQ=255)

                                                      gATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcag                                   >  1:2510075/1‑93 (MQ=255)

                                                       aTTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcagc                                  <  1:3557165/93‑1 (MQ=255)

                                                        ttttATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcagca                                 >  1:424965/1‑93 (MQ=255)

                                                         tttATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCAt                                >  2:6954539/1‑93 (MQ=255)

                                                         tttATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCAt                                >  1:6587856/1‑93 (MQ=255)

                                                           tATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCt                              >  2:309207/1‑93 (MQ=255)

                                                           tATCACGAGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCt                              <  1:7544873/93‑1 (MQ=255)

                                                            aTCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTg                             <  2:6794092/93‑1 (MQ=255)

                                                                 gcgAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTaa                        <  1:7916850/93‑1 (MQ=255)

                                                                 gcgAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGATCCCTGCAGCATCTGGGTaa                        <  1:6667129/93‑1 (MQ=255)

                                                                  cgAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAAt                       <  1:1109085/93‑1 (MQ=255)

                                                                  cgAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAAt                       <  1:1573965/93‑1 (MQ=255)

                                                                    aaTCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTa                     >  1:5597872/1‑93 (MQ=255)

                                                                     aTCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTAc                    >  2:5887678/1‑93 (MQ=255)

                                                                       cGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCg                  <  2:2878503/93‑1 (MQ=255)

                                                                        gAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGc                 >  1:1284710/1‑93 (MQ=255)

                                                                        gAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGc                 >  1:3514164/1‑93 (MQ=255)

                                                                         aaGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTAGCGcc                >  2:2210482/1‑93 (MQ=255)

                                                                         aaGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTAGGTAATTACCGcc                <  2:4437979/93‑1 (MQ=255)

                                                                          aGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCct               >  1:538035/1‑93 (MQ=255)

                                                                           gcgcATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctg              >  2:5794827/1‑93 (MQ=255)

                                                                           gcgcATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctg              <  1:6057875/93‑1 (MQ=255)

                                                                             gcATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctgct            >  2:6608995/1‑93 (MQ=255)

                                                                              cATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctgctg           <  1:1167476/93‑1 (MQ=255)

                                                                              cATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctgctg           >  2:4071865/1‑93 (MQ=255)

                                                                                tatTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGa         >  1:7646837/1‑93 (MQ=255)

                                                                                tatTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGa         <  2:3550104/93‑1 (MQ=255)

                                                                                tatTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGa         >  1:1968213/1‑93 (MQ=255)

                                                                                 atTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCCGCATCTGGGTAATTACCCCCTGCTGGAt        >  1:5921784/1‑93 (MQ=255)

                                                                                 atTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGAt        <  1:3641412/93‑1 (MQ=255)

                                                                                 atTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGAt        >  2:6822473/1‑93 (MQ=255)

                                                                                  ttGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGAtt       >  2:235128/1‑93 (MQ=255)

                                                                                   tGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGAttt      <  1:7329352/93‑1 (MQ=255)

                                                                                   tGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGAttt      <  1:5437698/93‑1 (MQ=255)

                                                                                   tGCCGCGTTCGAAAGCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGGGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGAttt      <  2:2528318/93‑1 (MQ=255)

                                                                                      cgcgTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTgcg   >  2:2024193/1‑93 (MQ=255)

                                                                                       gcgTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTgcgc  >  2:744327/1‑93 (MQ=255)

CTGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGATCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTGCGC  >  NC_000913/2954174‑2954353

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