BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,274,907	A→G	G45G (GGT→GGC)	yjcB ←	uncharacterized protein YjcB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,274,907	0	A	G	100.0%	1251.7 / NA	296	G45G (GGT→GGC)	yjcB	uncharacterized protein YjcB
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (153/143); total (153/143)

Only 100 of 296 total aligned reads displayed.

ATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCC  >  NC_000913/4274809‑4275004

atcatcACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGcc                                                                                                 >  1:6677838/1‑101 (MQ=255)

 tcatcaCCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCg                                                                                                <  1:6677082/101‑1 (MQ=255)

  catcaCCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCgc                                                                                               >  2:6677841/1‑101 (MQ=255)

    tcaCCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCgctg                                                                                             >  2:6677845/1‑101 (MQ=255)

     caCCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCgctgc                                                                                            >  1:6677846/1‑101 (MQ=255)

      aCCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCgctgct                                                                                           >  1:6677850/1‑101 (MQ=255)

      aCCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCgctgct                                                                                           <  1:6677470/101‑1 (MQ=255)

       ccGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCgctgct                                                                                           <  2:6677339/100‑1 (MQ=255)

       ccGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGctgctg                                                                                          >  1:6677853/1‑101 (MQ=255)

        cGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGa                                                                                         <  1:6677432/101‑1 (MQ=255)

         gACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGAt                                                                                        <  2:6677468/101‑1 (MQ=255)

         gACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGAt                                                                                        >  2:6677856/1‑101 (MQ=255)

          aCTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATc                                                                                       >  2:6677858/1‑101 (MQ=255)

             cGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCAcc                                                                                    >  1:6677862/1‑101 (MQ=255)

             cGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCAcc                                                                                    <  2:6677424/101‑1 (MQ=255)

              gATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCAGCTGCTGATCACCg                                                                                   >  2:6677863/1‑101 (MQ=255)

               aTGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGc                                                                                  >  1:6677867/1‑101 (MQ=255)

                tgttgtGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCa                                                                                 <  1:6677731/101‑1 (MQ=255)

                  ttgtGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAg                                                                               >  1:6677871/1‑101 (MQ=255)

                   tgtgCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCgct                                                                                              <  2:6677349/85‑1 (MQ=255)

                   tgtgCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGc                                                                              <  1:6677237/101‑1 (MQ=255)

                    gtgCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGcc                                                                             >  2:6677872/1‑101 (MQ=255)

                        cAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGAc                                                                         >  2:6677875/1‑101 (MQ=255)

                         aGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACa                                                                        >  1:6677880/1‑101 (MQ=255)

                           aTAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATa                                                                      <  1:6677322/101‑1 (MQ=255)

                            tAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATaa                                                                     >  2:6677883/1‑101 (MQ=255)

                               cGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCa                                                                  <  2:6677643/101‑1 (MQ=255)

                                ggTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATc                                                                   <  2:6677410/99‑1 (MQ=255)

                                ggTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAg                                                                 >  1:6677886/1‑101 (MQ=255)

                                 gTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGa                                                                >  2:6677888/1‑101 (MQ=255)

                                    aTGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACcg                                                             <  1:6677481/101‑1 (MQ=255)

                                     tGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACcga                                                            >  2:6677894/1‑101 (MQ=255)

                                      gTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACcgac                                                           >  1:6677896/1‑101 (MQ=255)

                                       tCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACcgacg                                                          >  2:6677899/1‑101 (MQ=255)

                                        ccGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACcgacga                                                         >  1:6677901/1‑101 (MQ=255)

                                         cGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGaa                                                        <  1:6677279/101‑1 (MQ=255)

                                         cGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGaa                                                        <  1:6677396/101‑1 (MQ=255)

                                          gCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGaaa                                                       >  2:6677907/1‑101 (MQ=255)

                                           ccATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGcc                                                                                                 >  1:6677909/1‑58 (MQ=255)

                                            cATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGAcc                                                              <  2:6677182/92‑1 (MQ=255)

                                             aTTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGc                                                                                  <  2:6677287/71‑1 (MQ=255)

                                             aTTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACg                                                    >  2:6677911/1‑101 (MQ=255)

                                              ttGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCgct                                                                                              >  1:6677915/1‑58 (MQ=255)

                                              ttGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACgg                                                   <  1:6677463/101‑1 (MQ=255)

                                                gtgATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGa                                                  >  2:6677918/1‑100 (MQ=255)

                                                gtgATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGAt                                                 <  1:6677184/101‑1 (MQ=255)

                                                 ggATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATg                                                >  1:6677921/2‑101 (MQ=255)

                                                   aTGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGtt                                              <  2:6677326/101‑1 (MQ=255)

                                                    tGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTg                                             <  1:6677300/101‑1 (MQ=255)

                                                     gCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGa                                            >  1:6677925/1‑101 (MQ=255)

                                                      ccAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAg                                           >  2:6677928/1‑101 (MQ=255)

                                                       cAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAgt                                          >  1:6677932/1‑101 (MQ=255)

                                                       cAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAgt                                          <  1:6677667/101‑1 (MQ=255)

                                                       cAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAgt                                          <  2:6677569/101‑1 (MQ=255)

                                                       cAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAgt                                          <  2:6677453/101‑1 (MQ=255)

                                                        aaGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAgtg                                         <  2:6677747/101‑1 (MQ=255)

                                                         aGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAgtgt                                        <  2:6677542/101‑1 (MQ=255)

                                                          gCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGtgtg                                       <  1:6677573/101‑1 (MQ=255)

                                                           cAACCCCATTGCGAACCAACTGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGc                                      <  2:6677257/101‑1 (MQ=255)

                                                           cAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGc                                      <  2:6677407/101‑1 (MQ=255)

                                                           cAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGCATGTTGAGTGTGc                                      <  1:6677548/101‑1 (MQ=255)

                                                            aaCCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCt                                     >  2:6677937/1‑101 (MQ=255)

                                                            aaCCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCt                                     <  1:6677097/101‑1 (MQ=255)

                                                            aaCCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCt                                     >  1:6677936/1‑101 (MQ=255)

                                                             aCCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTg                                    <  1:6677440/101‑1 (MQ=255)

                                                              ccccATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTg                                    <  2:6677587/100‑1 (MQ=255)

                                                               cccAATGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAAc                                                     <  2:6677444/82‑1 (MQ=255)

                                                                ccATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGcc                                 <  1:6677450/101‑1 (MQ=255)

                                                                 cATTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCa                                >  1:6677940/1‑101 (MQ=255)

                                                                  aTTGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACAGATGTTGAGTGTGCTGGCCAg                               >  2:6677943/1‑101 (MQ=255)

                                                                    tGCGAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCa                                <  2:6677510/98‑1 (MQ=255)

                                                                       gAACCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAAca                          >  1:6677947/1‑101 (MQ=255)

                                                                        aaCCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCt                                     >  2:6677953/1‑89 (MQ=255)

                                                                        aaCCAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAAcat                         >  2:6677950/1‑101 (MQ=255)

                                                                          ccAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAAcatca                       <  1:6677420/101‑1 (MQ=255)

                                                                           cAACAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAAcatcat                      >  2:6677957/1‑101 (MQ=255)

                                                                            aaCAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCAtt                     >  1:6677960/1‑101 (MQ=255)

                                                                              cAGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTAc                   >  2:6677963/1‑101 (MQ=255)

                                                                               aGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACg                  >  2:6677967/1‑101 (MQ=255)

                                                                               aGGCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACg                  >  1:6677965/1‑101 (MQ=255)

                                                                                 gCGGAAATAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGc                <  2:6677377/101‑1 (MQ=255)

                                                                                 gCGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGc                <  2:6677522/101‑1 (MQ=255)

                                                                                  cGGAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCa               >  2:6677972/1‑101 (MQ=255)

                                                                                   ggAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGt                                               >  1:6677974/1‑68 (MQ=255)

                                                                                   ggAAACAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCAc              <  1:6677457/101‑1 (MQ=255)

                                                                                     aaaCAACGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTa            >  2:6677977/1‑101 (MQ=255)

                                                                                      aacaacGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGAt                                                 <  1:6677342/63‑1 (MQ=255)

                                                                                      aacaacGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTaa           <  1:6677509/101‑1 (MQ=255)

                                                                                      aacaacGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTaa           >  1:6677980/1‑101 (MQ=255)

                                                                                        caacGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTaaga         <  2:6677363/101‑1 (MQ=255)

                                                                                        caacGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTaaga         >  2:6677982/1‑101 (MQ=255)

                                                                                        caacGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTaaga         <  2:6677536/101‑1 (MQ=255)

                                                                                         aacGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTaagaa        <  2:6677502/101‑1 (MQ=255)

                                                                                         aacGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTaagaa        <  1:6677529/101‑1 (MQ=255)

                                                                                         aacGCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTaagaa        >  2:6677985/1‑101 (MQ=255)

                                                                                            gCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGtt     >  1:6677991/1‑101 (MQ=255)

                                                                                            gCCCAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGtt     >  2:6677989/1‑101 (MQ=255)

                                                                                             cccAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTg    <  2:6677389/101‑1 (MQ=255)

                                                                                             cccAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTg    <  2:6677670/101‑1 (MQ=255)

                                                                                               cAGGCCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGcc  >  1:6677996/1‑101 (MQ=255)

ATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCC  >  NC_000913/4274809‑4275004

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