BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	147,141	A→G	G58G (GGA→GGG)	yadD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	147,141	0	A	G	100.0%	786.2 / NA	228	G58G (GGA→GGG)	yadD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (114/114); total (114/114)

Only 100 of 235 total aligned reads displayed.

GATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAA  >  NC_000913/147051‑147230

gaTACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGgca                                                                                         <  1:5055569/93‑1 (MQ=255)

gaTACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGgca                                                                                         >  1:2010622/1‑93 (MQ=255)

 aTACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGgcac                                                                                        <  1:3701091/93‑1 (MQ=255)

  tACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCGTCATTGTAGAGAGGCTGAAAGGGcaca                                                                                       >  1:3234881/1‑93 (MQ=255)

   aCATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAg                                                                                      <  1:2613973/93‑1 (MQ=255)

   aCATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGGTTCATTGAAGAGAGCCTGAAAGGGAACAg                                                                                      >  2:844731/1‑93 (MQ=255)

    cATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGc                                                                                     >  2:5046082/1‑93 (MQ=255)

      tttGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCAc                                                                                   <  1:3969048/93‑1 (MQ=255)

         gCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGa                                                                                >  1:4148632/1‑93 (MQ=255)

          ccAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGAc                                                                               <  2:3957656/93‑1 (MQ=255)

          ccAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGAc                                                                               <  1:299319/93‑1 (MQ=255)

           cAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACg                                                                              >  2:4971106/1‑93 (MQ=255)

            aGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGt                                                                             >  1:3353866/1‑93 (MQ=255)

            aGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGt                                                                             >  2:6052233/1‑93 (MQ=255)

              tGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGc                                                                           >  1:1976126/1‑93 (MQ=255)

               ggAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGct                                                                          >  1:4526966/1‑93 (MQ=255)

                gAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGctc                                                                         <  2:525767/93‑1 (MQ=255)

                gAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGctc                                                                         <  1:4285386/93‑1 (MQ=255)

                gAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGctc                                                                         <  1:1642065/93‑1 (MQ=255)

                 aaTTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGctct                                                                        >  1:2394597/1‑93 (MQ=255)

                 aaTTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGctct                                                                        >  1:5405034/1‑93 (MQ=255)

                  aTTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTa                                                                       <  2:386633/93‑1 (MQ=255)

                   ttACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTAt                                                                      <  1:2342263/93‑1 (MQ=255)

                   ttACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTAt                                                                      >  1:902751/1‑93 (MQ=255)

                   ttACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTAt                                                                      <  1:35856/93‑1 (MQ=255)

                     aCGCGAACTTTGTGACCTCAACACGGTTCATTTAGAGTCGGTGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTc                                                                    <  1:5309461/93‑1 (MQ=255)

                     aCGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTc                                                                    <  2:3726495/93‑1 (MQ=255)

                       gcgAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCg                                                                  >  2:2825260/1‑93 (MQ=255)

                       gcgAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCg                                                                  >  1:3554942/1‑93 (MQ=255)

                        cgAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGt                                                                 <  1:5260072/93‑1 (MQ=255)

                        cgAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGt                                                                 <  1:4424550/93‑1 (MQ=255)

                          aaCTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGc                                                               >  2:1367698/1‑93 (MQ=255)

                          aaCTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGAACAGCACGGACGTGCTCTATTCCGTGc                                                               >  2:4121552/1‑93 (MQ=255)

                            cTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCaa                                                             >  2:4650461/1‑93 (MQ=255)

                            cTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCaa                                                             >  1:4234174/1‑93 (MQ=255)

                             tttGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCaaa                                                            <  2:5481632/93‑1 (MQ=255)

                              ttGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAAt                                                           <  2:5122400/93‑1 (MQ=255)

                                 tgACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCa                                                        <  1:5666194/93‑1 (MQ=255)

                                 tgACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCa                                                        <  1:1560215/93‑1 (MQ=255)

                                   aCCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGCAAGGGCACAGCACGGACCTGCTCTATTCCGTGCAAATGCAgg                                                      >  1:193868/1‑93 (MQ=255)

                                    ccTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAggg                                                     <  2:1836371/93‑1 (MQ=255)

                                    ccTCAACACGCTGCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAggg                                                     <  2:4191507/93‑1 (MQ=255)

                                     cTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGc                                                    <  1:716594/93‑1 (MQ=255)

                                      tCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCGCTATTCCGTGCAAATGCAGGGCa                                                   >  2:4243278/1‑93 (MQ=255)

                                       cAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCaa                                                  >  2:94875/1‑93 (MQ=255)

                                        aaCACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAgggacac                                                 >  2:560252/1‑89 (MQ=255)

                                         acacGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATc                                                >  2:2067179/1‑93 (MQ=255)

                                         acacGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATc                                                >  2:6007096/1‑93 (MQ=255)

                                          cacGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATcc                                               <  2:2336596/93‑1 (MQ=255)

                                            cGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGGAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGAAGGGCAATCCCg                                             >  1:3420914/1‑93 (MQ=255)

                                            cGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCg                                             <  2:2121651/93‑1 (MQ=255)

                                             gCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCgg                                            >  1:6154297/1‑93 (MQ=255)

                                             gCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCgg                                            >  2:1404122/1‑93 (MQ=255)

                                                tCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                                         <  2:4050203/93‑1 (MQ=255)

                                                 cATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTAt                                        >  2:2488829/1‑93 (MQ=255)

                                                  aTTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATc                                       <  2:359217/93‑1 (MQ=255)

                                                   tttAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCt                                      <  1:3030102/93‑1 (MQ=255)

                                                   tttAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCt                                      <  1:2930652/93‑1 (MQ=255)

                                                   tttAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCt                                      <  1:1314799/93‑1 (MQ=255)

                                                    ttAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTg                                     <  1:4641463/93‑1 (MQ=255)

                                                    ttAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTg                                     >  1:6465313/1‑93 (MQ=255)

                                                      agagTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCGCAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCa                                   <  1:4888179/93‑1 (MQ=255)

                                                      agagTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCa                                   >  1:3475161/1‑93 (MQ=255)

                                                       gagTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAt                                  >  2:2906961/1‑93 (MQ=255)

                                                       gagTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAt                                  >  2:6242483/1‑93 (MQ=255)

                                                        agTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtg                                 <  2:1154103/93‑1 (MQ=255)

                                                        agTCGGGGAGTTTCATTGAAGAGAGCC‑GAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgt                                <  2:6015306/93‑1 (MQ=255)

                                                         gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgt                                <  1:175779/93‑1 (MQ=255)

                                                         gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgt                                <  2:4919518/93‑1 (MQ=255)

                                                            ggggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt                             <  2:5711702/93‑1 (MQ=255)

                                                            ggggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt                             >  1:5576897/1‑93 (MQ=255)

                                                             gggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg                            >  1:1496552/1‑93 (MQ=255)

                                                             gggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg                            >  1:824642/1‑93 (MQ=255)

                                                              ggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGa                           >  2:333361/1‑93 (MQ=255)

                                                               gAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGAt                          >  1:2901513/1‑93 (MQ=255)

                                                                aGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGAtt                         <  2:458898/93‑1 (MQ=255)

                                                                 gTTTCATGGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTg                        <  1:2739934/93‑1 (MQ=255)

                                                                  tttCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGa                       <  2:5352033/93‑1 (MQ=255)

                                                                  tttCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGa                       >  1:6270568/1‑93 (MQ=255)

                                                                   ttCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGaa                      <  2:3641629/93‑1 (MQ=255)

                                                                     cATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAaca                    >  2:4346192/1‑93 (MQ=255)

                                                                     cATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAaca                    >  1:4499994/1‑93 (MQ=255)

                                                                      aTTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAacac                   >  1:1919983/1‑93 (MQ=255)

                                                                      aTTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAacac                   <  1:6209362/93‑1 (MQ=255)

                                                                           agagagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaag              >  1:1119691/1‑93 (MQ=255)

                                                                            gagagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagc             >  2:4703244/1‑93 (MQ=255)

                                                                             agagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCCaagca            >  1:279945/1‑93 (MQ=255)

                                                                              gagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaa           >  2:3988423/1‑93 (MQ=255)

                                                                              gagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaa           <  2:759877/93‑1 (MQ=255)

                                                                              gagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaa           >  2:6469736/1‑93 (MQ=255)

                                                                              gagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaa           >  1:6024787/1‑93 (MQ=255)

                                                                                gCCTGAAAGGGCCCAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaagc         <  1:2114286/93‑1 (MQ=255)

                                                                                gCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaagc         <  1:5857346/93‑1 (MQ=255)

                                                                                gCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaagc         <  1:2441179/93‑1 (MQ=255)

                                                                                 ccTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGcc        <  1:4364443/93‑1 (MQ=255)

                                                                                   tGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCgg      <  2:1741724/93‑1 (MQ=255)

                                                                                     aaaGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGAt    >  2:2198889/1‑93 (MQ=255)

                                                                                      aaGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATa   <  1:6101669/93‑1 (MQ=255)

                                                                                      aaGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATa   >  2:299493/1‑93 (MQ=255)

                                                                                       aGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATaa  <  1:1255387/93‑1 (MQ=255)

GATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAA  >  NC_000913/147051‑147230

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