BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,987,807	A→G	I198I (ATT→ATC)	hemC ←	porphobilinogen deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,987,807	0	A	G	100.0%	860.2 / NA	264	I198I (ATT→ATC)	hemC	porphobilinogen deaminase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (142/122); total (142/122)

Only 100 of 278 total aligned reads displayed.

GTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC  >  NC_000913/3987717‑3987896

gtgATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAt                                                                                         >  1:5809399/1‑93 (MQ=255)

gtgATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAt                                                                                         <  1:5401100/93‑1 (MQ=255)

     tCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCgg                                                                                    <  1:4780021/93‑1 (MQ=255)

      cAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCggg                                                                                   <  1:6181466/93‑1 (MQ=255)

       aGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGt                                                                                  <  2:3361679/93‑1 (MQ=255)

       aGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGt                                                                                  >  1:2205724/1‑93 (MQ=255)

         cgcgGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTgg                                                                                >  2:2313042/1‑93 (MQ=255)

          gcgGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATGCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGc                                                                               <  2:855942/93‑1 (MQ=255)

            ggCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCaa                                                                             <  2:2775968/93‑1 (MQ=255)

            ggCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCaa                                                                             <  2:1709424/93‑1 (MQ=255)

             gCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCaac                                                                            <  1:1535976/93‑2 (MQ=255)

              cAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCaacg                                                                           >  1:1060115/1‑91 (MQ=255)

               aaGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCaacgc                                                                          >  2:1245272/1‑90 (MQ=255)

                agcagcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCaacgcg                                                                         >  2:805456/1‑89 (MQ=255)

                  cagcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGgc                                                                       >  2:1580107/1‑93 (MQ=255)

                    gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGgcgc                                                                     >  2:644175/1‑93 (MQ=255)

                    gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGgcgc                                                                     >  1:5111678/1‑93 (MQ=255)

                    gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGgcgc                                                                     >  2:1995736/1‑93 (MQ=255)

                        gcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAAt                                                                 <  1:4924000/93‑1 (MQ=255)

                        gcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAAt                                                                 >  2:4287521/1‑93 (MQ=255)

                          gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATAc                                                               <  1:3999724/93‑1 (MQ=255)

                          gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATAc                                                               <  1:1916178/93‑1 (MQ=255)

                           aGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACg                                                              <  2:690104/93‑1 (MQ=255)

                            gTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGt                                                             >  2:4466365/1‑93 (MQ=255)

                              gCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGa                                                           >  1:2060285/1‑93 (MQ=255)

                               cGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGAc                                                          <  2:2598591/93‑1 (MQ=255)

                                gTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACt                                                         >  2:3942365/1‑93 (MQ=255)

                                 tGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTc                                                        >  2:5415010/1‑93 (MQ=255)

                                  gAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcc                                                       >  1:4040457/1‑93 (MQ=255)

                                  gAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcc                                                       >  1:4475814/1‑93 (MQ=255)

                                   aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                      >  2:3443300/1‑93 (MQ=255)

                                   aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                      >  2:1805661/1‑93 (MQ=255)

                                   aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                      >  1:5033107/1‑93 (MQ=255)

                                    atcatcAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAg                                                     >  2:3874291/1‑93 (MQ=255)

                                    atcatcAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAg                                                     <  1:672221/93‑1 (MQ=255)

                                     tcatcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGa                                                    <  2:3692103/93‑1 (MQ=255)

                                      catcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCCGAc                                                   >  2:5372661/1‑93 (MQ=255)

                                        tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                 >  2:5835926/1‑93 (MQ=255)

                                        tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                 >  1:3846189/1‑93 (MQ=255)

                                          aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                               >  2:963260/1‑93 (MQ=255)

                                          aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                               <  1:4594173/93‑1 (MQ=255)

                                          aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                               <  1:6154800/93‑1 (MQ=255)

                                           aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                              >  2:3070350/1‑93 (MQ=255)

                                               ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                          >  2:255008/1‑93 (MQ=255)

                                                 cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                        <  1:6032217/93‑1 (MQ=255)

                                                 cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                        <  1:4607049/93‑1 (MQ=255)

                                                  aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                       <  1:3322776/93‑1 (MQ=255)

                                                  aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                       >  1:233232/1‑93 (MQ=255)

                                                      aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcc                                   >  2:5989731/1‑93 (MQ=255)

                                                       aTACCCACCGCACCTTGTCCTACCGCCGGGAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAAGTTTTAGTGCg                                  >  2:9160/1‑93 (MQ=255)

                                                       aTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCg                                  >  1:4313554/1‑93 (MQ=255)

                                                       aTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCg                                  >  2:609734/1‑93 (MQ=255)

                                                        tACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCgg                                 >  1:1013719/1‑93 (MQ=255)

                                                        tACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCgg                                 >  1:197812/1‑93 (MQ=255)

                                                         aCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGc                                <  1:3079712/93‑1 (MQ=255)

                                                          cccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCt                               >  2:3281580/1‑93 (MQ=255)

                                                           ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                              >  2:5166744/1‑93 (MQ=255)

                                                           ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                              <  1:4123132/93‑1 (MQ=255)

                                                           ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                              >  1:2551699/1‑93 (MQ=255)

                                                           ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                              <  1:1155048/93‑1 (MQ=255)

                                                             aCCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACg                            <  1:5980786/93‑1 (MQ=255)

                                                             aCCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACg                            >  1:3513340/1‑93 (MQ=255)

                                                             aCCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACg                            >  1:751193/1‑93 (MQ=255)

                                                             aCCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACg                            >  1:2830434/1‑93 (MQ=255)

                                                               cGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGc                          >  1:4739278/1‑93 (MQ=255)

                                                                gCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                         <  2:2103784/93‑1 (MQ=255)

                                                                  aCCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAg                       >  1:1643520/1‑93 (MQ=255)

                                                                  aCCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAg                       >  2:4560441/1‑93 (MQ=255)

                                                                  aCCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAg                       >  2:2952143/1‑93 (MQ=255)

                                                                  aCCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAg                       >  1:4195366/1‑93 (MQ=255)

                                                                   ccTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                      >  2:2186285/1‑93 (MQ=255)

                                                                   ccTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                      >  1:2375626/1‑93 (MQ=255)

                                                                   ccTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                      >  2:2364281/1‑93 (MQ=255)

                                                                    cTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                     >  1:908895/1‑93 (MQ=255)

                                                                    cTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                     <  2:1367326/93‑1 (MQ=255)

                                                                     ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                    >  1:3261418/1‑93 (MQ=255)

                                                                     ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                    <  2:4973704/93‑1 (MQ=255)

                                                                     ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                    <  1:3756248/93‑1 (MQ=255)

                                                                     ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                    >  1:10076/1‑93 (MQ=255)

                                                                      tGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                   <  1:3877389/93‑1 (MQ=255)

                                                                          cTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATgg               <  1:571704/93‑1 (MQ=255)

                                                                           tACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc              <  2:5805822/93‑1 (MQ=255)

                                                                           tACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc              <  2:1632911/93‑1 (MQ=255)

                                                                             ccgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCAt            <  2:4201502/93‑1 (MQ=255)

                                                                             ccgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCAt            <  2:4338027/93‑1 (MQ=255)

                                                                              cgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc           <  2:5658218/93‑1 (MQ=255)

                                                                              cgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc           <  1:1704210/93‑1 (MQ=255)

                                                                                 cgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTa        >  2:4870153/1‑93 (MQ=255)

                                                                                 cgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTa        <  2:3230746/93‑1 (MQ=255)

                                                                                 cgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTa        <  1:831436/93‑1 (MQ=255)

                                                                                  ggAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt       <  1:4239295/93‑1 (MQ=255)

                                                                                   gAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt      >  2:4915663/1‑93 (MQ=255)

                                                                                   gAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt      >  2:542198/1‑93 (MQ=255)

                                                                                   gAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt      >  1:5910852/1‑93 (MQ=255)

                                                                                    aaGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTc     <  2:1122517/93‑1 (MQ=255)

                                                                                     agagagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg    >  2:3647958/1‑93 (MQ=255)

                                                                                     agagagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg    >  1:1373766/1‑93 (MQ=255)

                                                                                      gagagaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTctc   <  2:5540592/93‑3 (MQ=255)

                                                                                      gagagaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTctc   <  1:2688392/93‑3 (MQ=255)

                                                                                       agagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGcc  >  1:515401/1‑93 (MQ=255)

GTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC  >  NC_000913/3987717‑3987896

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