BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,210,945	C→T	S86S (AGC→AGT)	rpoD →	RNA polymerase sigma(70)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,210,945	0	C	T	100.0%	627.9 / NA	178	S86S (AGC→AGT)	rpoD	RNA polymerase sigma(70)
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (87/91); total (87/91)

Only 100 of 191 total aligned reads displayed.

GAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGA  >  NC_000913/3210856‑3211034

gaagaaGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGtgtg                                                                                        >  1:1730059/1‑93 (MQ=255)

gaagaaGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACAACGCGGACGAAGATGATGCTGAAGGCGCCGCGCAGCTGCTTTACAGtgtg                                                                                        >  2:552796/1‑93 (MQ=255)

  agaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGa                                                                                      <  1:1064589/93‑1 (MQ=255)

  agaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGa                                                                                      <  2:3934151/93‑1 (MQ=255)

   gaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGGAGCCGCCGCGCAGGTGCTTTCCAGTGTggga                                                                                     >  2:1741654/1‑91 (MQ=255)

   gaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGaa                                                                                     >  1:3601766/1‑93 (MQ=255)

   gaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGaa                                                                                     <  2:224497/93‑1 (MQ=255)

   gaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTACAGTGTggga                                                                                     >  1:3456039/1‑91 (MQ=255)

   gaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGGGCTTTCCAGTGTGGaa                                                                                     >  2:4029652/1‑93 (MQ=255)

     agCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATc                                                                                   >  1:3493006/1‑93 (MQ=255)

     agCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATc                                                                                   >  2:1849043/1‑93 (MQ=255)

     agCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATc                                                                                   <  2:1437018/93‑1 (MQ=255)

     agCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATc                                                                                   >  2:137176/1‑93 (MQ=255)

     agCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGCGCTTTCCAGTGTGGGATc                                                                                   >  2:445720/1‑93 (MQ=255)

      gCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCt                                                                                  <  2:2478129/93‑1 (MQ=255)

      gCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCt                                                                                  >  2:1079309/1‑93 (MQ=255)

        aCCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGa                                                                                <  2:3700646/93‑1 (MQ=255)

         ccGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaa                                                                               >  1:440334/1‑93 (MQ=255)

          cGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaaa                                                                              <  1:2228296/93‑1 (MQ=255)

          cGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaaa                                                                              <  2:2270736/93‑1 (MQ=255)

          cGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaaa                                                                              >  2:1962178/1‑93 (MQ=255)

          cGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaaa                                                                              <  2:159040/93‑1 (MQ=255)

            gATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATc                                                                            <  2:656117/93‑1 (MQ=255)

               gCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCggg                                                                         <  2:1711138/93‑1 (MQ=255)

               gCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCggg                                                                         <  1:3170534/93‑1 (MQ=255)

                ccGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGgc                                                                        <  1:401946/93‑1 (MQ=255)

                  gatgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGgcgc                                                                      >  2:1103967/1‑93 (MQ=255)

                  gatgatCTGATGCTGGCTGAAAACAACGCGGACGACGATGCTGTCGAAGCCGACGCGCAGGAGCTTTCTAGTGTGCAATCTGTAATCGGgcgc                                                                      >  2:2507916/1‑93 (MQ=255)

                   atgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCa                                                                     >  1:1775859/1‑93 (MQ=255)

                    tgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCAc                                                                    <  1:2545236/93‑1 (MQ=255)

                     gatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACg                                                                   >  1:4207859/1‑93 (MQ=255)

                      atCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGa                                                                  <  1:1970570/93‑1 (MQ=255)

                       tCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGAc                                                                 <  1:2382729/93‑1 (MQ=255)

                         tGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTg                                                               >  1:261516/1‑93 (MQ=255)

                           aTGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCCCGACTGAc                                                             >  1:4673473/1‑93 (MQ=255)

                           aTGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGAc                                                             >  2:3393524/1‑93 (MQ=255)

                            tGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGAcc                                                            >  1:1539030/1‑93 (MQ=255)

                             gctggctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGAccc                                                           >  2:3259723/1‑93 (MQ=255)

                               tggctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCgg                                                         <  1:2316030/93‑1 (MQ=255)

                                   tgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACg                                                     >  1:3694452/1‑93 (MQ=255)

                                    gAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGc                                                    <  1:494535/93‑1 (MQ=255)

                                    caaaaCACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGc                                                    <  1:2030576/92‑1 (MQ=255)

                                     aaaaCACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCa                                                   <  2:2323845/93‑1 (MQ=255)

                                     aaaaCACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCa                                                   <  2:3927878/93‑1 (MQ=255)

                                      aaaCACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCAt                                                  >  2:1787070/1‑93 (MQ=255)

                                       aaCACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATg                                                 >  1:2347089/1‑93 (MQ=255)

                                        acacCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGt                                                >  2:3025254/1‑93 (MQ=255)

                                           ccGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACa                                             >  2:2647510/1‑93 (MQ=255)

                                           ccGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACa                                             <  2:4583937/93‑1 (MQ=255)

                                           ccGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACa                                             >  1:906713/1‑93 (MQ=255)

                                           ccGCGGACGAAGACGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACa                                             <  2:4276716/93‑1 (MQ=255)

                                            cgcgGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACAt                                            >  1:629550/1‑93 (MQ=255)

                                              cgGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGc                                          <  1:1373385/93‑1 (MQ=255)

                                               ggACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCg                                         >  2:2899248/1‑93 (MQ=255)

                                               ggACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCg                                         <  2:1513058/93‑1 (MQ=255)

                                               ggACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCg                                         >  1:2471522/1‑93 (MQ=255)

                                                 aCGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTg                                       <  2:4550163/93‑1 (MQ=255)

                                                  cGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGa                                      >  1:1121533/1‑93 (MQ=255)

                                                   gAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaa                                     >  1:64789/1‑93 (MQ=255)

                                                   gAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaa                                     >  1:542812/1‑93 (MQ=255)

                                                    aaGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                    <  2:2088924/93‑1 (MQ=255)

                                                    aaGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                    <  2:795783/93‑1 (MQ=255)

                                                    aaGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                    <  2:3178549/93‑1 (MQ=255)

                                                    aaGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                    >  1:3135689/1‑93 (MQ=255)

                                                    aaGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGAGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                    >  2:763808/1‑93 (MQ=255)

                                                      gATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATg                                  >  2:1489525/1‑93 (MQ=255)

                                                       aTGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATgg                                 >  1:865088/1‑93 (MQ=255)

                                                       aTGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATgg                                 >  1:2527455/1‑93 (MQ=255)

                                                        tgctgcCGAAGCCGCCGCGCAGGTGTTTTCCCGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATggg                                >  2:3694233/1‑93 (MQ=255)

                                                         gctgcCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGc                               >  1:702294/1‑93 (MQ=255)

                                                         gctgcCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGc                               <  1:3550466/93‑1 (MQ=255)

                                                          ctgcCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCa                              <  1:1475792/93‑1 (MQ=255)

                                                           tgcCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAc                             >  2:2443087/1‑93 (MQ=255)

                                                             ccGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCg                           <  1:2437389/93‑1 (MQ=255)

                                                                    cgccgcGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                    <  2:3082666/93‑1 (MQ=255)

                                                                     gccgcGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTg                   >  1:2978070/1‑93 (MQ=255)

                                                                      ccgcGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGt                  <  1:1156442/93‑1 (MQ=255)

                                                                        gcgcAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTg                >  1:171392/1‑93 (MQ=255)

                                                                        gcgcAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTg                <  2:2572335/93‑1 (MQ=255)

                                                                         cgcAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGa               <  1:2323617/93‑1 (MQ=255)

                                                                          gcAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGAc              <  1:1784906/93‑1 (MQ=255)

                                                                            aGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGAccc            <  1:1980177/93‑1 (MQ=255)

                                                                             ggTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcg           >  2:3561691/1‑93 (MQ=255)

                                                                              gTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgc          <  2:4064734/93‑1 (MQ=255)

                                                                              gTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgc          >  1:2204808/1‑93 (MQ=255)

                                                                               tGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgcg         >  1:4316950/1‑93 (MQ=255)

                                                                               tGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgcg         >  1:3736766/1‑93 (MQ=255)

                                                                               tGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgcg         >  1:3510633/1‑93 (MQ=255)

                                                                               tGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgcg         <  1:805213/93‑1 (MQ=255)

                                                                                gCTTTCCAGTGTGGAATCTGAAATCGGGTGCACGACTGACCCGGTACGCATGTACATGCGGGAAATGGGCACCGTTGAAGTGTTGACCCGCGa        >  2:3840140/1‑93 (MQ=255)

                                                                                gCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGa        <  1:1246309/93‑1 (MQ=255)

                                                                                 cTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCTTTGAACTGTTGACCCGCGaa       <  2:2666963/93‑1 (MQ=255)

                                                                                 cTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGaa       <  1:3283294/93‑1 (MQ=255)

                                                                                 cTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGaa       <  1:1637733/93‑1 (MQ=255)

                                                                                 cTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGaa       >  1:1748954/1‑93 (MQ=255)

                                                                                  tttCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg      <  2:4398318/93‑1 (MQ=255)

                                                                                  tttCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg      >  2:135664/1‑93 (MQ=255)

                                                                                   ttCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAgg     <  2:3783403/93‑1 (MQ=255)

                                                                                    tCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGc    <  2:301309/93‑1 (MQ=255)

                                                                                      cAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGa  >  1:2727683/1‑93 (MQ=255)

GAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGA  >  NC_000913/3210856‑3211034

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