BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,806,842	C→T	R226C (CGC→TGC)	kdtA →	3‑deoxy‑D‑manno‑octulosonic‑acid transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,806,842	0	C	T	100.0%	1571.7 / NA	427	R226C (CGC→TGC)	kdtA	3‑deoxy‑D‑manno‑octulosonic‑acid transferase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (182/245); total (182/245)

Only 100 of 464 total aligned reads displayed.

GCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTG  >  NC_000913/3806752‑3806931

gcAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACgctgcgctgc                                                                                         <  1:6484533/93‑1 (MQ=255)

 cAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCgctgcc                                                                                        <  1:4719415/93‑1 (MQ=255)

 cAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCgctgcc                                                                                        >  1:253079/1‑93 (MQ=255)

  aaaaaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCa                                                                                       <  2:8220521/93‑1 (MQ=255)

   aaaaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAg                                                                                      <  2:4109975/93‑1 (MQ=255)

     aaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTg                                                                                    >  1:8531772/1‑93 (MQ=255)

       ataatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTggg                                                                                  <  2:9394420/93‑1 (MQ=255)

       ataatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTggg                                                                                  <  1:4882945/93‑1 (MQ=255)

         aatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGca                                                                                <  1:2706867/93‑1 (MQ=255)

          atCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcac                                                                               >  1:9736165/1‑93 (MQ=255)

          atCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcac                                                                               <  1:3817472/93‑1 (MQ=255)

            cAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcacca                                                                             <  2:562968/93‑1 (MQ=255)

             aGGTGCCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTCCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcaccac                                                                            <  1:4541492/93‑1 (MQ=255)

              ggTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAGGCCGTGACGCTGCGCTGCCAGTGGGCACcaca                                                                           <  2:5119006/93‑1 (MQ=255)

               gTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCacac                                                                          <  1:4013236/93‑1 (MQ=255)

               gTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCacac                                                                          >  1:9256845/1‑93 (MQ=255)

               gTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCacac                                                                          >  1:9423603/1‑93 (MQ=255)

               gTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCacac                                                                          >  1:4623613/1‑93 (MQ=255)

                 gACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCg                                                                        >  2:9654943/1‑93 (MQ=255)

                  aCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCAACACACCGc                                                                       >  2:3377365/1‑93 (MQ=255)

                     gTTACCGGTAGCCTGAATTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCg                                                                    >  2:8981947/1‑93 (MQ=255)

                     gTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCg                                                                    >  2:6532769/1‑93 (MQ=255)

                      ttACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCgg                                                                   >  1:8273920/1‑93 (MQ=255)

                       tACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTTCGCTGCCAGTGGGCACCACACCCCCCGGt                                                                  <  2:9880396/93‑1 (MQ=255)

                        aCCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTa                                                                 <  1:7631164/93‑1 (MQ=255)

                         ccGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTag                                                                >  2:6782006/1‑92 (MQ=255)

                         ccGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTAt                                                                <  2:3550298/93‑1 (MQ=255)

                         ccGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTAt                                                                <  1:7333256/93‑1 (MQ=255)

                           ggTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATgg                                                              >  2:274656/1‑93 (MQ=255)

                           ggTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATgg                                                              >  2:8764760/1‑93 (MQ=255)

                            gTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCTCTGCCAGGGGGCACCACCCCGCCCGGTATggt                                                             >  2:7506770/1‑92 (MQ=255)

                            gTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGa                                                             >  1:1546453/1‑93 (MQ=255)

                            gTAGCCTAAAATTCGATATTTCTGTCACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGa                                                             <  2:3838846/93‑1 (MQ=255)

                             tAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGAt                                                            >  2:8288083/1‑93 (MQ=255)

                             tAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGAt                                                            <  1:80879/93‑1 (MQ=255)

                             tAGCCTAAATGTCGATATTTCTGTAACGCCGCAGTGGGTTTCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGAt                                                            <  2:9195193/93‑1 (MQ=255)

                              aGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGAtt                                                           <  1:9963262/93‑1 (MQ=255)

                               gCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTg                                                          >  2:2003981/1‑93 (MQ=255)

                               gCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTg                                                          <  1:2936154/93‑1 (MQ=255)

                                  tGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGcca                                                       <  2:7196418/93‑1 (MQ=255)

                                   gAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATCGATTGtcac                                                      >  2:6302236/1‑93 (MQ=255)

                                   gAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGATGCCAGTGGGCACCACACCGCCCGGTATGGATTGccac                                                      <  2:7914367/93‑1 (MQ=255)

                                    aaaTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccacc                                                     <  2:1276149/93‑1 (MQ=255)

                                     aaTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccacca                                                    <  1:3495911/93‑1 (MQ=255)

                                     aaTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccacca                                                    >  1:6294249/1‑93 (MQ=255)

                                     aaTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGAACCACACCGCCCGGCATGGATTGccacca                                                    >  2:2551603/1‑93 (MQ=255)

                                     aaTTCGATATTTCTGGAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccacca                                                    <  2:7595396/93‑1 (MQ=255)

                                       ttCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGc                                                  <  2:5392012/93‑1 (MQ=255)

                                         cGATATTTCTGTACCGCCAAAGATGGCTGCTAACGCCGTGACGCTGCGCTGCCAGTGGGCCCCACACCACACGATATGGATTGCCACCAGCAc                                                <  1:8951635/93‑1 (MQ=255)

                                         cGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCAc                                                <  1:6716033/93‑1 (MQ=255)

                                         cGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCAc                                                <  2:3973012/93‑1 (MQ=255)

                                         cGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCAc                                                <  2:8008644/93‑1 (MQ=255)

                                          gATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCAcc                                               >  1:1810750/1‑92 (MQ=255)

                                          gATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACt                                               <  2:4922843/93‑1 (MQ=255)

                                           atatTTATGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTc                                              >  2:2409064/1‑93 (MQ=255)

                                            tatTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGGCGCTGCGCTGCCAGTGGGCACCACGCCGCCCGGTATGGATTGCCACCAGCACTCa                                             <  2:3038054/93‑1 (MQ=255)

                                              tttCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACg                                           >  1:7481933/1‑93 (MQ=255)

                                               ttCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGa                                          <  2:8521787/93‑1 (MQ=255)

                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCCGCACCCACGaa                                         >  1:6928933/1‑93 (MQ=255)

                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                         >  2:4522921/1‑93 (MQ=255)

                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                         <  2:1469735/93‑1 (MQ=255)

                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                         <  1:9571388/93‑1 (MQ=255)

                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                         <  1:7203401/93‑1 (MQ=255)

                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                         <  1:1595200/93‑1 (MQ=255)

                                                 cTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAg                                        >  1:9072576/1‑93 (MQ=255)

                                                  tGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAgg                                       >  1:7936363/1‑93 (MQ=255)

                                                   gTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGc                                      <  2:7746059/93‑1 (MQ=255)

                                                   gTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGc                                      >  2:5948395/1‑93 (MQ=255)

                                                      aCGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGaa                                   <  2:4749824/93‑1 (MQ=255)

                                                        gccgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAaga                                 >  1:3264356/1‑93 (MQ=255)

                                                          cgcAGTTGGCTGCGAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagaga                               <  1:418685/93‑1 (MQ=255)

                                                           gcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagagag                              <  2:5854805/93‑1 (MQ=255)

                                                            cAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgt                             <  2:1886666/93‑1 (MQ=255)

                                                            cAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgt                             <  2:9568298/93‑1 (MQ=255)

                                                                tGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTgtgg                         >  1:1013958/1‑93 (MQ=255)

                                                                tGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTgtgg                         >  2:2181182/1‑93 (MQ=255)

                                                                  gctgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTgtggtg                       >  2:1796539/1‑93 (MQ=255)

                                                                   ctgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGa                      >  1:5819289/1‑93 (MQ=255)

                                                                    tgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGAt                     >  1:7096187/1‑93 (MQ=255)

                                                                     gctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATc                    <  1:1256942/93‑1 (MQ=255)

                                                                      ctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcg                   >  2:6959136/1‑93 (MQ=255)

                                                                       tAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgc                  >  2:4253869/1‑93 (MQ=255)

                                                                       tAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgc                  >  1:2201764/1‑93 (MQ=255)

                                                                       tAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgc                  <  2:5328658/93‑1 (MQ=255)

                                                                            ccGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGcac             <  2:6386120/93‑1 (MQ=255)

                                                                              gTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACAt           <  2:2805650/93‑1 (MQ=255)

                                                                              gTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACAt           <  1:6036964/93‑1 (MQ=255)

                                                                              gTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACAt           >  2:7428086/1‑93 (MQ=255)

                                                                               tGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc          >  1:5469354/1‑93 (MQ=255)

                                                                               tGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc          <  1:3897701/93‑1 (MQ=255)

                                                                               tGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc          >  1:2017765/1‑93 (MQ=255)

                                                                                gACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCa         >  1:5620943/1‑93 (MQ=255)

                                                                                gACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCa         <  2:4424576/93‑1 (MQ=255)

                                                                                  cgctgcgctgCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAgg       <  1:8886263/93‑1 (MQ=255)

                                                                                   gctgcgctgcCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGGGTGGTGATCGCCGCCCATCAgga      >  2:8642059/1‑92 (MQ=255)

                                                                                     tGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCAt    <  1:8350972/93‑1 (MQ=255)

                                                                                     tGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCAt    <  2:1577226/93‑1 (MQ=255)

                                                                                     tGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCAt    >  1:6577872/1‑93 (MQ=255)

                                                                                      gcgcTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCAtt   <  1:5361456/93‑1 (MQ=255)

                                                                                       cgctgcCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCCCCGCACACCAGGCAttt  >  1:3120441/1‑92 (MQ=255)

GCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTG  >  NC_000913/3806752‑3806931

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