BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,067,548	T→C	H137R (CAC→CGC)	b3879 ←	putative aldose‑1‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,067,548	0	T	C	100.0%	1346.7 / NA	419	H137R (CAC→CGC)	b3879	putative aldose‑1‑epimerase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (179/240); total (179/240)

Only 100 of 467 total aligned reads displayed.

TTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTG  >  NC_000913/4067458‑4067637

ttGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgt                                                                                         <  2:5210474/93‑1 (MQ=255)

 tGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCcgtgcgtg                                                                                        <  1:2183893/93‑1 (MQ=255)

    aGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTgcgtgcgt                                                                                     >  2:2659154/1‑93 (MQ=255)

     ggTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGtt                                                                                    >  2:7191828/1‑93 (MQ=255)

      gTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTg                                                                                   >  2:5769713/1‑93 (MQ=255)

      gTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTg                                                                                   <  2:1485865/93‑1 (MQ=255)

      gTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTg                                                                                   >  2:2586351/1‑93 (MQ=255)

        tAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAg                                                                                 >  1:8123176/1‑93 (MQ=255)

        tAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAg                                                                                 >  1:2075330/1‑93 (MQ=255)

           aTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGa                                                                              <  1:6932757/93‑1 (MQ=255)

           aTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGa                                                                              >  1:4759813/1‑93 (MQ=255)

            tAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGaa                                                                             >  2:7000751/1‑93 (MQ=255)

              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           <  1:2560744/93‑1 (MQ=255)

              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           <  2:2864149/93‑1 (MQ=255)

              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           >  1:5052590/1‑93 (MQ=255)

               ggATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaa                                                                          <  1:1471180/93‑1 (MQ=255)

                gATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaaa                                                                         >  2:8910592/1‑93 (MQ=255)

                 aTGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaac                                                                        <  2:2220210/93‑1 (MQ=255)

                 aTGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaac                                                                        <  1:3983834/93‑1 (MQ=255)

                  tGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaaca                                                                       >  1:7982836/1‑93 (MQ=255)

                     aTACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACt                                                                    <  2:4481183/93‑1 (MQ=255)

                       aCCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtg                                                                  <  1:5553556/93‑1 (MQ=255)

                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 <  2:4600469/93‑1 (MQ=255)

                          cACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTTCGTTGAGCGAATAAACAACTTGTGa                                                               >  2:7972879/1‑93 (MQ=255)

                          cACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGa                                                               <  1:10035914/93‑1 (MQ=255)

                          cACGACATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGa                                                               <  1:3275090/93‑1 (MQ=255)

                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              <  1:6689353/93‑1 (MQ=255)

                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              >  2:9579580/1‑93 (MQ=255)

                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              <  2:8472587/93‑1 (MQ=255)

                            cGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAg                                                             >  1:8068744/1‑93 (MQ=255)

                             gCCATAAGGCGCTGCGACAGTTCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGc                                                            <  2:4198362/93‑1 (MQ=255)

                             gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGc                                                            >  1:7102531/1‑93 (MQ=255)

                             gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGc                                                            >  2:3856531/1‑93 (MQ=255)

                             gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCATTGAGCGAATAAACAACTTGTGAAGc                                                            >  2:4881182/1‑93 (MQ=255)

                              ccATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGcc                                                           >  2:6920063/1‑93 (MQ=255)

                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGc                                                        >  2:5346620/1‑93 (MQ=255)

                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGc                                                        >  1:5223988/1‑93 (MQ=255)

                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGc                                                        >  1:8508949/1‑93 (MQ=255)

                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGGGAAGCCAGc                                                        >  2:8706321/1‑93 (MQ=255)

                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAAAAAACAACTTGTGAAGCCAGc                                                        >  2:7805453/1‑93 (MQ=255)

                                  aaGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCa                                                       >  2:2326169/1‑93 (MQ=255)

                                  aaGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCa                                                       <  1:8770669/93‑1 (MQ=255)

                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      >  1:3784652/1‑93 (MQ=255)

                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      <  2:6490829/93‑1 (MQ=255)

                                    ggCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATa                                                     >  2:6393936/1‑93 (MQ=255)

                                    ggCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATa                                                     >  1:7762396/1‑93 (MQ=255)

                                    ggCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATa                                                     >  1:2410565/1‑93 (MQ=255)

                                     gcgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaa                                                    <  1:338758/93‑1 (MQ=255)

                                      cgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  1:6341485/93‑1 (MQ=255)

                                      cgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  1:970762/93‑1 (MQ=255)

                                       ggtGCGACAGTGACGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAAc                                                  <  2:8056380/91‑1 (MQ=255)

                                       gcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaag                                                  >  1:7900935/1‑92 (MQ=255)

                                        cTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGCAAAACCCGCGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACg                                                 >  2:3065486/1‑93 (MQ=255)

                                        cTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACg                                                 >  2:3379668/1‑93 (MQ=255)

                                         tGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACgg                                                >  2:7333028/1‑93 (MQ=255)

                                         tGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACgg                                                <  2:8587264/93‑1 (MQ=255)

                                           cGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGAt                                              <  1:3002409/93‑1 (MQ=255)

                                              cAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGc                                           >  1:2852854/1‑93 (MQ=255)

                                               aGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGCGCGAATAAACAACTTGTGAAGCACGCATAAACGAATaaca                                          >  1:2312880/1‑90 (MQ=255)

                                                gTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCg                                         >  2:7761652/1‑93 (MQ=255)

                                                gTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCg                                         >  2:6770486/1‑93 (MQ=255)

                                                 tGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGt                                        >  1:9507105/1‑93 (MQ=255)

                                                  gCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTa                                       >  1:3558086/1‑93 (MQ=255)

                                                   ccGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTaa                                      <  2:200260/93‑1 (MQ=255)

                                                   ccGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTaa                                      <  1:1181463/93‑1 (MQ=255)

                                                    cGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAAc                                     <  2:108999/93‑1 (MQ=255)

                                                    cGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAAc                                     <  1:4520232/93‑1 (MQ=255)

                                                     gATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACt                                    >  1:830763/1‑93 (MQ=255)

                                                       tatTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  <  1:6573325/93‑1 (MQ=255)

                                                        atTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTggg                                 <  2:6144579/93‑1 (MQ=255)

                                                        atTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTggg                                 <  1:7327086/93‑1 (MQ=255)

                                                         ttCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTCACTgggg                                >  1:5736813/1‑92 (MQ=255)

                                                         ttCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggc                                <  1:2630254/93‑1 (MQ=255)

                                                         ttCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggc                                <  1:8652787/93‑1 (MQ=255)

                                                           cTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGCTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggcgg                              <  2:9301859/93‑1 (MQ=255)

                                                            tGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggcggc                             <  2:9117687/93‑1 (MQ=255)

                                                            tGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggcggc                             <  1:1074265/93‑1 (MQ=255)

                                                             tgCTTGCGATGTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCa                            <  1:6449633/92‑1 (MQ=255)

                                                              gCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggcggccg                           >  2:10018207/1‑91 (MQ=255)

                                                               cTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGa                          >  1:6027161/1‑93 (MQ=255)

                                                                 tGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaaa                        >  1:9239906/1‑93 (MQ=255)

                                                                    gATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGc                     <  1:549538/93‑1 (MQ=255)

                                                                    gATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGc                     <  2:7613310/93‑1 (MQ=255)

                                                                    gATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGc                     <  1:9155520/93‑1 (MQ=255)

                                                                     aTCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCt                    <  2:1753775/93‑1 (MQ=255)

                                                                         cAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCa                <  1:5118349/93‑1 (MQ=255)

                                                                         cAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCa                <  1:743482/93‑1 (MQ=255)

                                                                          aaCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAg               <  2:5600608/93‑1 (MQ=255)

                                                                          aaCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAAAGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAg               <  1:1837278/93‑1 (MQ=255)

                                                                            cGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAggg             <  2:5032704/93‑1 (MQ=255)

                                                                             ggAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGt            >  2:4055230/1‑93 (MQ=255)

                                                                               aaaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTaa          <  1:7295901/93‑1 (MQ=255)

                                                                               aaaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTaa          <  1:4238321/93‑1 (MQ=255)

                                                                                 aGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACg        >  1:3434562/1‑93 (MQ=255)

                                                                                  gACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGc       <  2:2059015/93‑1 (MQ=255)

                                                                                     ccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtgg    <  2:4759812/93‑1 (MQ=255)

                                                                                      cgtgcgtgCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggt   >  2:1184114/1‑93 (MQ=255)

                                                                                      cgtgcgtgCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggt   <  1:900290/93‑1 (MQ=255)

                                                                                      cgtgcgtgCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggg   >  2:3697988/1‑92 (MQ=255)

                                                                                       gtgcgtgcGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggtg  <  1:4652496/93‑1 (MQ=255)

TTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTG  >  NC_000913/4067458‑4067637

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