BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,608,034	A→G	V224A (GTG→GCG)	uxaB ←	altronate oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,608,034	0	A	G	100.0%	886.0 / NA	250	V224A (GTG→GCG)	uxaB	altronate oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (131/119); total (131/119)

Only 100 of 255 total aligned reads displayed.

GATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCT  >  NC_000913/1607944‑1608123

 aTTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCtt                                                                                        >  2:6292975/1‑93 (MQ=255)

 aTTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCtt                                                                                        <  1:261978/93‑1 (MQ=255)

  tttCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTc                                                                                       <  2:5135757/93‑1 (MQ=255)

   ttCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCa                                                                                      >  1:6137710/1‑93 (MQ=255)

     cGGTCCCTGAATCACAAACAGGTAAAAGTGTTCATCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATc                                                                                    <  2:6105813/93‑1 (MQ=255)

       gTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATcgc                                                                                  <  1:3853025/93‑1 (MQ=255)

       gTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATcgc                                                                                  >  1:3093109/1‑93 (MQ=255)

          ccTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATcgcgcg                                                                               >  2:227926/1‑93 (MQ=255)

            tGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGa                                                                             >  2:3775569/1‑93 (MQ=255)

             gAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGAt                                                                            >  2:929426/1‑93 (MQ=255)

              aaTCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATa                                                                           >  1:3013998/1‑93 (MQ=255)

                  acaAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCg                                                                       >  2:1872060/1‑93 (MQ=255)

                  acaAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCg                                                                       >  1:5695538/1‑93 (MQ=255)

                   caAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCgg                                                                      >  1:6457602/1‑93 (MQ=255)

                   caAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCgg                                                                      >  2:6394903/1‑93 (MQ=255)

                   caAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCgg                                                                      >  1:4149225/1‑93 (MQ=255)

                    aaaCAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGt                                                                     >  1:2870009/1‑93 (MQ=255)

                     aaCAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTa                                                                    >  2:462384/1‑93 (MQ=255)

                     aaCAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTa                                                                    >  2:3830836/1‑93 (MQ=255)

                      aCAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTaa                                                                   <  1:1012599/93‑1 (MQ=255)

                        aGGTAAAATTGTTCAGCGGGGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACg                                                                 >  1:6323625/1‑93 (MQ=255)

                        aGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACg                                                                 >  1:2126269/1‑93 (MQ=255)

                        aGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACg                                                                 >  2:1681457/1‑93 (MQ=255)

                         ggTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGa                                                                >  2:5364447/1‑93 (MQ=255)

                          gTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGAt                                                               <  2:3123125/93‑1 (MQ=255)

                            aaaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAATCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATAc                                                             <  1:555868/93‑1 (MQ=255)

                            aaaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATAc                                                             <  2:3486297/93‑1 (MQ=255)

                            aaaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATAc                                                             >  2:5764979/1‑93 (MQ=255)

                            aaaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATAc                                                             <  1:548036/93‑1 (MQ=255)

                             aaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACg                                                            >  2:3235987/1‑93 (MQ=255)

                             aaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACg                                                            >  1:4553009/1‑93 (MQ=255)

                             aaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACg                                                            >  1:4498910/1‑93 (MQ=255)

                             aaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACg                                                            <  2:3607338/93‑1 (MQ=255)

                              aaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACgg                                                           >  1:2768280/1‑93 (MQ=255)

                              aaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACgg                                                           >  2:1496315/1‑93 (MQ=255)

                               aGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGt                                                          >  2:217621/1‑93 (MQ=255)

                               aGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGt                                                          >  2:5917662/1‑93 (MQ=255)

                               aGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGt                                                          >  1:2227486/1‑93 (MQ=255)

                               aGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGAGGATAACCGGTAACGATACggc                                                          >  2:1998855/1‑92 (MQ=255)

                               aGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTGCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGt                                                          <  2:5860899/93‑1 (MQ=255)

                                gtgtTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTc                                                         >  2:2357930/1‑93 (MQ=255)

                                gtgtTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTc                                                         <  2:4772055/93‑1 (MQ=255)

                                 tgtTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTcc                                                        >  1:2497098/1‑93 (MQ=255)

                                   ttCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTccac                                                      <  2:5249848/93‑1 (MQ=255)

                                   ttCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTccaa                                                      >  2:1843312/1‑92 (MQ=255)

                                     cAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTccacca                                                    <  2:5668367/93‑1 (MQ=255)

                                      aGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTGATCGCGCGGATAACCGGTAACGATACGGTccacaac                                                   >  2:1590744/1‑90 (MQ=255)

                                       gCGGTGTCGAGAAAACCATCGTGATAACCCGGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGc                                                  >  2:705329/1‑93 (MQ=255)

                                       gCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGc                                                  <  2:2781801/93‑1 (MQ=255)

                                         ggTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGt                                                <  2:3370427/93‑1 (MQ=255)

                                          gtgtCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTa                                               >  1:1052614/1‑93 (MQ=255)

                                              cGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAAc                                           <  2:2521680/93‑1 (MQ=255)

                                               gagaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACa                                          >  2:3176050/1‑93 (MQ=255)

                                                agaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAg                                         >  2:726514/1‑93 (MQ=255)

                                                agaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAg                                         >  1:4192531/1‑93 (MQ=255)

                                                  aaaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGaa                                       >  2:1461625/1‑93 (MQ=255)

                                                  aaaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGaa                                       <  2:623132/93‑1 (MQ=255)

                                                  aaaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCAGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGaa                                       <  2:1702121/93‑1 (MQ=255)

                                                   aaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAg                                      >  2:1013608/1‑93 (MQ=255)

                                                   aaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAg                                      <  1:17135/93‑1 (MQ=255)

                                                   aaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAg                                      <  2:4958522/93‑1 (MQ=255)

                                                   aaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAg                                      >  1:3876/1‑93 (MQ=255)

                                                    aaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGc                                     <  2:1432446/93‑1 (MQ=255)

                                                     aCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCt                                    <  1:5316478/93‑1 (MQ=255)

                                                     aCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCt                                    >  1:159604/1‑93 (MQ=255)

                                                     aCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCt                                    <  1:2054871/93‑1 (MQ=255)

                                                     aCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCt                                    >  1:4857837/1‑93 (MQ=255)

                                                      ccATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTg                                   <  1:6469868/93‑1 (MQ=255)

                                                      ccATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTg                                   <  1:6343315/93‑1 (MQ=255)

                                                      ccATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTg                                   >  2:4457839/1‑93 (MQ=255)

                                                       cATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGt                                  <  2:4215596/93‑1 (MQ=255)

                                                       cATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGt                                  <  1:1528634/93‑1 (MQ=255)

                                                        aTCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGtt                                 >  1:3418614/1‑93 (MQ=255)

                                                           gTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGc                              >  2:1079177/1‑93 (MQ=255)

                                                            tGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGcc                             <  1:1195144/93‑1 (MQ=255)

                                                             gATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCt                            >  2:4376973/1‑93 (MQ=255)

                                                             gATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCt                            <  1:2547815/93‑1 (MQ=255)

                                                                aaCCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAt                         <  2:5398727/93‑1 (MQ=255)

                                                                aaCCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAt                         <  1:3647477/93‑1 (MQ=255)

                                                                  cccAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATcc                       <  2:5621935/93‑1 (MQ=255)

                                                                    cAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAg                     >  2:1315239/1‑93 (MQ=255)

                                                                    cAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAg                     >  2:4721454/1‑93 (MQ=255)

                                                                    cAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAg                     >  2:5822669/1‑93 (MQ=255)

                                                                     aGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGc                    <  1:283054/93‑1 (MQ=255)

                                                                      gTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGcc                   >  2:2686789/1‑93 (MQ=255)

                                                                      gTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGcc                   <  2:4103430/93‑1 (MQ=255)

                                                                        tcttcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCAt                 <  1:3292280/93‑1 (MQ=255)

                                                                         cttcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCAtt                <  1:1896657/93‑1 (MQ=255)

                                                                            cttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGaa             >  1:1317290/1‑93 (MQ=255)

                                                                              tcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATa           <  2:4006595/93‑1 (MQ=255)

                                                                               ccAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATaa          >  1:3622612/1‑93 (MQ=255)

                                                                                cAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATaaa         <  1:3529391/93‑1 (MQ=255)

                                                                                cAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATaaa         <  1:1724451/93‑1 (MQ=255)

                                                                                cAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATaaa         <  1:1682458/93‑1 (MQ=255)

                                                                                   tttAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGc      <  1:3796795/93‑1 (MQ=255)

                                                                                    ttAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCt     <  2:2990210/93‑1 (MQ=255)

                                                                                     tAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCtt    <  2:503632/93‑1 (MQ=255)

                                                                                      aGCCGCTTCATCGCGCGGATAACCGGTAAGGATACGGTCCACCAGCGTAGAACAGAAGCTGTTTGCCTGATCCAGCCATTGAATAAATGCTTc   >  2:687349/1‑93 (MQ=255)

                                                                                      aGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTc   <  2:3015760/93‑1 (MQ=255)

                                                                                       gCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCt  >  1:5880384/1‑93 (MQ=255)

GATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCT  >  NC_000913/1607944‑1608123

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