BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,582,538	A→G	G336G (GGT→GGC)	yffG ←	putative oxidoreductase Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,582,538	0	A	G	100.0%	429.1 / NA	140	G336G (GGT→GGC)	yffG	putative oxidoreductase Fe‑S subunit
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (54/86); total (54/86)

Only 100 of 145 total aligned reads displayed.

GAGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGACCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGAAAT  >  NC_000913/2582448‑2582627

gAGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGcc                                                                                         <  1:4012052/93‑1 (MQ=255)

gAGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGcc                                                                                         <  2:3850477/93‑1 (MQ=255)

  gCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTg                                                                                       <  2:4539591/93‑1 (MQ=255)

  gCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTg                                                                                       >  2:3531514/1‑93 (MQ=255)

   cAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGc                                                                                      >  1:1246682/1‑93 (MQ=255)

   cAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGc                                                                                      >  1:444989/1‑93 (MQ=255)

   cAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGc                                                                                      <  1:2430135/93‑1 (MQ=255)

     aGCCACCGATTTCTGGATGGCGCTCGTACACCGTCCCCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCAc                                                                                    <  1:840988/93‑1 (MQ=255)

     aGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCAc                                                                                    >  1:4595565/1‑93 (MQ=255)

     aGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCAc                                                                                    >  1:1936640/1‑93 (MQ=255)

      gCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGTCTGCAcc                                                                                   >  2:398151/1‑93 (MQ=255)

       ccACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCg                                                                                  >  1:3212424/1‑93 (MQ=255)

        cACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGa                                                                                 >  2:630597/1‑93 (MQ=255)

        cACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGa                                                                                 <  1:3982210/93‑1 (MQ=255)

         aCCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGAt                                                                                >  1:1363152/1‑93 (MQ=255)

         aCCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGAt                                                                                >  2:2206214/1‑93 (MQ=255)

            gATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAAt                                                                             >  1:2271861/1‑93 (MQ=255)

            gATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAAt                                                                             <  1:3766119/93‑1 (MQ=255)

             aTTTCTGGATGGCGATCGTGCAGGTTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGAGAATc                                                                            <  1:1531820/93‑1 (MQ=255)

             aTTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATc                                                                            >  2:1746909/1‑93 (MQ=255)

              tttCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCg                                                                           >  2:411564/1‑93 (MQ=255)

                    gATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCAccc                                                                     >  1:4322679/1‑93 (MQ=255)

                        gCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCtt                                                                 <  2:450059/93‑1 (MQ=255)

                          gATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGt                                                               <  1:2038750/93‑1 (MQ=255)

                          gATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCCGGCCTGCACCGATAATCGCCACCCGCTTGt                                                               >  2:4248011/1‑93 (MQ=255)

                           aTCGTACACCGTCTCCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTc                                                              <  2:1297148/93‑1 (MQ=255)

                           aTCGTACACCGTCACCCCCACGCCATTGCGGGTCAGGACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACAGATAGTCGCCACACGCGTGTc                                                              >  2:1064026/1‑93 (MQ=21)

                           aTCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTc                                                              >  1:4581663/1‑93 (MQ=255)

                           aTCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTc                                                              >  1:2402942/1‑93 (MQ=255)

                           aTCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCGTGTc                                                              >  2:1954822/1‑93 (MQ=255)

                            tCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTcc                                                             >  1:1659056/1‑93 (MQ=255)

                            tCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTcc                                                             <  1:2300296/93‑1 (MQ=255)

                               tACACCGCCACCCCCCCGCCATTGCGGGTCAGACCATCCGCACAGGCCAGACCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACt                                                          <  1:3308666/93‑1 (MQ=25)

                                 cacCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACttt                                                        <  1:4225433/93‑1 (MQ=255)

                                 cacCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACttt                                                        >  1:2773814/1‑93 (MQ=255)

                                  cccGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTg                                                       <  2:1888556/92‑1 (MQ=255)

                                  acCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTg                                                       >  2:1472432/1‑93 (MQ=255)

                                   ccGTCCCCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTgg                                                      <  2:3473660/93‑1 (MQ=255)

                                   ccGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTgg                                                      >  1:319841/1‑93 (MQ=255)

                                   ccGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTgg                                                      <  2:2481146/93‑1 (MQ=255)

                                     gTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTg                                                    >  1:2910963/1‑93 (MQ=255)

                                      tCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGa                                                   >  2:3252042/1‑93 (MQ=255)

                                      tCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGa                                                   <  1:1482105/93‑1 (MQ=255)

                                      ccccccccACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGa                                                   <  2:3834896/90‑1 (MQ=255)

                                        ccccccACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACa                                                 <  2:197048/92‑1 (MQ=255)

                                        ccccccACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACa                                                 <  1:3829774/92‑1 (MQ=255)

                                        aCCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTGCACTTTGGTGACa                                                 >  2:952039/1‑93 (MQ=255)

                                              aCGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCtt                                           <  2:2944094/93‑1 (MQ=255)

                                               cGCCATTGGGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATAGTCACCCGCGTGTTCACCTTGGTGACATTGCTTa                                          >  2:1726941/1‑93 (MQ=255)

                                               cGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTa                                          >  1:1194689/1‑93 (MQ=255)

                                                gCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTaa                                         <  1:4215791/93‑1 (MQ=255)

                                                 ccATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAg                                        >  1:2604540/1‑93 (MQ=255)

                                                   aTTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTc                                      >  1:4631851/1‑93 (MQ=255)

                                                   aTTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTc                                      <  1:3119930/93‑1 (MQ=255)

                                                    ttGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCa                                     <  2:2058908/93‑1 (MQ=255)

                                                     tGCGGGACAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAg                                    <  1:1256310/93‑1 (MQ=255)

                                                      gCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAgg                                   <  1:290680/93‑1 (MQ=255)

                                                       cGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGa                                  <  1:3722498/93‑1 (MQ=255)

                                                        gggTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGAc                                 <  1:1677496/93‑1 (MQ=255)

                                                             aGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTGAGGGCGCCaa                            >  1:1672201/1‑93 (MQ=255)

                                                             aGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCaa                            <  2:3607012/93‑1 (MQ=255)

                                                             aGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCaa                            <  2:3210052/93‑1 (MQ=255)

                                                              gAACATCCGCACAGGCCCGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAc                           <  2:235449/93‑1 (MQ=255)

                                                              gAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAc                           <  2:2815272/93‑1 (MQ=255)

                                                              gAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAc                           >  2:2148369/1‑93 (MQ=255)

                                                              gAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAAGCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAc                           >  1:2518029/1‑93 (MQ=255)

                                                                aCATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCt                         <  2:4140966/93‑1 (MQ=255)

                                                                 cATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCtt                        >  1:3699274/1‑93 (MQ=255)

                                                                  aTCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCttt                       >  1:3556390/1‑93 (MQ=255)

                                                                   tCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTc                      <  1:3491182/93‑1 (MQ=255)

                                                                   tCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTc                      <  1:3852715/93‑1 (MQ=255)

                                                                    ccGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCg                     <  2:2141942/93‑1 (MQ=255)

                                                                     cGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGc                    <  2:101824/93‑1 (MQ=255)

                                                                      gCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGcc                   >  1:3284420/1‑93 (MQ=255)

                                                                       cacaGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGGGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCa                  >  1:3441981/1‑93 (MQ=255)

                                                                        acaGGCCACCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCaa                 <  2:1786140/93‑1 (MQ=255)

                                                                         caGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAAc                <  1:2326131/93‑1 (MQ=255)

                                                                         caGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAAc                <  1:3151558/93‑1 (MQ=255)

                                                                            gCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGcc             >  1:3599500/1‑93 (MQ=255)

                                                                             ccAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCt            >  2:3280772/1‑93 (MQ=255)

                                                                               aGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGa          <  1:4532031/93‑1 (MQ=255)

                                                                               aGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGa          <  2:4306177/93‑1 (MQ=255)

                                                                               aGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGa          <  2:4382748/93‑1 (MQ=255)

                                                                                gCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGAt         <  1:4664612/93‑1 (MQ=255)

                                                                                gCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGAt         <  1:4289356/93‑1 (MQ=255)

                                                                                 cccTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATc        >  2:4160097/1‑93 (MQ=255)

                                                                                 cccTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATc        <  1:3099711/93‑1 (MQ=255)

                                                                                 cccTGCCCGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATc        <  2:3103553/93‑1 (MQ=255)

                                                                                  ccTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCt       <  1:2183015/93‑1 (MQ=255)

                                                                                  ccTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCt       <  1:1625815/93‑1 (MQ=255)

                                                                                   cTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTg      <  1:3064305/93‑1 (MQ=255)

                                                                                   cTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTg      >  2:888598/1‑93 (MQ=255)

                                                                                    tGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGa     <  1:116585/93‑1 (MQ=255)

                                                                                     gCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGaa    >  1:755622/1‑93 (MQ=255)

                                                                                      ccGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGaaa   <  1:1077103/93‑1 (MQ=255)

                                                                                      ccGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGaaa   <  2:501914/93‑1 (MQ=255)

                                                                                      ccGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGaaa   <  1:4380790/93‑1 (MQ=255)

                                                                                      ccGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGaaa   <  2:1191939/93‑1 (MQ=255)

                                                                                      ccGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGaaa   <  1:3081409/93‑1 (MQ=255)

                                                                                       cGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGAAAt  <  2:2715911/93‑1 (MQ=255)

GAGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGACCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGAAAT  >  NC_000913/2582448‑2582627

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