BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,007,730	T→C	W316R (TGG→CGG)	ygeY →	putative deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,007,730	0	T	C	83.4%	691.1 / 135.7	313	W316R (TGG→CGG)	ygeY	putative deacetylase
Reads supporting (aligned to +/- strand): ref base T (26/26); new base C (131/130); total (157/156)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

Only 100 of 313 total aligned reads displayed.

ACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTGGAC  >  NC_000913/3007639‑3007821

aCAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTg                                                                                            >  1:4994349/1‑93 (MQ=255)

    aaaGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGaaa                                                                                        <  2:6066981/93‑1 (MQ=255)

      aGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGt                                                                                      <  2:714861/93‑1 (MQ=255)

      aGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGt                                                                                      <  1:4610227/93‑1 (MQ=255)

       gCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTg                                                                                     <  2:763098/93‑1 (MQ=255)

       gCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTg                                                                                     <  2:2943594/93‑1 (MQ=255)

         tAACGCTGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGga                                                                                   >  2:4725860/1‑93 (MQ=255)

         tAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGga                                                                                   >  2:2674223/1‑93 (MQ=255)

         tAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGga                                                                                   >  1:3858298/1‑93 (MQ=255)

         tAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGga                                                                                   >  1:2472451/1‑93 (MQ=255)

          aaCGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGgaa                                                                                  >  2:4381596/1‑93 (MQ=255)

           aCGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGgaag                                                                                 <  1:795562/93‑1 (MQ=255)

              cgGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGaagaag                                                                              <  1:1120120/93‑1 (MQ=255)

               ggTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAagaaga                                                                             >  1:5812560/1‑93 (MQ=255)

               ggTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAagaaga                                                                             >  1:6416065/1‑93 (MQ=255)

               ggTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAagaaga                                                                             <  1:2827795/93‑1 (MQ=255)

                gttgttTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGAt                                                                            >  1:4519339/1‑93 (MQ=255)

                gttgttTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGAt                                                                            <  2:5199271/93‑1 (MQ=255)

                  tgttTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCa                                                                          <  1:3769180/93‑1 (MQ=255)

                    tttCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACt                                                                        >  2:1653071/1‑93 (MQ=255)

                     ttCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACtt                                                                       <  2:2597159/93‑1 (MQ=255)

                     ttCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACtt                                                                       <  1:2561283/93‑1 (MQ=255)

                     ttCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACtt                                                                       <  1:6697348/93‑1 (MQ=255)

                      tCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTc                                                                      <  1:2706523/93‑1 (MQ=255)

                        tATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCAc                                                                    <  2:5479446/93‑1 (MQ=255)

                        tATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCAc                                                                    >  1:5306188/1‑93 (MQ=255)

                         aTGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCAcc                                                                   <  1:4845822/93‑1 (MQ=255)

                           gTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGt                                                                 >  1:1886344/1‑93 (MQ=255)

                            tACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGtt                                                                <  1:6976005/93‑1 (MQ=255)

                              cAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTaa                                                              >  1:4354684/1‑93 (MQ=255)

                               aaCTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTaaa                                                             <  2:6563513/93‑1 (MQ=255)

                               aaCTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTaaa                                                             <  1:1061634/93‑1 (MQ=255)

                               aaCTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTaaa                                                             >  1:5073927/1‑93 (MQ=255)

                                  tACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCa                                                          >  1:2410657/1‑93 (MQ=255)

                                  tACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCa                                                          >  1:1368496/1‑93 (MQ=255)

                                  tACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCa                                                          >  2:5524899/1‑93 (MQ=255)

                                  tACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCa                                                          >  2:1908276/1‑93 (MQ=255)

                                   aCGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGAAAGAAGATCACATCACCGTTAAAGCAc                                                         <  1:5651676/93‑1 (MQ=255)

                                    cGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTTCTTCCCGACCCGGAAAGTGGAAGAAGATGACTTCACCGTTAAAGCACt                                                        >  1:3222110/1‑93 (MQ=255)

                                    cGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACt                                                        >  1:1674308/1‑93 (MQ=255)

                                     gACCGTCCGTCCTGGACTGGCCTGGTTTTCCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTg                                                       <  2:1190628/93‑1 (MQ=255)

                                     gACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTg                                                       >  2:2866604/1‑93 (MQ=255)

                                      aCCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTgg                                                      <  1:756789/93‑1 (MQ=255)

                                      aCCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTgg                                                      >  2:4260128/1‑93 (MQ=255)

                                      aCCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTgg                                                      <  1:1994413/93‑1 (MQ=255)

                                       ccgtccgtCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGt                                                     <  1:2633454/93‑1 (MQ=255)

                                         gtccgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGa                                                   >  2:1012576/1‑93 (MQ=255)

                                          tccgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGaa                                                  <  2:3931440/93‑1 (MQ=255)

                                          tccgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGaa                                                  >  1:3343757/1‑93 (MQ=255)

                                          tccgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGaa                                                  <  1:7191203/93‑1 (MQ=255)

                                          tccgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGaa                                                  >  2:404856/1‑93 (MQ=255)

                                          tccgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGaa                                                  >  2:4472412/1‑93 (MQ=255)

                                            cgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATg                                                <  2:6891127/93‑1 (MQ=255)

                                            cgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATg                                                >  1:6898047/1‑93 (MQ=255)

                                             gtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGc                                               >  2:3526722/1‑93 (MQ=255)

                                             gtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGc                                               >  2:2270630/1‑93 (MQ=255)

                                              tccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGcc                                              <  2:5011493/93‑1 (MQ=255)

                                              tccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGcc                                              >  2:6649592/1‑93 (MQ=255)

                                              tccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGcc                                              >  1:1442001/1‑93 (MQ=255)

                                                 tGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTAc                                           <  2:4933156/93‑1 (MQ=255)

                                                 tGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTAc                                           >  2:2158108/1‑93 (MQ=255)

                                                 tGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTAc                                           <  1:960221/93‑1 (MQ=255)

                                                 tGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTAc                                           <  2:6263698/93‑1 (MQ=255)

                                                 tGGACTGGCCTGGATTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTAc                                           >  1:633837/1‑93 (MQ=255)

                                                  ggACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACg                                          >  1:4165482/1‑93 (MQ=255)

                                                    aCTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGaa                                        <  2:7392505/93‑1 (MQ=255)

                                                    aCTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGaa                                        <  2:3411347/93‑1 (MQ=255)

                                                      tGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAgg                                      <  2:3845736/93‑1 (MQ=255)

                                                      tGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAgg                                      >  2:4147181/1‑93 (MQ=255)

                                                      tGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAgg                                      <  1:2144314/93‑1 (MQ=255)

                                                       ggCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAggg                                     >  1:6021076/1‑93 (MQ=255)

                                                       ggCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAggg                                     >  1:3660895/1‑93 (MQ=255)

                                                        gCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGc                                    <  2:5240365/93‑1 (MQ=255)

                                                        gCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGc                                    <  1:7402168/93‑1 (MQ=255)

                                                        gCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGc                                    >  2:5739313/1‑93 (MQ=255)

                                                         ccTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCt                                   <  2:468476/93‑1 (MQ=255)

                                                         ccTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCt                                   <  1:1766219/93‑1 (MQ=255)

                                                          cTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTg                                  <  2:2001932/93‑1 (MQ=255)

                                                          cTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTg                                  <  2:5561027/93‑1 (MQ=255)

                                                           tGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGt                                 >  2:5671855/1‑93 (MQ=255)

                                                            ggTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGtt                                <  1:2388947/93‑1 (MQ=255)

                                                            ggTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGtt                                <  1:6500878/93‑1 (MQ=255)

                                                              tttACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTg                              >  2:3180349/1‑93 (MQ=255)

                                                               ttACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTgg                             <  1:5468714/93‑1 (MQ=255)

                                                                     aaCCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAgc                       >  1:2277898/1‑93 (MQ=255)

                                                                        cGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGcc                    >  2:3297785/1‑93 (MQ=255)

                                                                        cGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGcc                    >  2:3663878/1‑93 (MQ=255)

                                                                        cGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGcc                    <  1:5390252/93‑1 (MQ=255)

                                                                         gAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCg                   <  2:4833605/93‑1 (MQ=255)

                                                                          aaTGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCgg                  <  1:2976680/93‑1 (MQ=255)

                                                                           aTGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGgt                 >  1:6260449/1‑93 (MQ=255)

                                                                             gCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGgttg               <  2:5391225/93‑1 (MQ=255)

                                                                                  ttCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGAt          >  2:700816/1‑93 (MQ=255)

                                                                                   tCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATa         <  1:3484333/93‑1 (MQ=255)

                                                                                    cccGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATaa        <  1:849258/93‑1 (MQ=255)

                                                                                    cccGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATaa        <  1:6621973/93‑1 (MQ=255)

                                                                                     ccGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAg       <  2:6177205/93‑1 (MQ=255)

                                                                                       gACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTg     >  1:4313637/1‑93 (MQ=255)

                                                                                       gACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTg     >  1:7481891/1‑93 (MQ=255)

                                                                                       gACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTg     >  1:5203538/1‑93 (MQ=255)

ACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTGGAC  >  NC_000913/3007639‑3007821

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