BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,809,113	G→A	pseudogene (798/1178 nt)	ygaY →	putative transporter YgaY, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,809,113	0	G	A	100.0%	829.4 / NA	203	pseudogene (798/1178 nt)	ygaY	putative transporter YgaY, N‑terminal fragment
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (88/115); total (88/115)

Only 100 of 217 total aligned reads displayed.

TTTTGCTTGCCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCGGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTCACACTTCCGT  >  NC_000913/2809015‑2809210

ttttGCTTGCCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGt                                                                                                 <  1:3630773/101‑1 (MQ=255)

ttttGCTTGCCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGt                                                                                                 <  2:3630955/101‑1 (MQ=255)

ttttGCTTGCCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGt                                                                                                 <  2:3630717/101‑1 (MQ=255)

 tttGCTTGCCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGtt                                                                                                >  1:3631296/1‑101 (MQ=255)

  ttgcttgcCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGttt                                                                                               >  1:3631299/1‑101 (MQ=255)

   tgcttgcCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGtttt                                                                                              <  1:3630980/101‑1 (MQ=255)

     cttgcCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGc                                                                                            <  1:3631020/101‑1 (MQ=255)

      ttgcCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGcc                                                                                           >  1:3631303/1‑101 (MQ=255)

        gcCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGa                                                                                         >  2:3631304/1‑101 (MQ=255)

        gcCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGa                                                                                         <  1:3631078/101‑1 (MQ=255)

         ccGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGAt                                                                                        >  2:3631306/1‑101 (MQ=255)

           gCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCg                                                                                          >  2:3631313/1‑97 (MQ=255)

           gCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATaa                                                                                      >  2:3631311/1‑101 (MQ=255)

           gCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATaa                                                                                      <  2:3630685/101‑1 (MQ=255)

           gCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATaa                                                                                      >  1:3631309/1‑101 (MQ=255)

            cTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAg                                                                                     <  1:3630710/101‑1 (MQ=255)

             tCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAgg                                                                                    >  2:3631316/1‑101 (MQ=255)

              ccACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAggg                                                                                   <  2:3631075/101‑1 (MQ=255)

              ccACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAggg                                                                                   <  1:3630957/101‑1 (MQ=255)

               cACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGc                                                                                  >  1:3631322/1‑101 (MQ=255)

               cACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGc                                                                                  >  1:3631319/1‑101 (MQ=255)

               cACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGc                                                                                  >  2:3631321/1‑101 (MQ=255)

                aCCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCa                                                                                 >  1:3631325/1‑101 (MQ=255)

                  cTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCaaa                                                                               <  2:3631143/101‑1 (MQ=255)

                  cTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCaaa                                                                               <  2:3631056/101‑1 (MQ=255)

                   ttttAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAAt                                                                              >  2:3631326/1‑101 (MQ=255)

                   ttttAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAAt                                                                              <  2:3631092/101‑1 (MQ=255)

                    tttAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATc                                                                             <  1:3631046/101‑1 (MQ=255)

                    tttAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATc                                                                             >  2:3631328/1‑101 (MQ=255)

                      tAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGc                                                                           <  2:3631030/101‑1 (MQ=255)

                      tAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGc                                                                           <  2:3630824/101‑1 (MQ=255)

                          tACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcacca                                                                       >  2:3631331/1‑101 (MQ=255)

                           aCAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcaccac                                                                      <  2:3630978/101‑1 (MQ=255)

                           aCAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcaccac                                                                      >  2:3631333/1‑101 (MQ=255)

                            cAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcaccaca                                                                     >  1:3631334/1‑101 (MQ=255)

                            cAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcaccaca                                                                     <  1:3631240/101‑1 (MQ=255)

                             aGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcaccacac                                                                    >  1:3631336/1‑101 (MQ=255)

                             aGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcaccacac                                                                    <  1:3630696/101‑1 (MQ=255)

                               cGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcaccacacca                                                                  >  1:3631339/1‑101 (MQ=255)

                               cGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcaccacacca                                                                  <  1:3630903/101‑1 (MQ=255)

                                gATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCaccacaccac                                                                 <  1:3630799/101‑1 (MQ=255)

                                gATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCaccacaccac                                                                 >  2:3631341/1‑101 (MQ=255)

                                 aTGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCAccacaccaca                                                                <  2:3630762/101‑1 (MQ=255)

                                   ggTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGGTAAGGGCAAATCGCACCACACCACAAc                                                              >  1:3631342/1‑101 (MQ=255)

                                    gtgtCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACt                                                             >  1:3631344/1‑101 (MQ=255)

                                     tgtCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACtt                                                            <  2:3631005/101‑1 (MQ=255)

                                     tgtCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACtt                                                            <  1:3630912/101‑1 (MQ=255)

                                      gtCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACttt                                                           <  1:3630931/101‑1 (MQ=255)

                                      gtCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACttt                                                           <  2:3630725/101‑1 (MQ=255)

                                        cATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCa                                                                                 >  1:3631349/1‑77 (MQ=255)

                                        cATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCg                                                         >  1:3631347/1‑101 (MQ=255)

                                         aTTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGcaccacacca                                                                  <  2:3631014/91‑1 (MQ=255)

                                          ttGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGt                                                       <  1:3630921/101‑1 (MQ=255)

                                          ttGGTCGGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGt                                                       >  1:3631352/1‑101 (MQ=255)

                                            ggTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTct                                                     <  1:3630995/101‑1 (MQ=255)

                                            ggTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTct                                                     <  1:3630987/101‑1 (MQ=255)

                                                tGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAggg                                                                                    >  2:3631354/1‑66 (MQ=255)

                                                tGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTctgctg                                                 <  1:3630740/101‑1 (MQ=255)

                                                     ggACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTgctgctgct                                            <  2:3631039/101‑1 (MQ=255)

                                                     ggACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTgctgctgct                                            >  2:3631356/1‑101 (MQ=255)

                                                      gACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCAccacaccaca                                                                >  2:3631358/1‑80 (MQ=255)

                                                       aCTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTAt                                          <  1:3631180/101‑1 (MQ=255)

                                                        cTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTc                                                      <  2:3630868/88‑1 (MQ=255)

                                                         ttGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTa                                        <  2:3630809/101‑1 (MQ=255)

                                                         ttGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTa                                        <  2:3630992/101‑1 (MQ=255)

                                                            cGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACtt                                     >  1:3631361/1‑101 (MQ=255)

                                                             gggAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACttt                                    <  1:3630928/101‑1 (MQ=255)

                                                             gggAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACttt                                    >  1:3631364/1‑101 (MQ=255)

                                                              ggAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTc                                   <  2:3631003/101‑1 (MQ=255)

                                                                aGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCAt                                 <  2:3631052/101‑1 (MQ=255)

                                                                aGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCAt                                 >  2:3631365/1‑101 (MQ=255)

                                                                  cTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCAtgg                               >  2:3631367/1‑101 (MQ=255)

                                                                  cTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCAtgg                               <  1:3630945/101‑1 (MQ=255)

                                                                  cTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCAtgg                               <  2:3631103/101‑1 (MQ=255)

                                                                    gCCGGAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCAtggct                             <  1:3630867/101‑1 (MQ=255)

                                                                       ggAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCAtggctggc                          >  2:3631370/1‑101 (MQ=255)

                                                                       ggAGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCAtggctggc                          >  1:3631371/1‑101 (MQ=255)

                                                                         aGCGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGa                        <  1:3631023/101‑1 (MQ=255)

                                                                           cGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTctgctgctg                                              <  1:3630937/77‑1 (MQ=255)

                                                                           cGTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATc                      <  2:3631202/101‑1 (MQ=255)

                                                                            gTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCt                     >  1:3631375/1‑101 (MQ=255)

                                                                            gTTGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCt                     >  2:3631376/1‑101 (MQ=255)

                                                                             ttGGGCGCTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTg                    <  2:3630926/101‑1 (MQ=255)

                                                                                 gcgcTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGttt                <  2:3630780/101‑1 (MQ=255)

                                                                                 gcgcTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGttt                >  1:3631378/1‑101 (MQ=255)

                                                                                   gcTCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTgg              <  2:3630848/101‑1 (MQ=255)

                                                                                     tCGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTc            <  1:3631036/101‑1 (MQ=255)

                                                                                      cGTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGtt                 <  1:3631008/95‑1 (MQ=255)

                                                                                       gTCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTcac          >  2:3631380/1‑101 (MQ=255)

                                                                                        tCCGGGGGGCAGTTTTGCCGATAAGGGCAAATCGcacca                                                                       >  2:3631385/1‑39 (MQ=37)

                                                                                        tCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACCCCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTcaca         <  2:3630951/101‑1 (MQ=255)

                                                                                        tCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCAtggct                             <  2:3630875/81‑1 (MQ=255)

                                                                                        tCCGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTcaca         >  2:3631382/1‑101 (MQ=255)

                                                                                         ccGGCGGGCAGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTCacac        >  1:3631386/1‑101 (MQ=255)

                                                                                           ggcgggcaGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTCACACtt      <  1:3631034/101‑1 (MQ=255)

                                                                                            gcgggcagTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTCACACTTc     >  1:3631391/1‑101 (MQ=255)

                                                                                            gcgggcagTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTCACACTTc     <  1:3631126/101‑1 (MQ=255)

                                                                                            gcgggcagTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTCACACTTc     >  2:3631389/1‑101 (MQ=255)

                                                                                               ggcagTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTCACACTTCCGt  <  2:3631179/101‑1 (MQ=255)

                                                                                               ggcagTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTCACACTTCCGt  <  2:3630962/101‑1 (MQ=255)

TTTTGCTTGCCGCTCCACCTTTTAACTACAGCGATGGTGTCATTGGTCTGTTTGGACTTGCGGGAGCTGCCGGAGCGTTGGGCGCTCGTCCGGCGGGCGGTTTTGCCGATAAGGGCAAATCGCACCACACCACAACTTTCGGTCTGCTGCTGCTATTACTTTCATGGCTGGCGATCTGGTTTGGTCACACTTCCGT  >  NC_000913/2809015‑2809210

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: `‑`