BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,958,239	G→A	Q216* (CAG→TAG)	ptrA ←	protease 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,958,239	0	G	A	100.0%	1428.2 / NA	336	Q216* (CAG→TAG)	ptrA	protease 3
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (172/164); total (172/164)

Only 100 of 346 total aligned reads displayed.

CGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCTGACC  >  NC_000913/2958141‑2958336

  ccATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCtactgc                                                                                               <  1:4968697/101‑1 (MQ=255)

  ccATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCtactgc                                                                                               >  1:4969263/1‑101 (MQ=255)

   cATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCa                                                                                              >  2:4969266/1‑101 (MQ=255)

    aTTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCAc                                                                                             <  2:4968483/101‑1 (MQ=255)

      ttttGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCg                                                                                           <  1:4968673/101‑1 (MQ=255)

        ttGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGa                                                                                         <  1:4968817/101‑1 (MQ=255)

         tGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGAt                                                                                        <  2:4968769/101‑1 (MQ=255)

          gCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGAtt                                                                                       <  1:4968703/101‑1 (MQ=255)

          gCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGAtt                                                                                       >  1:4969273/1‑101 (MQ=255)

           ccAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTa                                                                                      >  1:4969276/1‑101 (MQ=255)

              aCTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCa                                                                                   >  2:4969278/1‑101 (MQ=255)

              aCTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCa                                                                                   <  2:4968811/101‑1 (MQ=255)

              aCTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCa                                                                                   <  1:4968912/101‑1 (MQ=255)

               cTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAg                                                                                  <  2:4968936/101‑1 (MQ=255)

                tCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAgg                                                                                 <  1:4968501/101‑1 (MQ=255)

                tCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAgg                                                                                 >  2:4969282/1‑101 (MQ=255)

                  cGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGtt                                                                               >  1:4969285/1‑101 (MQ=255)

                  cGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGtt                                                                               <  1:4968637/101‑1 (MQ=255)

                   ggCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGttt                                                                              <  1:4968832/101‑1 (MQ=255)

                      aGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTc                                                                           >  2:4969289/1‑101 (MQ=255)

                       gCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCg                                                                          <  2:4969044/101‑1 (MQ=255)

                         ggTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCt                                                                        <  2:4968794/101‑1 (MQ=255)

                         ggTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCt                                                                        >  2:4969292/1‑101 (MQ=255)

                           tttATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGt                                                                                <  2:4968674/91‑1 (MQ=255)

                            ttattaCTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTaa                                                                     <  1:4968874/101‑1 (MQ=255)

                            ttattaCTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTaa                                                                     >  2:4969296/1‑101 (MQ=255)

                             tattaCTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTaaa                                                                    <  1:4969154/101‑1 (MQ=255)

                             tattaCTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTaaa                                                                    <  1:4969101/101‑1 (MQ=255)

                             tattaCTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTaaa                                                                    >  1:4969298/1‑101 (MQ=255)

                               ttaCTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGt                                                                  <  1:4968848/101‑1 (MQ=255)

                               ttaCTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGt                                                                  >  1:4969303/1‑101 (MQ=255)

                               ttaCTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGt                                                                  <  2:4968997/101‑1 (MQ=255)

                               ttaCTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGt                                                                  >  2:4969304/1‑101 (MQ=255)

                                  cTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTc                                                               >  1:4969308/1‑101 (MQ=255)

                                  cTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTc                                                               <  1:4969075/101‑1 (MQ=255)

                                    gTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGa                                                             <  2:4969083/101‑1 (MQ=255)

                                    gTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGa                                                             <  2:4968758/101‑1 (MQ=255)

                                     tAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAg                                                            <  1:4968893/101‑1 (MQ=255)

                                     tAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAg                                                            >  2:4969312/1‑101 (MQ=255)

                                       aaTAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGt                                                          >  1:4969316/1‑101 (MQ=255)

                                        aTAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGtt                                                         >  1:4969320/1‑101 (MQ=255)

                                        aTAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGtt                                                         <  1:4968886/101‑1 (MQ=255)

                                         tAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTa                                                        >  2:4969322/1‑101 (MQ=255)

                                           aCCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTacc                                                      >  2:4969327/1‑101 (MQ=255)

                                            ccGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTacca                                                     >  2:4969330/1‑101 (MQ=255)

                                               ccTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTAccacca                                                  >  2:4969336/1‑101 (MQ=255)

                                               ccTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTAccacca                                                  >  1:4969335/1‑101 (MQ=255)

                                                  tcatcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGaa                                               >  2:4969341/1‑101 (MQ=255)

                                                  tcatcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGaa                                               <  2:4968721/101‑1 (MQ=255)

                                                   catcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGaaa                                              <  1:4969244/101‑1 (MQ=255)

                                                   catcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGaaa                                              >  2:4969345/1‑101 (MQ=255)

                                                    atcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGaaaa                                             <  2:4968956/101‑1 (MQ=255)

                                                    atcaAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGaaaa                                             >  1:4969347/1‑101 (MQ=255)

                                                       aaaTTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACtt                                          >  1:4969352/1‑101 (MQ=255)

                                                       aaaTTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACtt                                          >  2:4969353/1‑101 (MQ=255)

                                                        aaTTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACttt                                         >  1:4969356/1‑101 (MQ=255)

                                                        aaTTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACttt                                         <  2:4968872/101‑1 (MQ=255)

                                                         aTTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTg                                        <  1:4968805/101‑1 (MQ=255)

                                                         aTTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTg                                        >  2:4969359/1‑101 (MQ=255)

                                                           tGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGaa                                      <  2:4968627/101‑1 (MQ=255)

                                                             gCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAAcc                                    <  2:4968897/101‑1 (MQ=255)

                                                              cGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCg                                   >  2:4969366/1‑101 (MQ=255)

                                                                gAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCggg                                 <  2:4969135/101‑1 (MQ=255)

                                                                gAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCggg                                 >  2:4969369/1‑101 (MQ=255)

                                                                 aaTAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTCGAACCgggg                                >  1:4969370/1‑101 (MQ=255)

                                                                  aTAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGgt                               <  2:4968844/101‑1 (MQ=255)

                                                                  aTAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGgt                               >  2:4969374/1‑101 (MQ=255)

                                                                   tAGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCg                                   >  1:4969379/1‑96 (MQ=255)

                                                                    aGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGgtgt                             >  2:4969383/1‑101 (MQ=255)

                                                                    aGTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGgtgt                             >  1:4969384/1‑101 (MQ=255)

                                                                     gTACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGtgtg                            <  2:4968990/101‑1 (MQ=255)

                                                                      tACTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGc                           <  1:4968944/101‑1 (MQ=255)

                                                                       aCTTCTCGTGGGAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGcc                          >  1:4969388/1‑101 (MQ=255)

                                                                       aCTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGcc                          <  1:4968735/101‑1 (MQ=255)

                                                                        cTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGtt                                                                 <  2:4968569/61‑1 (MQ=255)

                                                                        cTTCTCGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCg                         <  2:4968975/101‑1 (MQ=255)

                                                                          tctcGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCggg                       >  1:4969395/1‑101 (MQ=255)

                                                                          tctcGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCggg                       >  2:4969394/1‑101 (MQ=255)

                                                                            tcGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGtt                     <  2:4968903/101‑1 (MQ=255)

                                                                            tcGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGtt                     <  1:4969003/101‑1 (MQ=255)

                                                                             cGTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGa                                                             >  2:4969398/1‑60 (MQ=255)

                                                                              gTGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTaa                   <  1:4968926/101‑1 (MQ=255)

                                                                               tGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAAt                  <  2:4968919/101‑1 (MQ=255)

                                                                               tGGAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAAt                  >  1:4969403/1‑101 (MQ=255)

                                                                                 gAAATCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATgg                <  1:4968428/101‑1 (MQ=255)

                                                                                  aaaTCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGt               <  2:4969239/101‑1 (MQ=255)

                                                                                  aaaTCTTTCAGCGCCTACTGCACCGGATTACCA‑GTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGtt              >  1:4969407/1‑101 (MQ=255)

                                                                                     tCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTc            >  1:4969412/1‑101 (MQ=255)

                                                                                     tCTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTc            >  2:4969413/1‑101 (MQ=255)

                                                                                      cTTTCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAgg                                                           >  1:4969416/1‑53 (MQ=255)

                                                                                       tttCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCCCCAGAAAACTTTGAACCGGGgtg                              >  2:4969420/1‑81 (MQ=255)

                                                                                         tCAGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTgcg        <  1:4968965/101‑1 (MQ=255)

                                                                                           aGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCt      <  1:4968781/101‑1 (MQ=255)

                                                                                           aGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCt      <  2:4968779/101‑1 (MQ=255)

                                                                                           aGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCt      >  2:4969427/1‑101 (MQ=255)

                                                                                           aGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCt      >  1:4969426/1‑101 (MQ=255)

                                                                                           aGCGCCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCt      <  1:4968948/101‑1 (MQ=255)

                                                                                            gcgcCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCt           >  2:4969431/1‑95 (MQ=255)

                                                                                             cgcCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCTGa    <  1:4968575/101‑1 (MQ=255)

                                                                                              gcCTACTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCTGAc   >  1:4969434/1‑101 (MQ=255)

CGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCTTAAAGTTTCGAGGTTACCACCAGAAAACTTTGAACCGGGGTGTGCCGGGTTAATGGTTTCTGCGCTGACC  >  NC_000913/2958141‑2958336

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