BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,549,377	T→C	V58A (GTG→GCG)	uxuA →	mannonate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,549,377	0	T	C	89.2%	861.0 / 74.4	325	V58A (GTG→GCG)	uxuA	mannonate hydrolase
Reads supporting (aligned to +/- strand): ref base T (23/12); new base C (127/163); total (150/175)
Fisher's exact test for biased strand distribution p-value = 1.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 333 total aligned reads displayed.

TTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATCAGCAGACTC  >  NC_000913/4549287‑4549469

ttACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGgcgt                                                                                            <  2:2836417/93‑1 (MQ=255)

  aCCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTgg                                                                                          >  1:5513918/1‑93 (MQ=255)

  aCCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTgg                                                                                          <  2:778419/93‑1 (MQ=255)

   ccGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGt                                                                                         >  2:1449808/1‑93 (MQ=255)

    cgcgCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTc                                                                                        <  2:1788409/93‑1 (MQ=255)

    cgcgCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTc                                                                                        <  1:4577353/93‑1 (MQ=255)

     gcgcTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCt                                                                                       >  2:5876267/1‑93 (MQ=255)

     gcgcTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCt                                                                                       <  1:3632801/93‑1 (MQ=255)

       gcTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGt                                                                                     >  2:1643953/1‑93 (MQ=255)

       gcTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGt                                                                                     <  2:3803658/93‑1 (MQ=255)

       gcTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGt                                                                                     <  1:2094839/93‑1 (MQ=255)

        cTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCTGTc                                                                                    <  1:3753925/93‑1 (MQ=255)

        cTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTc                                                                                    <  1:6633768/93‑1 (MQ=255)

         tGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCg                                                                                   >  2:3696125/1‑93 (MQ=255)

          gCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGCCTGTCGt                                                                                  >  1:1668890/1‑93 (MQ=255)

           cACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCTGTCGTa                                                                                 <  2:5760206/93‑1 (MQ=255)

            aCCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAg                                                                                <  1:81708/93‑1 (MQ=255)

                tatCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAg                                                                            >  2:4738136/1‑93 (MQ=255)

                tatCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAg                                                                            <  2:1416981/93‑1 (MQ=255)

                 atCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGc                                                                           <  1:6282898/93‑1 (MQ=255)

                 atCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGc                                                                           >  1:1394052/1‑93 (MQ=255)

                  tCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCg                                                                          <  1:4756930/93‑1 (MQ=255)

                  tCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCg                                                                          >  1:5951980/1‑93 (MQ=255)

                   cccGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGt                                                                         <  2:3182603/93‑1 (MQ=255)

                      gAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGcc                                                                      >  1:7492747/1‑93 (MQ=255)

                       aaCGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCa                                                                     >  1:7044597/1‑93 (MQ=255)

                        aCGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCaa                                                                    <  1:1017153/93‑1 (MQ=255)

                         cGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAAt                                                                   <  1:2012606/93‑1 (MQ=255)

                          ggCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCTGTCGTAGAAAGCGTGCCAAtt                                                                  >  1:143590/1‑93 (MQ=255)

                          ggCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAAtt                                                                  <  1:3495775/93‑1 (MQ=255)

                          ggCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAAtt                                                                  <  2:6120776/93‑1 (MQ=255)

                            cGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCa                                                                <  2:1910182/93‑1 (MQ=255)

                             gAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCAc                                                               >  2:2426658/1‑93 (MQ=255)

                             gAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCAc                                                               >  2:5298372/1‑93 (MQ=255)

                              aaGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACg                                                              >  1:254548/1‑93 (MQ=255)

                               aGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGa                                                             >  1:6121762/1‑93 (MQ=255)

                               aGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGa                                                             >  1:2809073/1‑93 (MQ=255)

                                gTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGaa                                                            >  1:4935177/1‑93 (MQ=255)

                                 tATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAg                                                           <  1:7348317/93‑1 (MQ=255)

                                 tATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAg                                                           <  2:5369331/93‑1 (MQ=255)

                                  aTGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGa                                                          <  1:1214592/93‑1 (MQ=255)

                                  aTGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGa                                                          <  1:1582453/93‑1 (MQ=255)

                                  aTGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGa                                                          <  1:4811423/93‑1 (MQ=255)

                                    ggTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGata                                                        <  1:5048231/93‑1 (MQ=255)

                                      tCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATc                                                      >  2:415723/1‑93 (MQ=255)

                                        cGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCaa                                                    >  2:553031/1‑93 (MQ=255)

                                           tGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAAc                                                 <  2:6829812/93‑1 (MQ=255)

                                             aaGAGATCCTCAAACGCAAGGCGATCATTTAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAAccc                                               >  1:6206706/1‑93 (MQ=255)

                                              agagATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCca                                              >  2:2309235/1‑93 (MQ=255)

                                               gagaTCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCcac                                             >  1:3880020/1‑93 (MQ=255)

                                               gagaTCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCcac                                             <  2:834706/93‑1 (MQ=255)

                                                agaTCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCcaca                                            <  2:5223937/93‑1 (MQ=255)

                                                 gaTCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCacac                                           >  1:6366350/1‑93 (MQ=255)

                                                  aTCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACt                                          <  1:3296877/93‑1 (MQ=255)

                                                   tCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTg                                         <  2:6506875/93‑1 (MQ=255)

                                                   tCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTg                                         <  2:2383505/93‑1 (MQ=255)

                                                   tCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTg                                         >  2:1041287/1‑93 (MQ=255)

                                                   tCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTg                                         >  1:7414076/1‑93 (MQ=255)

                                                    ccTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTgg                                        >  1:3936721/1‑93 (MQ=255)

                                                    ccTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTgg                                        <  2:7266542/93‑1 (MQ=255)

                                                     cTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTTGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGc                                       >  1:341684/1‑93 (MQ=255)

                                                      tCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCa                                      >  2:1596069/1‑93 (MQ=255)

                                                       cAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCaa                                     <  1:5909421/93‑1 (MQ=255)

                                                        aaaCGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAAc                                    >  2:4496168/1‑93 (MQ=255)

                                                            gCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATg                                >  2:5616668/1‑93 (MQ=255)

                                                            gCAAGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATg                                >  1:2357171/1‑93 (MQ=255)

                                                               aGGCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGc                             >  1:943239/1‑93 (MQ=255)

                                                                ggCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCa                            <  2:2160540/93‑1 (MQ=255)

                                                                ggCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCa                            <  1:4253115/93‑1 (MQ=255)

                                                                ggCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCa                            <  1:4459475/93‑1 (MQ=255)

                                                                 gCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAg                           <  2:3343361/93‑1 (MQ=255)

                                                                 gCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAg                           <  2:7378582/93‑1 (MQ=255)

                                                                 gCGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCAACACTGGCAACTATGAGCAg                           <  1:5623353/93‑1 (MQ=255)

                                                                  cGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGt                          <  2:374675/93‑1 (MQ=255)

                                                                  cGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGt                          >  2:136210/1‑93 (MQ=255)

                                                                  cGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGt                          <  1:6827956/93‑1 (MQ=255)

                                                                  cGATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGt                          >  2:7128594/1‑93 (MQ=255)

                                                                   gATCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTg                         >  1:1841118/1‑93 (MQ=255)

                                                                    aTCATTGAAGACGCAGGCCTGGTGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTgg                        <  2:4295641/93‑1 (MQ=255)

                                                                    aTCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTgg                        >  1:6928016/1‑93 (MQ=255)

                                                                    aTCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTgg                        >  2:2088138/1‑93 (MQ=255)

                                                                     tCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGa                       <  2:6391187/93‑1 (MQ=255)

                                                                     tCATTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGa                       >  2:5005313/1‑93 (MQ=255)

                                                                       aTTGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATc                     >  2:6712187/1‑93 (MQ=255)

                                                                         tGAAGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGc                   <  2:872152/93‑1 (MQ=255)

                                                                          gAAGACGCAGGCCTGGCGTAGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCt                  <  2:3582558/93‑1 (MQ=255)

                                                                            aGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTaa                <  2:1222466/93‑1 (MQ=255)

                                                                            aGACGCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTaa                <  2:1106994/93‑1 (MQ=255)

                                                                                gCAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTAt            >  2:6587317/1‑93 (MQ=255)

                                                                                 cAGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATc           <  2:4012805/93‑1 (MQ=255)

                                                                                  aGGCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATca          >  2:3039761/1‑93 (MQ=255)

                                                                                   ggCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATcag         >  2:4680588/1‑93 (MQ=255)

                                                                                   ggCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATcag         >  1:5268691/1‑93 (MQ=255)

                                                                                   ggCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATcag         <  1:4047250/93‑1 (MQ=255)

                                                                                   ggCCTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATcag         >  1:317872/1‑93 (MQ=255)

                                                                                     ccTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATcagca       <  1:5679544/93‑1 (MQ=255)

                                                                                      cTGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATcagcag      >  2:2552878/1‑93 (MQ=255)

                                                                                       tGGTGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATCAGCAGa     >  1:435546/1‑93 (MQ=255)

                                                                                       tGGCGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATCAGCAGa     <  1:1287710/93‑1 (MQ=255)

                                                                                          tgtgGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATCAGCAGActc  >  2:7163115/1‑93 (MQ=255)

TTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCATTGAAGACGCAGGCCTGGTGTGGTCTGTCGTAGAAAGCGTGCCAATTCACGAAGATATCAAAACCCACACTGGCAACTATGAGCAGTGGATCGCTAACTATCAGCAGACTC  >  NC_000913/4549287‑4549469

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