BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,016,772	T→C	C898R (TGC→CGC)	b2878 →	putative oxidoreductase Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,016,772	0	T	C	100.0%	764.0 / NA	237	C898R (TGC→CGC)	b2878	putative oxidoreductase Fe‑S subunit
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (110/127); total (110/127)

Only 100 of 263 total aligned reads displayed.

GCGGAAATCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCTGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGAAC  >  NC_000913/3016682‑3016861

gCGGAAATCTATCAACGTAAAGGCTATATCTCTATCACGTTGGTGAACAGTGACGATCGGGACGAGTTTGTCGCCGAGGACGCCGCTcgccgc                                                                                         >  2:1038554/1‑93 (MQ=255)

   gAAATCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTc                                                                                      <  2:4294641/93‑1 (MQ=255)

    aaaTCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCg                                                                                     >  2:5063693/1‑93 (MQ=255)

    aaaTCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCg                                                                                     >  1:2984771/1‑93 (MQ=255)

    aaaTCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCg                                                                                     >  1:5030283/1‑93 (MQ=255)

     aaTCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGa                                                                                    <  2:2117272/93‑1 (MQ=255)

      aTCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGaa                                                                                   >  1:2211551/1‑93 (MQ=255)

      aTCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGaa                                                                                   <  1:5921342/93‑1 (MQ=255)

       tCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAAt                                                                                  >  1:1844812/1‑93 (MQ=255)

        cTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATg                                                                                 <  1:288609/93‑1 (MQ=255)

         tATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGt                                                                                <  1:5361484/93‑1 (MQ=255)

          aTCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTa                                                                               <  1:3427782/93‑1 (MQ=255)

          aTCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTa                                                                               <  1:2406542/93‑1 (MQ=255)

            cAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAAc                                                                             >  2:909548/1‑93 (MQ=255)

            cAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAAc                                                                             <  2:1585260/93‑1 (MQ=255)

             aaCGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACt                                                                            >  1:5666474/1‑93 (MQ=255)

             aaCGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACt                                                                            >  2:4103553/1‑93 (MQ=255)

             aaCGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACt                                                                            <  2:4698607/93‑1 (MQ=255)

              aCGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTa                                                                           <  1:3044582/93‑1 (MQ=255)

               cGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTAc                                                                          <  1:744158/93‑1 (MQ=255)

                gTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACg                                                                         >  2:5683698/1‑93 (MQ=255)

                 tAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGt                                                                        <  2:2974218/93‑1 (MQ=255)

                 tAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGt                                                                        <  1:1397475/93‑1 (MQ=255)

                  aaaGGCGATATCTCTATCACGCTGGTGTACAGTGACGATCTCGACGCGTTTGTCGCCCAGGAAGCCCGTCGCCGCCTCGATTGTAACTATGtt                                                                       >  2:5183572/1‑93 (MQ=255)

                  aaaGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGtt                                                                       >  1:2707703/1‑93 (MQ=255)

                  aaaGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGtt                                                                       <  1:1809416/93‑1 (MQ=255)

                  aaaGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGtt                                                                       >  1:840935/1‑93 (MQ=255)

                  aaaGGCGATAGCTCTAACACGCTGGTGAACAGTGACGATCGTGACGCGTTCGTCGCCCGGGAAGCCGCTCGCCGCCTCGAATGTAACTACGtt                                                                       <  2:1485701/93‑1 (MQ=255)

                   aaGGCGCTATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGttt                                                                      <  1:2977673/93‑1 (MQ=255)

                   aaGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGttt                                                                      <  2:2835196/93‑1 (MQ=255)

                    aGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTg                                                                     >  2:385469/1‑93 (MQ=255)

                      gCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTgca                                                                   <  1:1511373/93‑1 (MQ=255)

                       cGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTgcag                                                                  >  2:4535454/1‑93 (MQ=255)

                        gATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTgcagc                                                                 <  2:4198696/93‑1 (MQ=255)

                         atatCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTgcagca                                                                >  1:4454841/1‑93 (MQ=255)

                         atatCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTgcagca                                                                >  1:2510092/1‑93 (MQ=255)

                          tatCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCaa                                                               >  2:448975/1‑93 (MQ=255)

                          tatCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGACGCTCGCCGCCTCGAATGTAACTACGTTTgcagcag                                                               >  1:583007/1‑92 (MQ=255)

                            tctcTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgt                                                             <  1:3702001/93‑1 (MQ=255)

                             ttttATCACGCTGGTGAACAGTGACGATCGTGACACGTTTGTCGCCCATGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgtg                                                            <  1:4231653/90‑1 (MQ=255)

                             ctctATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgtg                                                            >  1:888694/1‑93 (MQ=255)

                             ctctATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgtg                                                            <  1:2998388/93‑1 (MQ=255)

                             ctctATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgtg                                                            <  1:4047299/93‑1 (MQ=255)

                              tctATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgtgt                                                           >  2:4171995/1‑93 (MQ=255)

                               ctATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgtgtg                                                          >  1:2855308/1‑93 (MQ=255)

                               ctATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgtgtg                                                          <  2:2568317/93‑1 (MQ=255)

                               ctATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgtgtg                                                          <  2:1956511/93‑1 (MQ=255)

                                tATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAgtgtgt                                                         >  2:1742369/1‑93 (MQ=255)

                                 aTCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGtgtgtg                                                        >  1:5325286/1‑93 (MQ=255)

                                    aCGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGAt                                                     <  1:4955840/93‑1 (MQ=255)

                                      gCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGt                                                   >  1:4192527/1‑93 (MQ=255)

                                       cTGGTGAGCAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTc                                                  <  2:3862998/93‑1 (MQ=255)

                                        tGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCt                                                 >  1:1062402/1‑93 (MQ=255)

                                          gTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGc                                               <  2:3038743/93‑1 (MQ=255)

                                           tGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGcc                                              <  1:5502712/93‑1 (MQ=255)

                                           tGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGcc                                              <  1:2481294/93‑1 (MQ=255)

                                           tGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGcc                                              >  1:5635217/1‑93 (MQ=255)

                                            gAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGccc                                             >  2:1361244/1‑93 (MQ=255)

                                            gAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGccc                                             >  1:5763802/1‑93 (MQ=255)

                                              aCAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGa                                           >  2:4847400/1‑93 (MQ=255)

                                                aGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAAc                                         <  1:3636159/93‑1 (MQ=255)

                                                 gTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGACTGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAAcc                                        >  2:3560616/1‑93 (MQ=255)

                                                 gTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAAcc                                        <  2:696315/93‑1 (MQ=255)

                                                 gTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAAcc                                        <  1:642876/93‑1 (MQ=255)

                                                 gTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAAcc                                        >  2:1807981/1‑93 (MQ=255)

                                                       aTCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCa                                  >  1:4147024/1‑93 (MQ=255)

                                                       aTCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCa                                  >  1:3296189/1‑93 (MQ=255)

                                                       aTCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCa                                  >  1:4714373/1‑93 (MQ=255)

                                                        tCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCaa                                 >  2:1248599/1‑93 (MQ=255)

                                                        tCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCaa                                 <  2:5258221/93‑1 (MQ=255)

                                                          gTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACg                               <  1:1672373/93‑1 (MQ=255)

                                                          gTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACg                               >  2:5979568/1‑93 (MQ=255)

                                                           tGACGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGt                              <  1:2638181/93‑1 (MQ=255)

                                                             aCGCGTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGtct                            >  2:1159002/1‑93 (MQ=255)

                                                               ccgTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTcc                          <  2:1617615/92‑1 (MQ=255)

                                                                cgTTTGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTcct                         >  1:4849188/1‑92 (MQ=255)

                                                                 ctctGTCGCGCATGAACCCTCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGAGTGTGAATGTCTGCCCGATCCGCGCCAACGTCTCCAt                        <  2:3084632/90‑1 (MQ=255)

                                                                  tttGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCAtt                       <  2:3491726/93‑1 (MQ=255)

                                                                  tttGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCAtt                       <  2:2256664/93‑1 (MQ=255)

                                                                   ttGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTg                      <  1:6248850/93‑1 (MQ=255)

                                                                    tgcggCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGc                     <  2:5380142/89‑1 (MQ=255)

                                                                    tGTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGc                     >  2:138146/1‑93 (MQ=255)

                                                                     gTCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCg                    <  1:5285803/93‑1 (MQ=255)

                                                                      tCGCCCAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCgg                   <  2:4395746/93‑1 (MQ=255)

                                                                         cccAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTcc                >  1:3477135/1‑93 (MQ=255)

                                                                         cccAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTcc                >  2:5882825/1‑93 (MQ=255)

                                                                          ccAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTccc               <  1:4977924/93‑1 (MQ=255)

                                                                           cAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCa              <  2:3660634/93‑1 (MQ=255)

                                                                           cAGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCa              <  1:4664831/93‑1 (MQ=255)

                                                                            aGGAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAg             >  1:2049869/1‑93 (MQ=255)

                                                                              gAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGc           <  2:644671/93‑1 (MQ=255)

                                                                              gAAGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGGGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGc           >  1:1243018/1‑93 (MQ=255)

                                                                                aGCCGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCtt         >  1:5150686/1‑93 (MQ=255)

                                                                                  ccGCTCGCCGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTGTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTcc       <  2:6108751/93‑1 (MQ=255)

                                                                                   cgctcgccGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCa      <  2:5572834/93‑1 (MQ=255)

                                                                                   cgctcgccGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCa      >  1:4285857/1‑93 (MQ=255)

                                                                                   cgctcgccGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCa      <  1:3818964/93‑1 (MQ=255)

                                                                                   cgctcgccGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCa      >  2:806762/1‑93 (MQ=255)

                                                                                     ctcgccGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGa    >  2:275798/1‑93 (MQ=255)

                                                                                      tcgccGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGaa   >  2:3988120/1‑93 (MQ=255)

GCGGAAATCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCTGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGAAC  >  NC_000913/3016682‑3016861

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