BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	281,747	C→T	T89T (ACC→ACT)	yagE →	putative lyase/synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	281,747	0	C	T	84.7%	681.3 / 74.7	249	T89T (ACC→ACT)	yagE	putative lyase/synthase
Reads supporting (aligned to +/- strand): ref base C (25/13); new base T (117/93); total (143/106)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 260 total aligned reads displayed.

GAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTAC  >  NC_000913/281655‑281837

gAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCAcc                                                                                            <  1:6856239/93‑1 (MQ=255)

  gTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACtgg                                                                                          >  2:7058422/1‑93 (MQ=255)

  gTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACtgg                                                                                          <  2:6679606/93‑1 (MQ=255)

   ttCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACtggc                                                                                         <  1:649906/93‑1 (MQ=255)

    tctcCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACtggcg                                                                                        >  1:4887656/1‑93 (MQ=255)

    tctcCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACtggcg                                                                                        <  1:2302319/93‑1 (MQ=255)

     ctcCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACcggcgg                                                                                       >  1:4770638/1‑93 (MQ=255)

      tcCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCggcggc                                                                                      >  2:5535912/1‑93 (MQ=255)

       cccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCa                                                                                     <  1:2852092/93‑1 (MQ=255)

        ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCAc                                                                                    <  2:1436496/93‑1 (MQ=255)

         cAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCAcc                                                                                   <  2:1381764/93‑1 (MQ=255)

          aGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCa                                                                                  >  2:7384159/1‑93 (MQ=255)

           gCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCaa                                                                                 >  1:3236638/1‑93 (MQ=255)

            cTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAAc                                                                                >  2:5421131/1‑93 (MQ=255)

               ggCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGcc                                                                             <  1:3495471/93‑1 (MQ=255)

                  gcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCgg                                                                          <  1:262820/93‑1 (MQ=255)

                    cGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGa                                                                        >  1:1618747/1‑93 (MQ=255)

                     gAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGaa                                                                       >  1:7422049/1‑93 (MQ=255)

                     gAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGaa                                                                       >  2:1579174/1‑93 (MQ=255)

                     gAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGaa                                                                       >  1:5090390/1‑93 (MQ=255)

                        gagCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAAcc                                                                    >  2:5574214/1‑93 (MQ=255)

                          gCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCAt                                                                  >  1:1734461/1‑93 (MQ=255)

                            gTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCg                                                                >  1:6294353/1‑93 (MQ=255)

                            gTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCg                                                                >  2:216038/1‑93 (MQ=255)

                              aaaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGaa                                                              <  1:5219934/93‑1 (MQ=255)

                                   cATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTcag                                                         <  2:3043707/93‑1 (MQ=255)

                                   cATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTcag                                                         >  1:1275379/1‑93 (MQ=255)

                                    aTTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTcagc                                                        >  2:828770/1‑93 (MQ=255)

                                     ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagcc                                                       <  2:1715238/93‑1 (MQ=255)

                                      tGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTcagcca                                                      >  2:5100511/1‑93 (MQ=255)

                                         ccGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCa                                                   >  1:727986/1‑93 (MQ=255)

                                         ccGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCa                                                   <  1:1166794/93‑1 (MQ=255)

                                         ccGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCa                                                   >  2:3918314/1‑93 (MQ=255)

                                         ccGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCa                                                   >  1:1688177/1‑93 (MQ=255)

                                          cGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAc                                                  >  2:4742949/1‑93 (MQ=255)

                                          cGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAc                                                  <  2:2747256/93‑1 (MQ=255)

                                           gCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAcg                                                 >  1:1957150/1‑93 (MQ=255)

                                           gCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAcg                                                 >  2:4451902/1‑93 (MQ=255)

                                           gCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGCGCTGATCGGCACAGGCGGCACCAACGCCCGGGAAAACATCGAACTCAGCCAGCAcg                                                 >  1:4691372/1‑93 (MQ=255)

                                             tttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAcgcg                                               >  1:141536/1‑93 (MQ=255)

                                              ttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACgcgc                                              >  1:2211723/1‑93 (MQ=255)

                                              ttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACgcgc                                              >  1:6371046/1‑93 (MQ=255)

                                              ttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACgcgc                                              >  2:6959723/1‑93 (MQ=255)

                                              ttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACgcgc                                              <  1:3895640/93‑1 (MQ=255)

                                                  tATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCgcagca                                          >  2:2440275/1‑93 (MQ=255)

                                                  tATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCgcagca                                          <  2:6549065/93‑1 (MQ=255)

                                                    tCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAgg                                        >  1:3092658/1‑93 (MQ=255)

                                                    tCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAgg                                        >  2:3320189/1‑93 (MQ=255)

                                                    tCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAgg                                        >  1:5277899/1‑93 (MQ=255)

                                                     cGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggc                                       <  2:5846622/93‑1 (MQ=255)

                                                       aTCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcgg                                     <  2:4869839/93‑1 (MQ=255)

                                                       aTCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcgg                                     >  2:4059690/1‑93 (MQ=255)

                                                          aTGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcgggcg                                  >  1:973911/1‑93 (MQ=255)

                                                          aTGTCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcgggcg                                  >  1:1774492/1‑93 (MQ=255)

                                                             tCGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCgg                               >  2:5834540/1‑93 (MQ=255)

                                                              cGATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGa                              <  1:541586/93‑1 (MQ=255)

                                                               gATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGAc                             <  2:6814133/93‑1 (MQ=255)

                                                               gATCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGAc                             >  1:4582716/1‑93 (MQ=255)

                                                               gATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAAAGCCCCGGAAACCATCGAACTCAGCCCGCCCGCGCAGCAGGCGGGGGCgggg                             >  1:2773480/1‑91 (MQ=255)

                                                                aTCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACg                            >  2:630328/1‑93 (MQ=255)

                                                                aTCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACg                            >  1:5754043/1‑93 (MQ=255)

                                                                aTCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACg                            <  2:6201513/93‑1 (MQ=255)

                                                                aTCGTCGCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACg                            >  1:3274140/1‑93 (MQ=255)

                                                                aTCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACg                            >  1:3597124/1‑93 (MQ=255)

                                                                aTCGTCGCGTGCCCGGGCTGGACGGCACCGGCGGCACCAACACCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGcgggcgcgggcg                            >  1:241189/1‑93 (MQ=255)

                                                                 tcgtcgCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACgg                           <  1:6188885/93‑1 (MQ=255)

                                                                  cgtcgCGTGCCGGTGCTGATCGGCACTTGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGc                          >  2:6418438/1‑93 (MQ=255)

                                                                   gtcgCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCa                         >  1:2724754/1‑93 (MQ=255)

                                                                    tcgCGTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCAt                        <  1:6937253/93‑1 (MQ=255)

                                                                       cgTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgt                     <  2:4215592/93‑1 (MQ=255)

                                                                       cgTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgt                     <  2:1495016/93‑1 (MQ=255)

                                                                       cgTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgt                     <  1:164519/93‑1 (MQ=255)

                                                                       cgTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgt                     <  1:5975864/93‑1 (MQ=255)

                                                                        gTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtg                    <  1:5507570/93‑1 (MQ=255)

                                                                        gTGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtg                    <  1:998053/93‑1 (MQ=255)

                                                                        gTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtg                    >  2:4505634/1‑93 (MQ=255)

                                                                         tGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                   >  2:6488050/1‑93 (MQ=255)

                                                                         tGCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                   <  1:378836/93‑1 (MQ=255)

                                                                          gCCGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtggt                  <  1:1040404/93‑1 (MQ=255)

                                                                           ccGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtggtg                 <  2:3718906/93‑1 (MQ=255)

                                                                            cGGTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGa                <  1:2426044/93‑1 (MQ=255)

                                                                             ggTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGAt               <  1:2820102/93‑1 (MQ=255)

                                                                              gTGCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATc              >  2:3954899/1‑93 (MQ=255)

                                                                                gCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCaa            >  2:4954592/1‑93 (MQ=255)

                                                                                gCTGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCaa            <  2:193979/93‑1 (MQ=255)

                                                                                gCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCGAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCaa            >  1:7298999/1‑93 (MQ=255)

                                                                                  tGATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAcc          >  2:7500442/1‑93 (MQ=255)

                                                                                  tGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAcc          >  2:2306566/1‑93 (MQ=255)

                                                                                   gATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAccc         >  2:2912198/1‑93 (MQ=255)

                                                                                   gATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAccc         >  2:597331/1‑93 (MQ=255)

                                                                                   gATCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAccc         >  2:4325672/1‑93 (MQ=255)

                                                                                    aTCGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAcccc        <  2:6878826/93‑1 (MQ=255)

                                                                                    aTCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAcccc        >  1:5648419/1‑93 (MQ=255)

                                                                                      cGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCcta      <  2:6015388/93‑1 (MQ=255)

                                                                                      cGGCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCcta      <  1:6666686/93‑1 (MQ=255)

                                                                                       ggCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctac     <  1:3908452/93‑1 (MQ=255)

                                                                                       ggCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctac     <  1:4464072/93‑1 (MQ=255)

                                                                                        gCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctact    >  2:519168/1‑93 (MQ=255)

                                                                                        gCACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctact    >  2:6261579/1‑93 (MQ=255)

                                                                                         cACTGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctacta   >  2:1519734/1‑93 (MQ=255)

GAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTAC  >  NC_000913/281655‑281837

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