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%	1259.0 / NA	393	I198I (ATT→ATC)	hemC	porphobilinogen deaminase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (203/190); total (203/190)

Only 100 of 418 total aligned reads displayed.

GTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC  >  NC_000913/3987717‑3987896

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

 tgATTCAGCGCGTCAAGCAGCTCGCGAGTGTGTGAATCATCAAGGCGTCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtc                                                                                        <  1:4862273/93‑1 (MQ=255)

  gATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtct                                                                                       <  2:456495/93‑1 (MQ=255)

   aTTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtctc                                                                                      <  2:6219252/93‑1 (MQ=255)

   aTTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtctc                                                                                      <  1:952866/93‑1 (MQ=255)

    ttCAGCGCGGCAAGCAGCTCGCGAGTGCGTTAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCg                                                                                     <  1:5514055/93‑1 (MQ=255)

    ttCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCg                                                                                     <  1:4556733/93‑1 (MQ=255)

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

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

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

         cgcgGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTgg                                                                                <  2:828868/93‑1 (MQ=255)

          gcgGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGc                                                                               >  1:6342479/1‑93 (MQ=255)

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

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

                   agcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACgcggcg                                                                      >  1:5650593/1‑93 (MQ=255)

                   agcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACgcggcg                                                                      >  1:7592102/1‑93 (MQ=255)

                   agcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACgcggcg                                                                      >  1:1487119/1‑93 (MQ=255)

                     cTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGcgcg                                                                    <  2:5660720/93‑1 (MQ=255)

                      tCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGa                                                                   <  1:2904600/93‑1 (MQ=255)

                       cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGaa                                                                  <  2:4864375/93‑1 (MQ=255)

                       cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGaa                                                                  >  1:6927059/1‑93 (MQ=255)

                        gcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGTGTGGCAACTCGGCGCGagt                                                                 >  2:5076649/1‑91 (MQ=255)

                         cgcgTGCGTGAATCATCAAGGCAGCATTCACTGCCCCCCGCCCCTTGTCCTTCCGCCGGAAGAGAGATCTCGCGTGGCAACGCGGCGCGAAtc                                                                <  2:6835267/90‑2 (MQ=255)

                         cgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATa                                                                <  1:30064/93‑1 (MQ=255)

                           aGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACg                                                              >  1:3347260/1‑93 (MQ=255)

                              gCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGa                                                           <  1:3886613/93‑1 (MQ=255)

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

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

                                 tGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTc                                                        <  1:6540343/93‑1 (MQ=255)

                                  gAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcc                                                       <  2:1444019/93‑1 (MQ=255)

                                   aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGCGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                      <  1:619217/93‑1 (MQ=255)

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

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

                                     tcatcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGa                                                    >  2:3969569/1‑93 (MQ=255)

                                       ctcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGAcc                                                  >  1:3751737/2‑93 (MQ=255)

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

                                        tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                 <  1:1233878/93‑1 (MQ=255)

                                        tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                 <  1:1014727/93‑1 (MQ=255)

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

                                             gcggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACg                                            <  2:3367003/93‑1 (MQ=255)

                                              cggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGt                                           <  2:2676583/93‑1 (MQ=255)

                                               ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                          <  2:7182801/93‑1 (MQ=255)

                                               ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                          <  2:6636110/93‑1 (MQ=255)

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

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

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

                                                   ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                      >  1:2239873/1‑93 (MQ=255)

                                                   ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                      >  2:4479438/1‑93 (MQ=255)

                                                   ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                      >  1:1772362/1‑93 (MQ=255)

                                                   ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGGTTTAg                                      >  2:4662943/1‑93 (MQ=255)

                                                    tCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGt                                     <  1:914034/93‑1 (MQ=255)

                                                     cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTTGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                    <  2:4804457/93‑1 (MQ=255)

                                                     cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                    >  2:3223250/1‑93 (MQ=255)

                                                     cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                    <  1:2763690/93‑1 (MQ=255)

                                                     cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                    <  2:689526/93‑1 (MQ=255)

                                                     cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                    <  2:1310590/93‑1 (MQ=255)

                                                     cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                    >  1:4367620/1‑93 (MQ=255)

                                                      aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcc                                   <  1:154955/93‑1 (MQ=255)

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

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

                                                        tACCCACCGCACCTTGTCCTACCACCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCgg                                 >  2:1844802/1‑93 (MQ=255)

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

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

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

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

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

                                                           ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGCCTCCAGAACTAAACGTTTTAGTCCGGCTa                              <  1:1410806/93‑1 (MQ=255)

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

                                                              ccGCACCTTGTCCTACCGCCGGAAGAGAGATCTGGGTTGGCAACGCGGCGCGAATACGTGACTCCGGCGCTAAAGGTTTTAGTCCGGCTACgg                           >  1:2067049/1‑93 (MQ=255)

                                                              ccGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACgg                           >  1:5745294/1‑93 (MQ=255)

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

                                                                gCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                         >  1:6642827/1‑93 (MQ=255)

                                                                gCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                         >  1:7760245/1‑93 (MQ=255)

                                                                cgcccTTGGTCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                         <  2:4500783/90‑1 (MQ=255)

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

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

                                                                  aCCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTCTAGTCCGGCTACGGCAAg                       >  1:1886974/1‑93 (MQ=255)

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

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

                                                                   ccTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                      <  1:312467/93‑1 (MQ=255)

                                                                   ccTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCAGCTACGGCAAGa                      <  2:1923805/93‑1 (MQ=255)

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

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

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

                                                                       gTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatga                  <  2:2486645/93‑1 (MQ=255)

                                                                        tCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                 >  1:4793894/1‑93 (MQ=255)

                                                                          cTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATgg               >  2:5441472/1‑93 (MQ=255)

                                                                               gccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCg          >  2:7874049/1‑93 (MQ=255)

                                                                                ccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGt         <  1:6993158/93‑1 (MQ=255)

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

                                                                                  ggAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt       >  1:2035750/1‑93 (MQ=255)

                                                                                  ggAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt       >  2:7694136/1‑93 (MQ=255)

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

                                                                                  ggAAGAGAGATCTCGGGTGGCAACGCGGAGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt       >  2:2149077/1‑93 (MQ=255)

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

                                                                                     agagagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg    <  2:5758086/93‑1 (MQ=255)

                                                                                     agagagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg    <  1:5349816/93‑1 (MQ=255)

                                                                                      gagagaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGc   <  2:7417380/93‑1 (MQ=255)

                                                                                       agagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGcc  <  1:879144/93‑1 (MQ=255)

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

GTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC  >  NC_000913/3987717‑3987896

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