BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,581,520	C→T	S510N (AGC→AAC)	hsdM ←	type I restriction enzyme EcoKI methyltransferase component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,581,520	0	C	T	100.0%	1054.4 / NA	247	S510N (AGC→AAC)	hsdM	type I restriction enzyme EcoKI methyltransferase component
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (130/117); total (130/117)

Only 100 of 252 total aligned reads displayed.

GGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGC  >  NC_000913/4581422‑4581616

    cGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGccc                                                                                            >  2:6802557/1‑101 (MQ=255)

    cGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGccc                                                                                            <  2:6802221/101‑1 (MQ=255)

     gATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGcccc                                                                                           >  1:6802559/1‑101 (MQ=255)

     gATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGcccc                                                                                           <  1:6802247/101‑1 (MQ=255)

      aTAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGccccc                                                                                          <  2:6802102/101‑1 (MQ=255)

      aTAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGccccc                                                                                          >  2:6802563/1‑101 (MQ=255)

       tAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCa                                                                                         >  1:6802565/1‑101 (MQ=255)

        aaCCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAg                                                                                        >  2:6802568/1‑101 (MQ=255)

         acccacccCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGt                                                                                       >  2:6802571/1‑101 (MQ=255)

           ccacccCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTc                                                                                     >  1:6802575/1‑101 (MQ=255)

            cacccCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCa                                                                                    >  2:6802576/1‑101 (MQ=255)

              ccccTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACg                                                                                  >  2:6802579/1‑101 (MQ=255)

                ccTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCa                                                                                >  1:6802581/1‑101 (MQ=255)

                 cTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCAt                                                                               <  2:6802056/101‑1 (MQ=255)

                   ccGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCa                                                                             <  1:6802026/101‑1 (MQ=255)

                     ggCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAg                                                                                        <  2:6801931/88‑1 (MQ=255)

                      gCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAgcg                                                                          >  1:6802584/1‑101 (MQ=255)

                       cAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAgcgc                                                                         >  1:6802587/1‑101 (MQ=255)

                         aTTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCAt                                                                       >  1:6802590/1‑101 (MQ=255)

                            tCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCa                                                                    <  2:6802327/101‑1 (MQ=255)

                            tCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCa                                                                    >  2:6802593/1‑101 (MQ=255)

                            tCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCa                                                                    <  2:6802087/101‑1 (MQ=255)

                              cccGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGt                                                                  <  2:6802216/101‑1 (MQ=255)

                              cccGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGt                                                                  >  2:6802596/1‑101 (MQ=255)

                              cccGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGt                                                                  <  1:6802305/101‑1 (MQ=255)

                               ccGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGtt                                                                 >  2:6802598/1‑101 (MQ=255)

                                  cACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAg                                                              <  1:6802308/101‑1 (MQ=255)

                                   aCTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGa                                                             >  1:6802601/1‑101 (MQ=255)

                                   aCTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGACCCCAGTTCACGCATCAGCGCATCCAGTTCAGa                                                             >  2:6802604/1‑101 (MQ=255)

                                    cTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGAc                                                            >  2:6802606/1‑101 (MQ=255)

                                       aTTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAgc                                                         >  1:6802609/1‑101 (MQ=255)

                                        ttCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAgcg                                                        <  2:6802332/101‑1 (MQ=255)

                                          ccTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGcc                                                      >  1:6802613/1‑101 (MQ=255)

                                          ccTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGcc                                                      >  2:6802612/1‑101 (MQ=255)

                                            ttCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTg                                                    <  2:6802368/101‑1 (MQ=255)

                                            ttCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTg                                                    >  2:6802616/1‑101 (MQ=255)

                                             tCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGt                                                   >  1:6802617/1‑101 (MQ=255)

                                              cgccccACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGtt                                                                                      <  2:6802233/63‑1 (MQ=255)

                                              ccccccACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTa                                                  <  2:6801968/99‑1 (MQ=255)

                                              cACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTa                                                  <  1:6802135/101‑1 (MQ=255)

                                                  ccaccaAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAg                                              <  1:6802033/101‑1 (MQ=255)

                                                  ccaccaAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAg                                              >  1:6802622/1‑101 (MQ=255)

                                                   caccaAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGt                                             >  1:6802624/1‑101 (MQ=255)

                                                   caccaAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGt                                             <  2:6802590/101‑1 (MQ=255)

                                                    accaAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGtt                                            <  2:6802208/101‑1 (MQ=255)

                                                     ccaAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTc                                           <  1:6802387/101‑1 (MQ=255)

                                                     ccaAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTc                                           >  1:6802627/1‑101 (MQ=255)

                                                      caAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCg                                          <  1:6802199/101‑1 (MQ=255)

                                                      caAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCg                                          <  2:6802226/101‑1 (MQ=255)

                                                       aaaCGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGc                                         >  1:6802630/1‑101 (MQ=255)

                                                       aaaCGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGc                                         >  2:6802629/1‑101 (MQ=255)

                                                         aCGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGccc                                       <  1:6802414/101‑1 (MQ=255)

                                                          cGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCa                                      >  2:6802633/1‑101 (MQ=255)

                                                           gCTTCTTCCAGCAACTGACGCTGCAAACCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCAt                                     <  2:6802200/101‑1 (MQ=255)

                                                            cttcttCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATc                                    <  1:6802324/101‑1 (MQ=255)

                                                            cttcttCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATc                                    <  2:6802137/101‑1 (MQ=255)

                                                            cttcttCCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATc                                    <  2:6802130/101‑1 (MQ=255)

                                                             ttcttcCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCTCCCATCg                                   <  1:6802188/101‑1 (MQ=255)

                                                              tcttcCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGc                                  >  2:6802637/1‑101 (MQ=255)

                                                               cttcCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCt                                 >  2:6802640/1‑101 (MQ=255)

                                                               cttcCAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCt                                 >  1:6802639/1‑101 (MQ=255)

                                                                  ccAGCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCt                              <  2:6802112/101‑1 (MQ=255)

                                                                     gCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGcc                           <  1:6802206/101‑1 (MQ=255)

                                                                     gCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGcc                           >  1:6802645/1‑101 (MQ=255)

                                                                     gCAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGcc                           >  2:6802644/1‑101 (MQ=255)

                                                                      cAACTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCg                          <  2:6802246/101‑1 (MQ=255)

                                                                       aaCTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGc                         >  2:6802650/1‑101 (MQ=255)

                                                                       aaCTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGc                         >  1:6802649/1‑101 (MQ=255)

                                                                        aCTGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCt                        <  1:6802144/101‑1 (MQ=255)

                                                                          tGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTaa                      >  2:6802654/1‑101 (MQ=255)

                                                                          tGACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTaa                      <  1:6802195/101‑1 (MQ=255)

                                                                           gACGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAAt                     >  1:6802656/1‑101 (MQ=255)

                                                                            aCGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATa                    >  2:6802659/1‑101 (MQ=255)

                                                                             cGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATAc                   >  2:6802662/1‑101 (MQ=255)

                                                                             cGCTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATAc                   <  1:6802172/101‑1 (MQ=255)

                                                                               cTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGtt                                                                 <  2:6802479/53‑1 (MQ=255)

                                                                               cTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAt                 <  2:6802380/101‑1 (MQ=255)

                                                                               cTGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAt                 <  1:6802281/101‑1 (MQ=255)

                                                                                tGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtc                <  1:6802261/101‑1 (MQ=255)

                                                                                tGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtc                >  2:6802664/1‑101 (MQ=255)

                                                                                tGCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtc                <  1:6801903/101‑1 (MQ=255)

                                                                                 gCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtcc               <  1:6802168/101‑1 (MQ=255)

                                                                                 gCAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtcc               >  2:6802667/1‑101 (MQ=255)

                                                                                  cAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCtt                                <  1:6801970/83‑1 (MQ=255)

                                                                                  cAAATCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccg              >  1:6802668/1‑101 (MQ=255)

                                                                                   aaaTCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccgg             >  1:6802671/1‑101 (MQ=255)

                                                                                   aaaTCGGCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccgg             >  2:6802672/1‑101 (MQ=255)

                                                                                        ggCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccggctccg        <  2:6802302/101‑1 (MQ=255)

                                                                                        ggCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccggctccg        >  1:6802676/1‑101 (MQ=255)

                                                                                        ggCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccggctccg        <  1:6802089/101‑1 (MQ=255)

                                                                                         gCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccggctccgg       >  1:6802678/1‑101 (MQ=255)

                                                                                         gCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccggctccgg       <  1:6802300/101‑1 (MQ=255)

                                                                                         gCCTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccggctccgg       >  2:6802680/1‑101 (MQ=255)

                                                                                          ccTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACAtccggctccggc      >  1:6802683/1‑101 (MQ=255)

                                                                                           cTCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCa     <  2:6802285/101‑1 (MQ=255)

                                                                                            tCATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAg    <  1:6802104/101‑1 (MQ=255)

                                                                                             cATCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAgg   >  1:6802686/1‑101 (MQ=255)

                                                                                              aTCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGc  <  2:6802241/101‑1 (MQ=255)

                                                                                              aTCGTTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGc  <  1:6802229/101‑1 (MQ=255)

                                                                                               tcgttcgcCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCg                                          <  2:6802079/60‑1 (MQ=255)

GGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGC  >  NC_000913/4581422‑4581616

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