BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,950,224	A→G	I40V (ATC→GTC)	ilvE →	branched‑chain amino‑acid aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,950,224	0	A	G	100.0%	789.2 / NA	245	I40V (ATC→GTC)	ilvE	branched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (109/136); total (109/136)

Only 100 of 251 total aligned reads displayed.

TGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCGC  >  NC_000913/3950134‑3950313

 ggatCAATGGGCAGATGGGTCGCAAGGAAGACGCGAAAGTGCATGTGATGTCGCACGCGCTGCACAATGGCTCTTCGGTTTTTGAAGGCGTcc                                                                                        <  2:4769572/90‑1 (MQ=255)

 ggTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTcc                                                                                        <  2:6236564/93‑1 (MQ=255)

 ggTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTcc                                                                                        <  1:5098659/93‑1 (MQ=255)

 ggTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTcc                                                                                        >  2:173343/1‑93 (MQ=255)

  gTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCg                                                                                       <  1:5555908/93‑1 (MQ=255)

   ttCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGt                                                                                      <  2:1626691/93‑1 (MQ=255)

    tCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGtt                                                                                     <  2:2719458/93‑1 (MQ=255)

      aaTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGc                                                                                   <  2:34956/93‑1 (MQ=255)

      aaTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGc                                                                                   <  1:310054/93‑1 (MQ=255)

      aaTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGc                                                                                   >  2:2737626/1‑93 (MQ=255)

       aTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCt                                                                                  <  1:1853991/93‑1 (MQ=255)

       aTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCt                                                                                  <  1:2768357/93‑1 (MQ=255)

        tGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTa                                                                                 >  2:5103102/1‑93 (MQ=255)

         ggggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTAc                                                                                <  1:419496/93‑1 (MQ=255)

            gagaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGAc                                                                             <  2:4938730/93‑1 (MQ=255)

            gagaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGAc                                                                             >  2:1896926/1‑93 (MQ=255)

            gagaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATCGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGAc                                                                             >  1:5004342/1‑93 (MQ=255)

             agaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACt                                                                            <  1:4101280/93‑1 (MQ=255)

               aTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCg                                                                          >  2:1182936/1‑93 (MQ=255)

               aTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGCTTTTGAAGGCGTCCGTTGCTACGACTCg                                                                          >  2:3040093/1‑93 (MQ=255)

                tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGc                                                                         >  2:2382848/1‑93 (MQ=255)

                     cGCTGGGAAGACGCGAATGTGCATGTGATGTCGCACGCGCTGCTCTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACaaa                                                                    <  2:4079862/93‑1 (MQ=255)

                     cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACaaa                                                                    <  2:4374981/93‑1 (MQ=255)

                     cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACaaa                                                                    <  2:4427936/93‑1 (MQ=255)

                      gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAg                                                                   <  1:6055210/93‑1 (MQ=255)

                      gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAg                                                                   >  1:3728769/1‑93 (MQ=255)

                      gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAg                                                                   <  1:4809482/93‑1 (MQ=255)

                        tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGa                                                                 >  2:5498333/1‑93 (MQ=255)

                        tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGa                                                                 <  1:2514956/93‑1 (MQ=255)

                         gggAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGAc                                                                <  2:4600530/93‑1 (MQ=255)

                         gggAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGAc                                                                <  2:2556663/93‑1 (MQ=255)

                          ggAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGAcc                                                               <  1:116565/93‑1 (MQ=255)

                          ggAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGAcc                                                               >  2:730765/1‑93 (MQ=255)

                            aaGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCgg                                                             >  1:904239/1‑93 (MQ=255)

                             aGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGt                                                            <  2:5969825/93‑1 (MQ=255)

                             aGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGt                                                            >  1:5604334/1‑93 (MQ=255)

                               ccgcgAAGGTGCATGTGATGTCGCACGCGCTGCGCTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                                          <  2:197440/92‑1 (MQ=255)

                               aCGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                                          <  2:108732/93‑1 (MQ=255)

                                cgcgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGt                                                         >  1:1939079/1‑93 (MQ=255)

                                 gcgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTa                                                        <  1:3370492/93‑1 (MQ=255)

                                  ggAAGTTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAt                                                       <  2:3856396/92‑1 (MQ=255)

                                   gAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAtt                                                      >  1:1023940/1‑93 (MQ=255)

                                   gAAGGTGCATGTGATGTCGCACGAGATGCACTATGGCACTACGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAtt                                                      <  1:4955637/93‑1 (MQ=255)

                                    aaGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTc                                                     <  1:213254/93‑1 (MQ=255)

                                       gTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                                                  >  1:3342424/1‑93 (MQ=255)

                                       gTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                                                  >  2:2207183/1‑93 (MQ=255)

                                        ggCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGcc                                                 <  1:3832754/92‑1 (MQ=255)

                                         gCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCa                                                >  2:5896862/1‑93 (MQ=255)

                                           aTGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATc                                              <  2:74133/93‑1 (MQ=255)

                                           aTGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATc                                              <  2:2102304/93‑1 (MQ=255)

                                           aTGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATc                                              <  2:1472495/93‑1 (MQ=255)

                                              tgATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTg                                           >  1:4382835/1‑93 (MQ=255)

                                               gATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGa                                          <  2:5790555/93‑1 (MQ=255)

                                                 tGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGc                                        <  1:397377/93‑1 (MQ=255)

                                                 tGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGc                                        >  2:5400643/1‑93 (MQ=255)

                                                   tCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCat                                      >  1:2320814/1‑93 (MQ=255)

                                                    cGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCata                                     <  2:3136390/93‑1 (MQ=255)

                                                    cGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTACGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCata                                     >  2:5134627/1‑93 (MQ=255)

                                                     gCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCatat                                    >  2:4174419/1‑93 (MQ=255)

                                                     gCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCatat                                    <  1:467453/93‑1 (MQ=255)

                                                      cACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATg                                   <  1:4902540/93‑1 (MQ=255)

                                                       aCGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGc                                  <  1:2270160/93‑1 (MQ=255)

                                                        cgcgCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCa                                 <  1:3016003/93‑1 (MQ=255)

                                                        cgcgCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCa                                 <  1:5775689/93‑1 (MQ=255)

                                                        cgcgCTGCACTATGGCACGTCGGTTTTTGAAGGCGTCCTTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCa                                 <  1:130071/93‑1 (MQ=255)

                                                          cgcTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGc                               <  2:1793141/93‑1 (MQ=255)

                                                          cgcTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGc                               <  1:1697368/93‑1 (MQ=255)

                                                           gcTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg                              >  1:4419727/1‑93 (MQ=255)

                                                           gcTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg                              <  2:167747/93‑1 (MQ=255)

                                                            cTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGt                             <  2:5589163/93‑1 (MQ=255)

                                                             tGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTc                            >  2:3724202/1‑93 (MQ=255)

                                                              gCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCt                           >  1:2021857/1‑93 (MQ=255)

                                                               cACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTg                          <  1:269591/93‑1 (MQ=255)

                                                                aCTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTATGc                         >  1:2441068/1‑93 (MQ=255)

                                                                 cTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCa                        >  1:1520313/1‑93 (MQ=255)

                                                                  tATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCAt                       >  2:3183455/1‑93 (MQ=255)

                                                                  tATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCAt                       >  1:5281384/1‑93 (MQ=255)

                                                                   aTGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATg                      <  2:4284409/93‑1 (MQ=255)

                                                                    tGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGa                     <  1:3131514/93‑1 (MQ=255)

                                                                    tGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGa                     <  1:6242487/93‑1 (MQ=255)

                                                                    tGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGa                     <  2:4235566/93‑1 (MQ=255)

                                                                    tGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGa                     <  1:3698159/93‑1 (MQ=255)

                                                                     ggCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGAc                    >  1:4546177/1‑93 (MQ=255)

                                                                        aCTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTcc                 <  2:1496473/93‑1 (MQ=255)

                                                                        aCTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTcc                 <  2:4856423/93‑1 (MQ=255)

                                                                         cTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCg                >  2:3566509/1‑93 (MQ=255)

                                                                              gTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCaaa           >  1:630004/1‑93 (MQ=255)

                                                                               tttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCaaaa          >  1:4347342/1‑93 (MQ=255)

                                                                               tttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCaaaa          >  2:6212856/1‑93 (MQ=255)

                                                                               tttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCaaaa          >  1:3213566/1‑93 (MQ=255)

                                                                               tttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCaaaa          <  1:1865973/93‑1 (MQ=255)

                                                                                ttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAAt         >  2:2918914/1‑93 (MQ=255)

                                                                                ttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAAt         >  1:5385607/1‑93 (MQ=255)

                                                                                ttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAAt         <  1:5189646/93‑1 (MQ=255)

                                                                                   tGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTa      <  2:803360/93‑1 (MQ=255)

                                                                                      aGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCg   <  1:5244905/93‑1 (MQ=255)

                                                                                       ggCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCGc  <  1:854325/93‑1 (MQ=255)

                                                                                       ggCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCGc  >  1:590491/1‑93 (MQ=255)

                                                                                       ggCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCGc  <  1:4260108/93‑1 (MQ=255)

                                                                                       ggCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCGc  >  1:3800719/1‑93 (MQ=255)

TGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCGC  >  NC_000913/3950134‑3950313

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