BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,870,631	Δ1 bp	coding (610/1764 nt)	dgoA ←	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,870,630	0	C	.	100.0%	1367.0 / NA	293	coding (611/1764 nt)	dgoA	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (152/141); total (152/141)

Only 100 of 313 total aligned reads displayed.

CCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA  >  NC_000913/3870538‑3870718

cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                          >  1:2918671/1‑93 (MQ=255)

cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                          <  2:6722995/93‑1 (MQ=255)

cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                          <  2:3601157/93‑1 (MQ=255)

 ccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAct                                                                                         >  2:7457876/1‑92 (MQ=255)

 ccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAct                                                                                         <  1:145994/93‑2 (MQ=255)

  cAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActg                                                                                        >  1:6039872/1‑91 (MQ=255)

    aCCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActgcc                                                                                      <  1:1079631/93‑5 (MQ=255)

    aCCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActgcc                                                                                      <  1:4626453/93‑5 (MQ=255)

     ccGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCct                                                                                    >  1:2738749/1‑93 (MQ=255)

        aCCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCtctc                                                                                 <  1:3524780/93‑1 (MQ=255)

        aCCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCtctc                                                                                 >  1:2616096/1‑93 (MQ=255)

        aCCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCtctc                                                                                 <  2:1520115/93‑1 (MQ=255)

        aCCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCtctc                                                                                 >  1:1402085/1‑93 (MQ=255)

          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGa                                                                               <  1:7418963/93‑1 (MQ=255)

          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGa                                                                               <  1:917844/93‑1 (MQ=255)

          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGa                                                                               <  1:3071862/93‑1 (MQ=255)

           aCGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGAt                                                                              >  2:6277287/1‑93 (MQ=255)

           aCGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGAt                                                                              >  1:4392581/1‑93 (MQ=255)

           aCGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGAt                                                                              >  1:4903978/1‑93 (MQ=255)

             gCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATAc                                                                            <  1:4084493/93‑1 (MQ=255)

               cTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGc                                                                          >  2:5156607/1‑93 (MQ=255)

                ttCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGcc                                                                         >  1:5237650/1‑93 (MQ=255)

                 tCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCt                                                                        >  2:2354672/1‑93 (MQ=255)

                   aTCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTa                                                                      >  2:3017326/1‑93 (MQ=255)

                   aTCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTa                                                                      <  1:5840352/93‑1 (MQ=255)

                    tCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTaa                                                                     >  1:1745365/1‑93 (MQ=255)

                      ggTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACa                                                                   >  1:596372/1‑93 (MQ=255)

                      ggTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACa                                                                   >  2:4985476/1‑93 (MQ=255)

                         ttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAAt                                                                >  2:7241730/1‑93 (MQ=255)

                          tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAAtg                                                               >  2:926349/1‑93 (MQ=255)

                          tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAAtg                                                               >  2:1361465/1‑93 (MQ=255)

                           cAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAAtgc                                                              <  1:3156230/93‑1 (MQ=255)

                            aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAAtgct                                                             >  1:6932407/1‑93 (MQ=255)

                            aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAAtgct                                                             <  1:702682/93‑1 (MQ=255)

                              ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAAtgctgc                                                           <  2:2441825/93‑1 (MQ=255)

                              ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAAtgctgc                                                           <  1:5523648/93‑1 (MQ=255)

                               tttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAAtgctgct                                                          <  2:5718598/93‑1 (MQ=255)

                                    ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                                     >  1:3750513/1‑93 (MQ=255)

                                    ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                                     >  2:6625779/1‑93 (MQ=255)

                                     gAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGctgctgcctgc                                                    <  2:4247371/93‑1 (MQ=255)

                                        cATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctgc                                                 <  1:4279/93‑1 (MQ=255)

                                          tCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                                               <  1:3884974/93‑1 (MQ=255)

                                            cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCgg                                             <  2:3904849/93‑1 (MQ=255)

                                             aGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGt                                            <  2:7333506/93‑1 (MQ=255)

                                               cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                                          <  1:5113694/93‑1 (MQ=255)

                                               cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                                          <  1:6632420/93‑1 (MQ=255)

                                                ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcg                                         >  2:7032208/1‑93 (MQ=255)

                                                 gggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcgc                                        >  2:1906481/1‑93 (MQ=255)

                                                 gggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcgc                                        >  2:1103798/1‑93 (MQ=255)

                                                   gAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACATATGCTGCTGCCTGCTGCGCGGTGCGctc                                      <  2:1213428/93‑1 (MQ=255)

                                                   gAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctc                                      <  1:4111455/93‑1 (MQ=255)

                                                   gAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctc                                      >  2:5819741/1‑93 (MQ=255)

                                                     ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTa                                    <  2:6803188/93‑1 (MQ=255)

                                                     ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTa                                    >  2:7187035/1‑93 (MQ=255)

                                                      gTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGGCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                                   <  1:6214816/93‑1 (MQ=255)

                                                      gTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                                   <  2:537644/93‑1 (MQ=255)

                                                       tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                                  >  2:5329251/1‑93 (MQ=255)

                                                       tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                                  <  2:2808390/93‑1 (MQ=255)

                                                       tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                                  >  1:7343566/1‑93 (MQ=255)

                                                       tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                                  >  2:3684198/1‑93 (MQ=255)

                                                       tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                                  >  2:7430196/1‑93 (MQ=255)

                                                        aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                                 >  2:1762016/1‑93 (MQ=255)

                                                        aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                                 <  2:1417960/93‑1 (MQ=255)

                                                        aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                                 >  1:3646335/1‑93 (MQ=255)

                                                        aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                                 >  2:3284266/1‑93 (MQ=255)

                                                         aaCGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGa                                >  2:4583228/1‑93 (MQ=255)

                                                          aCGATACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt                               <  2:5497544/93‑1 (MQ=255)

                                                             aTACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGc                            >  1:468315/1‑93 (MQ=255)

                                                              tACGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc                           <  1:1291494/93‑1 (MQ=255)

                                                               aCGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGccc                          <  2:2582950/93‑1 (MQ=255)

                                                               aCGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGccc                          <  1:1620850/93‑1 (MQ=255)

                                                               aCGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGccc                          >  1:6340270/1‑93 (MQ=255)

                                                                cGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCg                         >  1:6156893/1‑93 (MQ=255)

                                                                cGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCg                         >  1:7191971/1‑93 (MQ=255)

                                                                cGTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCg                         <  2:6546887/93‑1 (MQ=255)

                                                                 gTGGTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg                        >  2:6096091/1‑93 (MQ=255)

                                                                   ggTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcg                      >  1:2165019/1‑93 (MQ=255)

                                                                   ggTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcg                      >  1:994122/1‑93 (MQ=255)

                                                                    gTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc                     >  1:1929555/1‑93 (MQ=255)

                                                                    gTAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc                     <  2:4745452/93‑1 (MQ=255)

                                                                     tAATTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGcgcg                    <  1:2303465/93‑1 (MQ=255)

                                                                       aTTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGAt                  >  2:378494/1‑93 (MQ=255)

                                                                       aTTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGAt                  >  1:2498949/1‑93 (MQ=255)

                                                                       aTTTTGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGAt                  <  1:2059204/93‑1 (MQ=255)

                                                                         tttGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATag                <  1:5607574/93‑1 (MQ=255)

                                                                          ttGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCTGTGCGCTCTACGGATTGCCCGGCGCGATaga               <  1:3932232/93‑1 (MQ=255)

                                                                          ttGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga               >  2:5560637/1‑93 (MQ=255)

                                                                          ttGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga               <  1:4498483/93‑1 (MQ=255)

                                                                          ttGGTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga               >  2:4048592/1‑93 (MQ=255)

                                                                            ggTGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAgaga             <  2:2129467/93‑1 (MQ=255)

                                                                              tGATTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATc           <  1:2868324/93‑1 (MQ=255)

                                                                                aTTTTCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGc         <  1:1154644/93‑1 (MQ=255)

                                                                                 ttttCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCt        <  2:3972498/93‑1 (MQ=255)

                                                                                 ttttCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCt        >  2:3180752/1‑93 (MQ=255)

                                                                                 ttttCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCt        <  1:6080910/93‑1 (MQ=255)

                                                                                  tttCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTg       >  1:5472162/1‑93 (MQ=255)

                                                                                    tCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGcc     <  1:6722659/93‑1 (MQ=255)

                                                                                    tCATTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGcc     >  1:3282446/1‑93 (MQ=255)

                                                                                      aTTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTa   >  1:2188008/1‑93 (MQ=255)

                                                                                      aTTGCA‑CTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTa   <  1:1013260/93‑1 (MQ=255)

CCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA  >  NC_000913/3870538‑3870718

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